Domino 状压DP

Domino (domino.pas)

某天，wangjy和DaoThree在玩多米诺游戏。。。

DaoThree提出一个令wangjy抓头的问题：

有一个m行n列的矩阵，用1*2的骨牌(可横放或竖放)完全覆盖，骨牌不能重叠，有多少种不同的覆盖的方法？

请你帮助wangjy解决:)

wangjy不想看到长串的高精度数，你只需要求出覆盖方法总数mod p 的值即可

Input(domino.in)

三个整数数n,m,p

m<=5

p<=10000

Output(domino.out)

一个整数：总数模p的结果

我们发现m很小，于是可以进行状态压缩，用一个二进制串表示一个状态，0表示没有被覆盖，1表示被覆盖。我们可以先预处理出能够横放的所有状态，那么我们对于每个状态，如果我们知道了它哪些是横放的，那么我们也就可以推出它上一行的状态，因为当前行必须要把上一行补满。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 10005
int f[maxn][64];
int n,m,p,ans,list[65],tot;

bool check(int x)
{
	int ok=1;
	for (int i=0;i<m;i++) 
	{
		if ((!(x&(1<<i)))&&!ok) return 0;
		if ((x&(1<<i))) ok=ok^1;		
	}
	if (!ok) return 0;
	return 1;
}

int main()
{
	freopen("domino.in","r",stdin);
	freopen("domino.out","w",stdout);
	scanf("%d%d%d",&n,&m,&p);
	int r=1<<m;
	for (int i=0;i<r;i++) 
		if (check(i)) list[++tot]=i;
	for (int i=1;i<=tot;i++) f[1][list[i]]=1;
	for (int i=2;i<=n;i++) 
		for (int j=0;j<r;j++) 
			for (int k=1;k<=tot;k++) 
			{
				if ((j|list[k])==j) f[i][j]=(f[i][j]+f[i-1][list[k]^j^(r-1)])%p;	
			}
	printf("%d\n",f[n][r-1]);
	fclose(stdin);
	fclose(stdout);
	return 0;	
}