poj 3254 状态压缩DP

状态转移：

dp[state][i]=sigma(dp[state'][i-1]);state和state'都是合法状态

更多的注释在代码中

View Code

#include<cstdio>
#include<cstring>
int dp[13][1<<13];
int cur[15];
int state[1<<13];
int m,n,tot;
const int mod = 100000000;
bool ok(int state)
{
    if(state&(state<<1)) return false;
    return true;
}
void init()
{
    for(int i=0;i<(1<<n);i++)
    {
        if(ok(i))  state[++tot]=i;
    }
}
bool fit(int state,int id)//判断某种没有两个相邻的状态与地图是否适合
{
    if(state&(cur[id])) return false;
    return true;
}
int main()
{
    int i,j,num;
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        tot=0;
        init();    
        for(i=1;i<=m;i++)
        {
            cur[i]=0;
            for(j=1;j<=n;j++)
            {
                scanf("%d",&num);
                if(num==0) cur[i]+=(1<<(n-j));//给地图为0的地方都标记为1
            }
        }
        memset(dp,0,sizeof(dp));
        for(i=1;i<=tot;i++)
        {
            if(fit(state[i],1))
                dp[1][i]=1;
        }
        for(i=2;i<=m;i++)
        {
            for(j=1;j<=tot;j++)
            {
              if(!fit(state[j],i)) continue;
              for(int k=1;k<=tot;k++)
              {
                  if(!fit(state[k],i-1)) continue;
                  if(state[k]&state[j]) continue;//上下两行不能有相邻的
                  dp[i][j]=(dp[i][j]+dp[i-1][k])%mod;
              }
            }
        }
        int ans=0;
        for(i=1;i<=tot;i++)
            ans=(ans+dp[m][i])%mod;
        printf("%d\n",ans);
    }
    return 0;
}