POJ 2411 - 状压DP详解

poj:http://poj.org/problem?id=2411

这个题属于状态压缩的一个入门题,在讲解该题之间我们先简单的回顾一下状压DP。

注意再看这个题之前,一定要对动态规划有一定的掌握

状压DP:

与传统的DP一样,状压DP也是定义状态,通过前一阶段的状态,推出当前阶段的状态。但是有时候状态特别多,但是每个状态的决策又非常少,我们如果用多维数组来表示状态会造成空间浪费,也可能会空间爆炸。有时候我们就可以考虑用状压DP来做。比如每个状态的决策只有两个,但是状态的维度很多。这道题就是属于这种。

为了让大家更加深刻的理解,这里简单的讲解一下利用状压DP来解决背包问题,背包问题的决策也非常的少,就是放与不放的问题,其决策也只有两个。但是背包问题使用状压DP来解决有点得不偿失,但是可以帮助大家更好的理解背包问题。

背包问题,状压DP解决

现在假设背包容量为10,物品个数5。那么我们就可以使用5的二进制位来表示这5个物品的放置情况。

如果5个物品都不放入背包,那么5个二进制位就为00000

如果5个物品都放入背包,那么5个二进制位就为11111

那么所有物品的放置情况就在00000-11111这个区间之中。我们就用这些二进制数表示完了所有物品的放置情况,如果我们将这些二进制数转换为10进制的话,那么就可以用这个10进制数来表示所有物品的放置情况了,这就是状态压缩的和核心。

上面简单的介绍了状态的表示,那么现在就来介绍状态的转移。

dp[00001]这个状态能由哪个状态转换过来?肯定就是dp[00000]这个状态了,就可以得出dp[00001] = dp[00000]+w1

注意:dp这个数组表示的是某种物品放置状态的价值。dp[00001]就表示放置了第1个物品的最大价值。

dp[01001]可以由dp[00001] + w4 或者 dp[01000] + w1转换过来,没什么区别这两个。这样就得到了状态转移方程,就能得到最终的答案了。可能有人要问了,背包容量限制呢?其实与价值的状态转移是一样的,容量的状态转移也可以用一个数组来进行表示。

例如:我们用dp1这个数组来表示容量的状态转移,dp1[00001]就表示放置了第1个物品所占容量,最后可以这两个数组,取出满足背包容量的,背包中物品价值最大的那一个就为结果。

讲的不明白?看看代码吧

代码实现:

#include
#define max(x,y) x > y ? x : y
// 记录不同状态的价值,注意根据物品的多少来决定容量
int dp[1000];
// 记录不同状态所占容量
int dp1[1000] ;
int main()
{
	int n,m;
	while(scanf("%d%d",&n,&m) == 2)
	{
		int W[10],C[10];
		for(int i = 0; i < n; i++){
            scanf("%d %d",&W[i],&C[i]) ;
		}
		for(int i = 0; i < (1 << n); i++)
		{
			for(int j = 0; j < n; j++)
			{
				if(!(i&(1 << j)))		//状态转移
				{
					int temp = i | (1<

我们也可以看出,状压DP解决的问题一般来说规模较小,状态一般只有两种状态

 

好了铺垫的差不多了,接下来进入我们的正题吧!

POJ 2411

简单的翻译一下题目吧

给你一个n*m的方格矩阵,要求用1*2的多米诺骨牌去填充,问有多少种填充方式。

这个题就是一个经典的状压DP题目,题中的状态非常的多,但是状态的决策又很少。进入正题吧

首先如果(n*m)%2!=0,那么填充方式为0种,这个很好分析。

给定状态转移方程:

f(i,j) = (f(i-1,k)之和) k 是与j状态对应的上一行的不同状态的集合。

f(i,j)表示第i行的状态为j时,前i行的总方案数。

现在给定一个矩阵,我们在矩阵上用0或者1来进行填充,当然这个填充不能是随便填充,我们必须要满足题目中的条件

首先我们定义如下这种填充表示方式:如果一个骨牌是横着放的,那么它所在的两个方格都填充1.如果它是竖着放的,那么它所在的两个格子中,上面的那个填0,下面的这个填1.如下图所示:

我们用1来表示矩阵当前位置的骨牌是竖着放置的,用0来表示骨牌横着放置的或者当前位置没有放置骨牌。

我们可以看出,每一行的格式我们可以将其看成一个二进制数,然后就可以使用一个10进制的数来表示其状态了。

我们把给定的矩阵两行两行的来看待,因为两行之间是相互影响的,我们只需要关注两个行之间不同的状态。

例如:图中第1行为110011 通过分析我们知道它的下一行只能是001111或者111111或001100或111100,根据我上面给出的条件这个很好分析。这里是通过上一行分析当前这一行,其实如果我们确定了当前这一行,他的上一行也会有多种状态与之对应

只要我们能得出当前行的所有状态对应的前一行的状态是什么,我们就能使用上面的状态转移方程了。最终就能解决这个题了

怎么得出这个状态关系?

我们知道当前行的状态确定了之后,他的前一行会有多个与之对用的状态。

由当前行来确定其前一行是可以用暴力的方式来解决,暴力枚举出合法的状态,能得到答案。但是其实有更简单的方法。

我们不确定当前行的状态,而是将其看待成对当前行的当前列的操作(操作:即为放置骨牌,对当前位置来说骨牌是横放还是竖放),寻找出与之对应的前一行的当前列的操作。

对于每个两两相邻的两行来说,我们可以使用三个变量,c表示当前所在列,now表示当前行的状态,pre表示前一行的状态。

然后从每一列的开始位置开始进行操作,直到向后模拟到最后一位,模拟到最后一位之后,就确定了一种状态。

我们的目的就是对当前行的当前列实施不同的操作,使其枚举出所有状态,并且找到前一行对应的状态

对于每一种操作我们都能唯一确定一个与之对应的前一行的操作,最终能得到对应的行,具体的操作在代码中进行说明

注意:这种操作只是关注当前行和其前一行的当前操作的位,而不去关注其他的位的,也不去关注其他行。千万不要想的太多

 

代码实现:

#include
#include
#define ll long long
#define N 12

int w,h,t;
ll dp[N][2100];
int path[14000][2];

void dfs(int c,int now,int pre){
    if(c > w){
        return;
    }
    if(c == w){
        // 得到了一种状态。
        path[t][0] = pre;
        path[t++][1] = now;
        return;
   
    // 当前位置横放,那么他的前一行也是横放,因为没有其他操作能一次填充前一行的两个位置
    dfs(c+2,(now<<2)|3,(pre<<2)|3);
    // 当前位置竖放,那么他的前一行肯定就是不放,因为前一行的当前位是由当前行的当前位竖放的骨牌进行填充的
    dfs(c+1,(now<<1)|1,pre<<1);
    // 当前位置不放,那么他的前一行肯定是竖放,因为前一行只有一个位的位置,只有竖放才满足
    dfs(c+1,now<<1,(pre<<1)|1);
}

int main(){
    while(scanf("%d %d",&h,&w) == 2){
        if(h == 0 && w == 0){
            return 0;
        }
        if(h < w){
            int temp = h;
            h = w;
            w = temp;
        }
        t = 0;
        dfs(0,0,0);
        memset(dp,0,sizeof(dp));
        dp[0][(1<

 

你可能感兴趣的:(算法)