【Acwing.291 蒙德里安的梦想】状态压缩dp(动态规划篇)

AC wing.291 蒙德里安的梦想

求把 N×M 的棋盘分割成若干个 1×2 的的长方形，有多少种方案。

例如当 N=2，M=4 时，共有 5 种方案。当 N=2，M=3 时，共有 3 种方案。

如下图所示：

输入格式
输入包含多组测试用例。

每组测试用例占一行，包含两个整数 N 和 M。

当输入用例 N=0，M=0 时，表示输入终止，且该用例无需处理。

输出格式
每个测试用例输出一个结果，每个结果占一行。

数据范围
1≤N,M≤11
输入样例：
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

核心思想：
状态压缩dp
把题意读明白了可以想到如果把横着的摆上了，那么竖着的就已经能确定了，所以就只需要看有多少横着的就行，看下面这张图，可以枚举上一列和这一列的所有情况，就是如果某一行有的话，那么二进制中相对应的位置就应该填1，否则就要填0，那么怎么从上一个状态转移到下一个状态呢，就是把上一列的状态用k变量枚举出来，这一列的状态用变量j枚举出来，那么有效的方案数必须满足两个条件，一个是j & k == 0 , 另一个是j | k 不存在连续的奇数次0，可以先打个表，把每个状态有无连续奇数次0给打出来，然后进行dp的核心部分，枚举每一列，最后答案是f[m][0],因为只有第m-1列横着的小方格全是0的时候才是合法方案，至于为什么从第一列开始，因为第0列他不可能放开啊，他就1列，不可能放得开2列的东西，注意咱们是从第0列到第m-1列是主角

具体看代码：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N = 12;
long long f[N][1 << N], n, m;
bool st[1 << N];

int main()
{
    while (cin >> n >> m, n || m)
    {
        memset(f, 0, sizeof f);
        for (int i = 0; i < 1 << n; i++)
        {
            st[i] = true;
            int cnt = 0;
            for (int j = 0; j < n; j++)
            {
                if (i >> j & 1)
                {
                    if (cnt & 1)
                        st[i] = false;
                    cnt = 0;
                }
                else
                    cnt++;
            }
            if(cnt & 1)
                st[i] = false;
        }
        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 << m; k++)
                {
                    if ((j & k) == 0 && st[j | k])
                        f[i][j] += f[i - 1][k];
                }
            }
        }
        cout << f[m][0] << endl;
    }
    return 0;
}