POJ 3254 状压DP

一个n*m的矩阵，每个格子是0或者1，1表示土壤肥沃可以种植草地，0则不可以。在种草地的格子可以放牛，但边相邻的两个格子不允许同时放牛，问总共有多少种放牛的方法？（不放牛也算一种情况）

状压DP ，POJ1185 HDU4539升级版

#include "stdio.h"
#include "string.h"
int b[]={1,2,4,8,16,32,64,128,256,512,1024,2048,4096,8192};
int mod=100000000;
int a[13][13],sum[13][430],dp[13][430],mark[430];

int main()
{
    int ans,n,m,cnt,aim,i,j,k;

    while (scanf("%d%d",&n,&m)!=EOF)
    {
        for (i=0;i<n;i++)
            for (j=0;j<m;j++)
            scanf("%d",&a[i][j]);

        cnt=0;
        aim=1<<m;
        for (i=0;i<aim;i++)
            if (i&(i<<1));else  mark[cnt++]=i; // 可行的方案数



        memset(sum,0,sizeof(sum));
        for (i=0;i<n;i++)
            for (j=0;j<cnt;j++)
                for (k=0;k<m;k++)
                    if (a[i][k]==1 && (b[k]&mark[j])==b[k])
                        sum[i][j]+=b[k]; // sum[i][j]==mark[j]  表示第i行的j方案是可行的

        memset(dp,0,sizeof(dp));
       for (i=0;i<cnt;i++)
        if (sum[0][i]==mark[i]) dp[0][i]=1;

        for (i=1;i<n;i++)
            for (j=0;j<cnt;j++)
                for (k=0;k<cnt;k++)
                {
                    if (mark[k]&mark[j]) continue; // 和上个方案没有重复点
                    if (dp[i-1][k]==0) continue; // 上个方案可行
                    if (sum[i][j]!=mark[j]) continue; //本方案可行
                    dp[i][j]=(dp[i][j]+dp[i-1][k])%mod; //累加
                }
        ans=0;
        for (i=0;i<cnt;i++)
            ans=(ans+dp[n-1][i])%mod;
        printf("%d\n",ans);


    }
    return 0;
}