[POJ 2411] Mondriaan's Dream 多米诺问题

http://poj.org/problem?id=2411

题意：用1*2的骨牌填满一个n*m的矩形，问有多少种填法。

思路：状态压缩，然后dp，横着的骨牌定义为0 0，竖着的定义为0 1，如下图

将每一行都压成一个数，图中第一行就是 6 第二行就是 7，每一行的状态数就是 2^m 个, 例如m = 3, 那么没一行就可以是：
0 0 0
0 0 1
0 1 0
0 1 1
1 0 0
1 0 1
1 1 0
1 1 1
转移方程：如果状态 k 和状态 l 兼容而且状态 l 合法则：dp[i][k] += dp[i-1][l] ；
所谓的兼容就是如果 上一次某个位置是 0 则下一层对应位置就必须是 1 (因为上一层放了一个竖着的骨牌)，如果上一层某个位置是 1 而且下一层对应位置也是 1 则下一层左边那个也必须是 1 ( 下一层这个肯定是横着的骨牌）。

#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

typedef long long loli;

int n, m, len;
loli dp[12][1<<11];

bool fun(int x, int y)
{
    if((x | y) != len) //不能存在上一层是 0 下一层也是 0。
        return false;
    int temp = m;
    y = y & x;   //用于判断上一层是 1 下一层也是 1 
    for(int i = 0; i < temp; i++)  //有 temp 位
    {
        if(y & 1) //上一层是 1 下一层也是 1
        {
            y >>= 1;  //判断左边
            temp -= 1;
            if(i == temp || !(y & 1))  //到了边界或者左边这位是 0 则不兼容
                return false;
        }
        y >>= 1;
    }
    return true;
}

int main()
{
    while(cin>>n>>m && (n || m)){
        if(m > n){swap(m, n);}
        memset(dp, 0, sizeof(dp));
        len = (1 << m) - 1;
        dp[0][len] = 1; //第一行的前面一行定义为全 1 
        for(int i = 1; i <= n; i++) // 有 n 层
        {
            for(int k = 0; k <= len; k++)//上一层的len种状态
            {
                for(int l = 0; l <= len && dp[i-1][k]; l++) //下一层的len中状态
                {
                    if(fun(k, l)){
                        dp[i][l] += dp[i-1][k];  //如果兼容就累加
                    }
                }
            }
        }
        cout<<dp[n][len]<<endl;  //最后一层肯定是全 1
    }
    return 0;
}