poj 3254(状态压缩dp专题)

点击打开链接


题意:

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


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

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

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

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


#include"stdio.h"
#include"string.h"
#define mod 100000000
int 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写成n
		if(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;
}


你可能感兴趣的:(dp,状态压缩dp)