动态规划：蒙德里安的梦想

题目链接

输入样例：
1 2
1 3
1 4
2 2
2 3
2 4
2 11
4 11
0 0

输出样例：
1
0
1
2
3
5
144
51205

解题核心： 将矩阵看成N×M个小正方形， 如果将横的1×2举证嵌入其中，那么竖得2×1矩阵也就确定，方案数也就确定。总的来说，横矩阵的摆放方案数就是总方案数。

解题方法：

状态压缩dp

状态表示：

f[i][j]表示前i-1列都摆放完毕，并且从i-1列伸到第i列的横块状态为j的方案数。

举例： 对于5×5的矩阵， f[2][25]的矩阵是
注意：25的二进制表示为11001（1代表该位置有横块， 0代表没有）
行和列都从0开始计数

状态转移：
对于f[i][j]，是所有合法的f[i-1][k] k = 0, 1, 2, 3 … (1 << n-1)的总和

那么如何判断是否合法呢， 先给出结论：

1. j&k == 0
2. j|k的二进制表示中只要偶数个连续的0

对于上述例子中5×5的矩阵， f[2][25]而言， f[1][1]是不合法的。
如图

其中k代表的黑色横块和j代表的红色横块重叠。
所有其实j&k == 0的含义很简单， 所代表的就是没有重叠的横块。

那么对于f[1][2]也是不合法的
如图

其中第i-1列剩余的空位只要一个，无法摆放竖块。

所以，j|k只有偶数个0代表的就是可以合法放下竖块。

最后重复一遍， 对于f[i][j]，是所有合法的f[i-1][k] k = 0, 1, 2, 3 … (1 << n-1)。

注意事项：最后答案的表示是f[m][0], 所以在题目中的N的范围应该是12， 而不是11。
方案数比较大， 用long long存储。

AC代码：

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int N = 12;
const int M = 1 << 12;
ll f[N][M];
int n, m;
bool st[M];//st[i]表示i是否有偶数个零

void init()
{
    memset(f, 0, sizeof f);
    //预处理st[]
    for(int i = 0; i < 1 <<n; i++)
    {
        st[i] = true;
        int cnt = 0;//记录连续0的个数
        for(int j = 0; j < n; j++)
        {
            int k = i >> j & 1;
            if(k == 0) cnt++;
            else
            {
                if(cnt&1) st[i] = false;
                cnt = 0;
            }
        }
        if(cnt&1) st[i] = false;
    }
}

int main()
{


    while(scanf("%d%d", &n, &m) != EOF)
    {
        if(n == 0 && m == 0) break;
        init();
        f[0][0] = 1;//初始化（什么都不放也是一种方案）
        for(int i = 1; i <= m; i++)
        {
            for(int j = 0; j < 1<<n; j++)
            {
                for(int k = 0; k < 1<<n; k++)
                {
                    //加上合法方案
                    if((j&k) == 0 && st[j|k]) f[i][j] += f[i-1][k];
                }
            }
        }
        printf("%lld\n", f[m][0]);

    }
}