Poj 2411 Mondriaan's Dream 状态压缩（方案数）

题目链接：http://poj.org/problem?id=2411

题目大意：给你一块N*M的地M，让你用1*2的地砖去铺满，求总的方案数。

解题思路：可以用深搜得到每行的合法状态，但本人弱菜，以后学习补上，直接暴力合法状态，具体解释见代码。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m;
long long  dp[12][1<<12],sta[1<<12],vis[12];
void init()
{
    memset(dp,0,sizeof(dp));
    memset(sta,0,sizeof(sta));
    memset(vis,0,sizeof(vis));
}
void solve()
{
    int num=(1<<m),i,j,k,p,ct=0;
    for(i=0;i<num;i++)   //筛选一行的合法状态数，其中1代表横放，则两个1相邻才合法
    {
        memset(vis,0,sizeof(vis));
        for(j=1; j<m; j++)           //j-1,故j从 1开始
            if( (i&(1<<j)) &&  (i&(1<<(j-1))) && !vis[j-1])
                vis[j-1]=vis[j]=1;
        for(j=0;j<m;j++)
            if( (i&(1<<j)) && !vis[j])
            break;
        if(j==m) sta[ct++]=i;
    }

    for(i=0;i<ct;i++)            //初始化，求方案数初始化为1
        dp[0][ sta[i] ]=1;

    for(i=1;i<n;i++)         //行数
    {
        for(j=0;j<num;j++)    //枚举i-1行状态
        {
            for(k=0;k<ct;k++)   //枚举第i行合法状态
            {
                int t=0;
                for(p=m-1;p>=0;p--)  //i行与i-1行状态对比，得到第i行状态
                {
                    //上下同为1，表示上下均横放，同为0 ，表示上下均竖放，均合法，且该位置已填满为1；
                    if( (j&(1<<p))==(sta[k]&(1<<p)) ) t=(t<<1)+1;
                    //上面为1，下面为0，表示上面横放，下面竖放，说明该位置未填为0；
                    else if( ((j&(1<<p))) &&  !(sta[k]&(1<<p)) )  t=(t<<1);
                    //上为0，下为1，说明下面横放，上面未放满，不合法，跳出
                    else break;
                }
                if(p==-1)  
                    dp[i][t]+=dp[i-1][j];
            }
        }
    }
    printf("%I64d\n",dp[n-1][num-1]);//最终状态为第n-1行铺满，即dp[n-1][num-1]
}
int main()
{
    while(scanf("%d%d",&n,&m)&&(n+m))
    {
        if(n<m) swap(n,m);
        init();
        solve();
    }
    return 0;
}