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;
}



posted @ 2012-02-18 21:08  Because Of You  Views(283)  Comments(0Edit  收藏  举报