POJ 2411 Mondriaan's Dream

最新推荐文章于 2021-03-30 20:02:40 发布

原创最新推荐文章于 2021-03-30 20:02:40 发布 · 780 阅读

1 ·

CC 4.0 BY-SA版权

DP 专栏收录该内容

39 篇文章

订阅专栏

http://poj.org/problem?id=2411

比较简单经典的状态压缩DP，题意是一个N*M(n,m<=11)的方格组成的块,要求用1*2的块去填充,问有填充满，有多少种方案。

dp[i][j]， i 表示到摆到第i行,此行的状态为j的方案数，j的状态是此行全部填充完的情况下，突出到下一行的状态，

为了方便，让dp[1][]表示摆放第一行之前的状态，最后求的就是dp[n+1][0]。

dp[i][j]=dp[i][j]+dp[i-1][k], （如果j可以由k转移过来，就是i ，j没有任何相同突出的列，而且其他空出的格子每块都为偶数个）.

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
#include<set>
using namespace std;
#define INF 1000000000
typedef __int64 LL;
#define N 13
int n,m;
LL dp[N][2050];
int test(int x,int y)
{
	if(x&y) return 0;
	int now=x|y,last=0;
	for(int i=0;i<m;i++)
	{
		if((now>>i)&1)
		{
			if(last&1) return 0;
			last=0;
		}
		else last++;
	}
	if(last&1) return 0;
	return 1;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen ("in.txt" , "r" , stdin);
#endif
	while(scanf("%d%d",&n,&m),n||m)
	{
		if((n&1) && (m&1))
		{
			printf("0\n"); continue;
		}
		memset(dp,0,sizeof(dp));
		dp[1][0]=1;
		for(int i=2;i<=n+1;i++)
		{
			for(int j=0;j<(1<<m);j++)
			{
				for(int k=0;k<(1<<m);k++)
				{
					if(dp[i-1][k]==0) continue;
					int tmp=test(j,k);
					if(!tmp) continue;
					dp[i][j]+=dp[i-1][k];
				}
			}
		}
		printf("%I64d\n",dp[n+1][0]);
	}
    return 0;
}