Mondriaan's Dream

最新推荐文章于 2021-11-18 20:39:39 发布

原创最新推荐文章于 2021-11-18 20:39:39 发布 · 459 阅读

0 ·

CC 4.0 BY-SA版权

DP_状态压缩dp 专栏收录该内容

11 篇文章

订阅专栏

本文介绍了一种使用状压动态规划方法解决矩形填充问题的经典案例。该问题要求计算不同方式填充给定尺寸的大矩形，利用2x1的小矩形，考虑所有可能的排列组合。通过按列进行递归搜索并采用状态压缩技巧，文章提供了一个高效的解决方案。

这道题堪称经典的状压dp的题目了。

题目：

quares and rectangles fascinated the famous Dutch painter Piet Mondriaan. One night, after producing the drawings in his 'toilet series' (where he had to use his toilet paper to draw on, for all of his paper was filled with squares and rectangles), he dreamt of filling a large rectangle with small rectangles of width 2 and height 1 in varying ways.

Expert as he was in this material, he saw at a glance that he'll need a computer to calculate the number of ways to fill the large rectangle whose dimensions were integer values, as well. Help him, so that his dream won't turn into a nightmare!

Input

The input contains several test cases. Each test case is made up of two integer numbers: the height h and the width w of the large rectangle. Input is terminated by h=w=0. Otherwise, 1<=h,w<=11.

Output

For each test case, output the number of different ways the given rectangle can be filled with small rectangles of size 2 times 1. Assume the given large rectangle is oriented, i.e. count symmetrical tilings multiple times.

Sample Input

Sample Output

题目就像图中的意思一样。就是有多少种方案。

题目分析：每一个位置都与上一个位置有关，而且也有可能影响其他的状态。

其实仔细分析就是三种：

 //如果这个位置被上一列占用，直接跳过
    if(((1<<j)&sta)>0)
    dfs(i,j+1,sta,next);
    //如果这个位置是空的，尝试放一个1*2的
    if(((1<<j)&sta)==0)
    {
        dfs(i,j+1,sta,next|(1<<j));
    }
    //如果这个位置和下个位置都是空的，尝试放一个2*1的
    if(j+1<N&&((1<<j)&sta)==0&&((1<<j+1)&sta)==0)
    {
        dfs(i,j+2,sta,next);
    }

这道题我是按列来做的。实际上按行也可以。
dp[i][j]用来表示第i列的状态为j。用j二进制表示各个位置放与不放

注意数组的大小。刚开始开成了dp[12][1<<12]反正就是小了

状态的储存和转移

#include<iostream>
#include<string.h>
using namespace std;
int N,M;
int i,j;
//long long dp[12][34];
long long dp[13][1<<12];
void dfs(int i,int j,int sta,int next)//第i列。第j行。sta当前列的状态，next对下一列的影响
{
    if(j==N)//j到了最后一行
    {
        dp[i+1][next]+=dp[i][sta];
        return ;
    }
    //如果这个位置被上一列占用，直接跳过
    if(((1<<j)&sta)>0)
    dfs(i,j+1,sta,next);
    //如果这个位置是空的，尝试放一个1*2的
    if(((1<<j)&sta)==0)
    {
        dfs(i,j+1,sta,next|(1<<j));
    }
    //如果这个位置和下个位置都是空的，尝试放一个2*1的
    if(j+1<N&&((1<<j)&sta)==0&&((1<<j+1)&sta)==0)
    {
        dfs(i,j+2,sta,next);
    }
    return;

}
int main()
{
    while(cin>>N>>M)
    {
        memset(dp,0,sizeof(dp));
        if(N==0&&M==0)
            break;
        dp[1][0]=1;//第一列开始不铺
        for(int i=1;i<=M;i++)
        {
            for(j=0;j<(1<<N);j++)
            {
                if(dp[i][j])
                {
                    dfs(i,0,j,0);
                }
            }
        }
        cout<<dp[M+1][0]<<endl;
    }
}