[poj2411] 蒙德里安的梦

---

题目描述

有一个W行H列的棋盘，你可以在里放1*2的长方形，长方形之间互相不重叠，问放满整个棋盘有多少种方案数。

---

输入格式

输入文件有多组数据，每组数据只有一行为两个整数W和H(1<=W,H<=11)。

---

输出格式

每组数据一行为方案总数，若不能够放满整个棋盘输出0。

---

样例数据

样例输入

1 2
1 3
1 4
2 2
2 3
2 4
2 11
4 11

样例输出

1
0
1
2
3
5
144
51205

---

题目分析

骨牌问题。
想法1：搜索
搜索会有2^121种状态，无论Dfs还是Bfs都超时无疑
想法2：动规
状态是什么？观察数据范围较小，使用状态压缩动态规划，简称状压Dp
以行为阶段，每个格子若已经被铺了则为1，否则为0，那么行的状态就可以用01串表示出来，对于每一个长方形可以横铺或竖铺，现在的问题是如何从上一行的i状态转移到下一行的j状态，这两个状态是否兼容
两行的状态无非是这4种情况：
上一行 0 0 1 1
下一行 0 1 0 1
观察若上一行是0，这一行也是0，那么上一行的0将永远不会被填，所以这种状态不合法
若上一行是0，这一行是1，那么可以竖着放填上上一行的0，是合法的
若上一行是1，这一行是0，可以交给下一行填上0，合法
若上一行是1，这一行也是1，那么这一行只能横着放，于是要检查后两位是否也是1，只有后两位全是1才合法
若状态转移合法，方案数累加。
故枚举状态一一检查即可，初始状态：第一行只能横着放，故有1必须以偶数出现

---

源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
inline const int Get_Int() {
    int num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
int n,m,cnt=0;
long long f[15][1<<12];
bool Check(int a,int b) { //检查当前状态为a，前一行为b是否兼容 
    for(int i=1; i<=m; i++) {
        if((a&(1<<(i-1)))==0&&(b&(1<<(i-1)))==0)return false; //都为0 
        if((a&(1<<(i-1)))&&(b&(1<<(i-1)))) { //都为1，检查下一位 
            if((a&(1<<i))&&(b&(1<<i)))i++;
            else return false;
        }
    }
    return true;
}
bool Pre(int x) {
    int cnt=0;
    while(x) {
        if(x&1)cnt++;
        else if(cnt&1)return 0;
        x>>=1;
    }
    return !(cnt&1);
}
int main() {
    while(scanf("%d%d",&n,&m)!=EOF) {
        if((n*m)&1) {
            puts("0");
            continue;
        }
        if(n<m)swap(n,m);
        cnt=1<<m; //总状态数 
        memset(f,0,sizeof(f));
        for(int i=0; i<cnt; i++)
            if(Pre(i))f[1][i]=1;
        for(int i=2; i<=n; i++) //行为阶段 
            for(int j=0; j<cnt; j++) //枚举第i行状态j 
                for(int k=0; k<cnt; k++) //枚举i-1行状态k 
                    if(Check(j,k))f[i][j]+=f[i-1][k];
        printf("%lld\n",f[n][cnt-1]);
    }
    return 0;
}

---