POJ 2411 Mondriaan's Dream

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 1000000000typedef __int64 LL;#define N 13int 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);#endifwhile(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;}