P1879 [USACO06NOV]Corn Fields G（状压dp）

题目传送门

题意： 有一张大小为N*M的由0/1组成的地图，你可以在标有1的地方种草，但是你种的任意两块草不能有公共边，问你有多少种不同的方案（不种也是一种方案）？（N<=12,M<=12）

思路： 看到这个数据范围我们很容易想到用二进制去压缩每一行，在同一行合法的情况去考虑每一行同时合法。我们先预处理出一行内所有可能的情况，然后根据输入的地图，用二进制将0标出，用于判断一个情况是否在该行合法。然后枚举所有情况，f[i][j]表示前i行，第i行是状态j时的方案数，同时需要枚举一个状态k，表示第i-1行的状态，在合法情况下由f[i-1][k]转移到f[i][j]即可。

代码：

#include<bits/stdc++.h>
#define endl '\n'
#define mp make_pair
#define pb push_back
#define ll long long
#define int long long
#define pii pair<int,int>
#define sz(x) (int)(x).size()
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read()
{
    int x=0,f=1;
    char ch=gc();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            f=-1;
        ch=gc();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=gc();
    }
    return x*f;
}
using namespace std;
const int N=1e3+10;
const int inf=0x3f3f3f3f;
const int mod=100000000;

int sta[N],tot=0,n,m,line[N],a[N],f[N][N];
void dfs(int cur,int pos,int x)
{
    if(pos > m)
    {
        sta[++tot] = cur;
        return ;
    }
    dfs(cur,pos+1,x);
    dfs(cur|(1<<pos),pos+2,x+1);
}

int solve()
{
    cin>>n>>m;
    dfs(0,1,0);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            cin>>a[j];
            if(a[j]==0)
                line[i] |= (1<<j);
        }
    }
    for(int i=1;i<=tot;i++)
    {
        if(sta[i]&line[1])
            continue;
        f[1][i] = 1;
    }
    for(int i=2;i<=n;i++)
    {
        for(int j=1;j<=tot;j++)
        {
            if(sta[j]&line[i])
                continue;
            for(int k=1;k<=tot;k++)
            {
                if(sta[k]&line[i-1])
                    continue;
                if(sta[k]&sta[j])
                    continue;
                f[i][j] = (f[i][j]+f[i-1][k])%mod;
            }
        }
    }
    int ans = 0;
    for(int i=1;i<=tot;i++)
        ans = (ans + f[n][i])%mod;
    cout<<ans<<endl;


    return 0;
}
signed main()
{
    solve();
    return 0;
}