POJ2411-Mondriaan's Dream（轮廓线DP）

题目链接

题意

用2*1的骨牌填充满m*n的矩阵，问有多少种不同的填法。

思路

有、没有用1、0表示，状压DP的思想，但是区别是用的是2*1的骨牌，填一行后可能会有不填的情况，所以轮廓线DP（或插头DP）就是解决这种每次一行不是整个的问题。记录的状态，不是一整行的状态，而可能是一个斜着的状态，状态的第一个位置并不一定对应着每行的第一个位置。

轮廓线DP详解：https://blog.csdn.net/u013480600/article/details/19499899

图示：https://blog.csdn.net/lvmaooi/article/details/79702273

每次对（i,j）位置只有三种情况，不放，横放，竖放。其中横放是（i-1,j）和（i,j），竖放是（i,j-1）和（i,j），即不能占右边和下边的位置。

对于状态K5K4K3K2K1和O，O上面的K5如果是0，只能够竖放。

每次的状态转移，是对每一个小格子（i,j），不同的状态K5K4K3K2K1相应的有多少种情况。所以轮廓线DP的状态转移不是对每一行，而是把行拆开成轮廓，对于每一个格子，相应的轮廓状态的情况。

dp[2][1<<N]，1<<N是状态个数，2是滚动数组，因为只需要记录前一个格子对应的状态，然后推出下一个格子的状态，其他的都是中间需要的，相当于滚动数组了。

2*1骨牌填矩阵的问题，也可以用状态转移矩阵解决+快速幂解决，特别是对于m,n有一个较小，有一个很大（1e9）。

轮廓线DP时间复杂度O(mn  *  2^m) ，如果mn有一个特别大就不适用了。

---

#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
using namespace std;
const int N = 11;
ll dp[2][1<<N];
int m,n,cur;

void update(int a,int b)
{
    if(b&(1<<m)) dp[cur][b^(1<<m)] += dp[cur^1][a];//如果第一个是1，才更新
}

int main()
{
    while(~scanf("%d%d",&m,&n))
    {
        if(m==0&&n==0) break;
        if(n<m) swap(n,m);
        cur = 0;
        memset(dp[cur],0,sizeof(dp[cur]));
        dp[cur][(1<<m)-1] = 1;
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++){
                cur ^= 1;
                memset(dp[cur],0,sizeof(dp[cur]));
                for(int k=0;k<(1<<m);k++){
                    update(k,k<<1);//不放
                    if(i && ((k&(1<<(m-1)))==0) ) update(k,(k<<1)^(1<<m)^1);//竖放
                    if(j && ((k&1)==0) ) update(k,(k<<1)^3);//横放
                }
            }
        }
        printf("%lld\n",dp[cur][(1<<m)-1]);
    }
    return 0;
}