Mondriaan's Dream【状态压缩DP】

题目链接 POJ-2411

  题中很好的控制了N、M<=11所以挺明显的状压，又如何状压？我令横放的长方块为横着{1 1}令竖着的小方块为竖着{0 1}，然后找限制它们状态的因素，其中横着的，需要{11}这样连着的状态，而这一行为{0}限制的是下一行的状态，所以考虑为{0}位得向上考虑，或者这一位为0即可直接忽略，而看到下一位，用数学归纳法思想类推下去，最后一行得是什么个情况呢？假如有一个数比特位为0都不行，说明还没排完，所以最后一位全为1。

 

我写了两遍，所以附上两份代码，均有注释的：

 

代码一：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
using namespace std;
typedef long long ll;
const int maxN=11;
int N, M;
ll dp[maxN+2][(1<<maxN)+5];     //第0～M-1列
bool init(int state)        //第一列处理
{
    int i=0;
    while(i<M)
    {
        if(state & (1<<i) )        //第i列情况
        {
            if(i==M-1) return false;    //如果该列是最后一列，那么它没有后缀了，是多余的1
            if(state & (1<<(i+1)) ) i+=2;     //它的下一个也是1，说明是横着放置的，可行
            else return false;      //剩下的状态就是该列的下一列为0，不可行
        }
        else i++;       //这一列是0，竖放情况不用急于处理
    }
    return true;
}
bool check(int now, int last)       //第i列与第i-1列的状态是否满足
{
    int i=0;
    while(i<M)
    {
        int u1=last & (1<<i), u2=last & (1<<(i+1)), v1=now & (1<<i), v2=now & (1<<(i+1));
        if(u1 && v1)
        {
            if(i==M-1 || !u2 || !v2) return false;
            else i+=2;
        }
        else if(u1 && !v1) i++;
        else if(!u1 && v1) i++;
        else return false;
    }
    return true;
}
void solve()
{
    int tot=(1<<M)-1;       //满情况
    memset(dp, 0, sizeof(dp));
    for(int i=0; i<=tot; i++)
    {
        if(init(i)) dp[1][i]=1;
    }
    for(int i=2; i<=N; i++)
    {
        for(int j=0; j<=tot; j++)       //第i列情况
        {
            for(int k=0; k<=tot; k++)       //第i-1列情况
            {
                if(check(j, k)) dp[i][j]+=dp[i-1][k];
            }
        }
    }
    printf("%lld\n", dp[N][tot]);
}
int main()
{
    while(scanf("%d%d", &N, &M)!=EOF)
    {
        if(!N && !M) break;
        if((N*M)%2) { printf("0\n"); continue; }
        if(N<M) swap(N, M);
        solve();
    }
    return 0;
}

 

第二份：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
using namespace std;
typedef long long ll;
const int maxN=11;
int N, M, tot;
ll dp[maxN+1][(1<<maxN)+2];
bool init(int state)        //state为状态，目前这一行的状态，此间用以初始化第一排使用
{
    int i=0;
    while(i<M)
    {
        if(state & (1<<i))      //第i位z情况
        {
            if(i==M-1) return false;        //已经是最后一位了，没有后继无法成立
            if(state & (1<<(i+1)) ) i+=2;       //符合条件横放
            else return false;      //下一位是比特位0位，不行
        }
        else i++;
    }
    return true;        //遍历至最后，说明可行
}
bool check(int now, int pre)
{
    int i=0;
    while(i<M)
    {
        int u1=pre & (1<<i), u2=pre & (1<<(i+1)), v1=now & (1<<i), v2=now & (1<<(i+1));
        if(u1 && v1)        //前两列均为1
        {
            if(i==M-1) return false;        //到底，所以不行
            if(u2 && v2) i+=2;      //两个均横放
            else return false;      //其他都不行
        }
        else if(u1) i++;     //前一列为1，这一列为0
        else if(v1) i++;
        else return false;
    }
    return true;
}
void solve()
{
    for(int i=0; i<=tot; i++)
    {
        if(init(i)) dp[1][i]=1;
    }
    for(int i=2; i<=N; i++)     //第i行
    {
        for(int j=0; j<=tot; j++)       //目前行的状态
        {
            for(int k=0; k<=tot; k++)       //上一行的状态
            {
                if(check(j, k)) dp[i][j]+=dp[i-1][k];
            }
        }
    }
    printf("%lld\n", dp[N][tot]);
}
int main()
{
    while(scanf("%d%d", &N, &M) && (M | N))
    {
        if(N*M%2) { printf("0\n"); continue; }
        if(N<M) swap(N, M);
        tot=(1<<M)-1;       //求得满时上限
        memset(dp, 0, sizeof(dp));
        solve();
    }
    return 0;
}