poj 3254（状态压缩dp专题）

来源：互联网发布：银行家算法实验报告编辑：程序博客网时间：2024/05/29 02:48

点击打开链接

题意：

给你一个n*m的矩阵，由0或1表示，要在1的位置放牛，问有多少种放牛的方式。

人生第一道状态压缩dp。。。。

state，求出所有可能出现的状态

cur【i】，求出第i行0所表示的二进制数。。

fit（x,k）函数表示第x行是否满足第k种状态。。。。

#include"stdio.h"#include"string.h"#define mod 100000000int n,m;int top;int state[1<<13];int cur[13];int ok(int i){if(i&(i<<1))return 0;return 1;}void init(){int i;top=0;for(i=0;i<(1<<m);i++)//刚开始一直把m写成nif(ok(i))state[++top]=i;}int fit(int x,int k){if(cur[x]&state[k])return 0;return 1;}int main(){int i,j,k;int dp[13][1<<13];while(scanf("%d%d",&n,&m)!=-1){init();for(i=1;i<=n;i++){cur[i]=0;for(j=1;j<=m;j++){scanf("%d",&k);if(k==0)cur[i]+=(1<<(m-j));}}memset(dp,0,sizeof(dp));for(i=1;i<=top;i++){if(fit(1,i))dp[1][i]=1;}for(i=2;i<=n;i++){for(j=1;j<=top;j++){if(!fit(i,j))continue;for(k=1;k<=top;k++){if(!fit(i-1,k))continue;//第i-1行不适合第k中状态if(state[j]&state[k])continue;//状态j和k同一列相邻dp[i][j]=(dp[i][j]+dp[i-1][k])%mod;}}}int ans=0;for(i=1;i<=top;i++)ans=(ans+dp[n][i])%mod;printf("%d\n",ans);}return 0;}