hiho一下 第九周 状压dp

题目1 : 状态压缩·二

时间限制:10000ms

单点时限:1000ms

内存限制:256MB

描述

历经千辛万苦，小Hi和小Ho终于到达了举办美食节的城市！虽然人山人海，但小Hi和小Ho仍然抑制不住兴奋之情，他们放下行李便投入到了美食节的活动当中。美食节的各个摊位上各自有着非常多的有意思的小游戏，其中一个便是这样子的：

小Hi和小Ho领到了一个大小为N*M的长方形盘子，他们可以用这个盒子来装一些大小为2*1的蛋糕。但是根据要求，他们一定要将这个盘子装的满满的，一点缝隙也不能留下来，才能够将这些蛋糕带走。

这么简单的问题自然难不倒小Hi和小Ho，于是他们很快的就拿着蛋糕离开了~

但小Ho却不只满足于此，于是他提出了一个问题——他们有多少种方案来装满这个N*M的盘子呢？

值得注意的是，这个长方形盘子的上下左右是有区别的，如在N=4, M=3的时候，下面的两种方案被视为不同的两种方案哦！

提示：我们来玩拼图吧！不过不同的枚举方式会导致不同的结果哦！

输入

每个测试点（输入文件）有且仅有一组测试数据。

每组测试数据的第一行为两个正整数N、M，表示小Hi和小Ho拿到的盘子的大小。

对于100%的数据，满足2<=N<=1000, 3<=m<=5。<>

输出

考虑到总的方案数可能非常大，只需要输出方案数除以1000000007的余数。

 

样例输入

2 4

样例输出

5

比赛已经结束，去题库提交。

解题思路：每一行只跟上一行有关。

所以只要找到相邻两行所有相容的组合情况即可。

有两种方法第一种是枚举所有的情况，即压缩后的二进制数，然后判断是否矛盾。

第二种方法较为简单就是dfs枚举所有的可行情况。

对于相邻的两行有一下两种情况。

第二行第j位为0，则上一行该位置一定为1.

第二行第j位为1，则上一行该位置为0

第二行第j位为1，第j+1位也为1，上一行也是两个1.

三种情况可分别理解为向下方，向上放和横放。

设dp[i][s]为前i-1行放置满且第i行的状态为s的方案数。

所以可以设dp[0][(1<<m)-1]=1。

然后向下递推就可以。

要求的就是dp[n][1<<m-1]；

刚开始对相容的情况枚举的不太对，认为必须装满两行的情况才是相容。

其实不然。因为计算的时候我们只关心相邻两行的情况。

比如一下这种情况。

00001

11110

这种情况我本来认为是不相容的，因为第一行只有1个1.不可能存在这种情况。

实则不然

因为上上一行有可能是

11110

00001

11110

像这种，就可以使第二行为00001。

所以这种情况其实是相容的。

poj2411

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define LL long long
#define sca(x) scanf("%d",&x)

const LL mod = 1e9+7;
bool can[1<<11][1<<11];
LL dp[1<<11][1<<11];

void dfs(int col,int pre,int now,int m)
{
    if(col>m)return ;
    if(col==m)
    {
        can[pre][now]=1;
        return ;
    }
    dfs(col+1,pre<<1,now<<1|1,m);
    dfs(col+1,pre<<1|1,now<<1,m);
    dfs(col+2,(pre<<2)|3,(now<<2)|3,m);
}

int main()
{
    int n,m;
    while(cin>>n>>m){
    if(n==0&&m==0)break;
    memset(can,0,sizeof(can));
    memset(dp,0,sizeof(dp));
    dfs(0,0,0,m);
    dp[0][(1<<m)-1]=1;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<(1<<m);j++)
        {
            for(int k=0;k<(1<<m);k++)
            {
                if(can[j][k])
                {
                    dp[i][k]+=dp[i-1][j];
                }
            }
        }
    }
    printf("%lld\n",dp[n][(1<<m)-1]);
    }
}