Mondriaan‘s Dream(状压dp)

题目描述

正方形和长方形着迷荷兰著名画家彼得·蒙德里安。一天晚上，他做了个梦，想在厕纸上画出来，他梦中用小长方形以不同的方式填充一个大长方形，而每个小长方形的大小是2x1。

他是在这个材料中的专家，他一眼就看出他需要一台计算机来计算填充长方形的方法，这些长方形的尺寸也是整数值。帮助他，使他的梦不会变成一场噩梦！

输入

输入包含多个测试案例。每个测试案例是由两个整数：高度H和W的大矩形的宽度，且1 <= h,w <=11。当输入H = W = 0时表示终止

输出

对于每个测试用例，输出不同的填充方法数量

样本输入

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 << 11; 对于每一个数字1表示竖着放置。0表示横放置， 但横着有要占连续的两格， 故在0 ~ （1<<11)中有许多不合法的状态，经过打表可选出合法状态。第二点就是如果上一行某一列为1， 那么这一一行这一列只能为0而如果上一行某一列为0，那么这一行可能为0，也可能为1。设x 为第i - 1行， y为第i 行，那么必须满足 x & y == 0. 而且 x | y 中的0都为连续的偶数个。细节见代码
枚举法

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 1e5 + 10;
#define ll long long
int n, m;
ll f[12][1<<11];
bool in_s[1<<11];

int main()
{
    while(cin >> n >> m && n)
    {
        memset(in_s, 1, sizeof in_s);
        for(int i = 0; i < (1<<m); i++)
        {
            int num = 0;
            for(int j = 0; j < m; j++)
            {
                if(i>>j & 1)
                {
                    if(num % 2 == 1)
                    {
                        in_s[i] = 0;
                        break;
                    }
                    num = 0;
                }
                else num ++;
            }
            if(num % 2 == 1) in_s[i] = 0;
        }
        f[0][0] = 1;
        for(int i = 1; i <= n; i++)
        {
            for(int j = 0; j < (1 << m); j++)
            {
                f[i][j] = 0;
                for(int k = 0; k < 1 << m; k ++)
                {
                    if((j & k) == 0 && in_s[j|k])
                        f[i][j] += f[i - 1][k];
                }
            }
        }
        cout << f[n][0] << endl;
    }
}

生成法

#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>

using namespace std;
vector<int> next_states;
int n, m;
long long f[12][1<<11];

void dfs(int state, int x, int ns)
{
    if(x == m)
    {
        next_states.push_back(ns);
        return;
    }
    if(state >> x & 1)
    {
        dfs(state, x + 1, ns);
        return;
    }
    dfs(state, x + 1, ns | 1 << x);
    if(x < m - 1 && (state>>(x + 1) & 1) == 0)
        dfs(state, x + 2, ns);
}

int main()
{
    while(cin >> n >> m && n)
    {
        memset(f, 0, sizeof f);
        f[0][0] = 1;
        for(int i = 0; i < n; i++)
        {
            for(int j = 0; j < 1 << m; j++)
            {
                if(f[i][j] > 0)
                {
                    next_states.clear();
                    dfs(j, 0, 0);
                    for(int k = 0; k < next_states.size(); k++)
                        f[i + 1][next_states[k]] += f[i][j];
                }
            }
        }
        cout << f[n][0] << endl;
    }
}