【状压DP】BZOJ 1725 Usaco Corn Fields牧场的安排

#include <cstdio>
#include <algorithm>

const int State = (1<<12)+5;
const int MOD = 100000000;

int n,m,Count;
int dp[20][State],Map[20];
int num[State],now[State];
bool Can[State][State];
long long ans;


//预处理 
void Run(){
    for(int i=0;i<(1<<m);i++){
        if(i & (i>>1) ) continue;
        now[++Count] = i;
        for(int j=i;j;j>>=1) num[Count]+=(j&1);//当前状态会有几个地方有草 
    }
    for(int i=0;i<=Count;i++)
        for(int j=0;j<=Count;j++)
            if( !(now[i]&now[j]) ) Can[i][j]=1;//这两种方案是否冲突 
}

int main(){
    scanf("%d%d",&n,&m);
    Run();
    for(int i=1;i<=n;i++){
        int x ;
        for(int j=0;j<m;j++){
            scanf("%d",&x);
            Map[i]|=(x<<j);
            //将本行的状态压缩  
        }
    }
    for(int i=1;i<=Count;i++) if( now[i] == (now[i]&Map[1]) ) dp[1][i]=1;
    for(int i=2;i<=n;i++){
        for(int j=1;j<=Count;j++){
            if( (now[j]&Map[i]) != now[j] ) continue;
            for(int k=1;k<=Count;k++){
                if( Can[j][k] ) dp[i][j] = (dp[i][j]+dp[i-1][k])%MOD;
            }
        }
    }
    for(int i=1;i<=Count;i++) ans=(ans+dp[n][i])%MOD;
    printf("%lld\n",ans);
    return 0;
}