码迷,mamicode.com
首页 > 其他好文 > 详细

动态规划之背包问题

时间:2019-12-15 14:22:12      阅读:97      评论:0      收藏:0      [点我收藏+]

标签:可见   max   经典   思路   模板   它的   its   百度百科   程序设计   

背包问题是动态规划的一个分支,这里先简单介绍一下动态规划的思想。

动态规划程序设计是对解最优化问题的一种途径、一种方法,而不是一种特殊算法。不像搜索或数值计算那样,具有一个标准的数学表达式和明确清晰的解题方法。动态规划程序设计往往是针对一种最优化问题,由于各种问题的性质不同,确定最优解的条件也互不相同,因而动态规划的设计方法对不同的问题,有各具特色的解题方法,而不存在一种万能的动态规划算法,可以解决各类最优化问题。因此在学习时,除了要对基本概念和方法正确理解外,必需具体问题具体分析处理,以丰富的想象力去建立模型,用创造性的技巧去求解。我们也可以通过对若干有代表性的问题的动态规划算法进行分析、讨论,逐渐学会并掌握这一设计方法。(以上来自百度百科)

简单来说动态规划是一个思想,而不是一个固定的算法模板,我们需要通过这种思想确定状态转移方程(一个好的状态),然后再求解。

动态规划适用于很多情况,其中一些情况被统一划分归类,而背包问题则是一种最为常见的动态规划问题。

背包问题又分为许多种,如:01背包,完全背包等

我们这里介绍最为常见的三种,分别是01背包,完全背包,和多重背包。

1.01背包

这是一个经典的动态规划问题,另外在贪心算法里也有背包问题,至于二者的区别在此就不做介绍了。

题目一般都是有N件物品和一个容量为V的背包。第i件物品的体积是v[i],价值是c[i],将哪些物品装入背包可使价值总和最大?

分析起来很简单,每一种物品都有两种可能,即放入背包或者不放入背包。可以用dp[i][j]表示第i件物品放入容量为j的背包所得的最大价值,则状态转移方程可以推出如下:

dp[i][j]=max{dp[i-1][j-v[i]]+c[i],dp[i-1][j]};

先放下代码:

for (int i = 1;i <= N;i++) //枚举物品  
{  
    for (int j = 0;j <= V;j++) //枚举背包容量  
    {  
        f[i][j] = f[i - 1][j];
        if (j >= v[i])
        {  
            f[i][j] = Max(f[i - 1][j],f[i - 1][j - v[i]] + c[i]);
        }
    }
}

 

我们可以发现0-1背包的状态转移方程 dp[i][j] = max{dp[i-1][j-w[i]]+v[i],dp[i-1][j]}的特点,当前状态仅依赖前一状态的剩余体积与当前物品体积v[i]的关系。根据这个特点,我们可以将dp降到一维即dp[j] = max{dp[j],dp[j-w[i]]+v[i]}。从这个方程中我们可以发现,有两个dp[j],但是要区分开。等号左边的dp[j]是当前i的状态,右边中括号内的dp[j]是第i-1状态下的值。

所以为了保证状态的正确转移,我们需要先更新等号左边中的dp[j](当前状态的dp[j])。

#include <iostream>
using namespace std;

#define MAXSIZE 100
int w[MAXSIZE];
int v[MAXSIZE];
int maxv;
int n;
int dp[MAXSIZE];

int max(int a, int b)
{
    if (a > b)
        return a;
    else
        return b;
}

int main()
{
    cin >> n >> maxv;
    for (int i = 1; i <= n; i++)
    {
        cin >> w[i] >> v[i];
    }
    for (int i = 0; i <= maxv; i++)
        dp[i] = 0;

    for (int i = 1; i <= n; i++)
    {
        //只有当j >= w[i],dp[j]才能进行选取最大值,否则dp[j]将不作更新,等于dp[i-1][j]。
        for (int j = maxv; j >= w[i]; j--)
        {
            dp[j] = max(dp[j], dp[j - w[i]] + v[i]);
        }
        
    }

    cout << dp[maxv] << endl;
    return 0;
}

对于01背包问题,我们通常使用的是第二种方法,相比之第一种,第二种的空间复杂度无疑是减小了许多。

现在我们来看这个状态转移方程:dp[j] = max(dp[j], dp[j - w[i]] + v[i]);

我们先分析上一个状态也就是dp[j-w[i]],我们将上一个状态转移到现在的状态dp[j]有三种情况:

首先在现在的状态dp[j],在不转移的情况下这也是一种情况;

第二种是dp[j-w[i]],即当前物品不放入上一个状态;

第三种就是dp[j-w[i]]+v[i],即当前物品放入上一个状态;

状态转移的过程则是要求我们在这三者当中取最大值,我们知道v[i]必定是大于等于0的,也就是说dp[j-w[i]]+v[i]必定是大于等于dp[j-w[i]]的,我们可以把两者合并为一种情况,然后是dp[j],为什么dp[j]也算是一种情况,我们回到for循环可以发现第一个物品枚举完之后有一些状态dp[k*w[1]]更新了,它们的值都变成v[i],然后我们在接下来的过程当中dp[j]可能是已经有了值的,通过状态转移方程我们可以知道这个值是当物品总重量达到j时刻物品总价值的最大值。所以我们需要将dp[j]也看作是一种情况。

2.完全背包

对于01背包问题,它的特点是:每一件物品之多只能选择一件,即在背包中该物品数量只有0和1两种情况。

现在扩展一下,有一个容积为V的背包,同时有n种物品,每种物品均有无数多个,并且每种物品的都有自己的体积和价值。求使用该背包最多能够装的物品价值总和。

这就是完全背包问题。

如果按照0-1背包的思路求解该问题,可设当前物品的体积为w,价值为v,考虑到背包中最多存放V/w件该物品,那么该物品的可选数量就为V/w件。依次可以对所有的物品进行拆分,最后对拆分的所有物品做0-1背包即可得到答案。但是,这样拆分会使物品数量大大增加,其时间复杂度为:O(V*∑ni=1(V/wi))。

可见,当V较大时每个物品的体积较小时,其时间复杂度会显著增大。所以将完全背包问题转化为0-1背包问题去解决的方法不可靠。

在0-1背包的解决算法中,其中一段代码是该算法的核心算法,如下:

struct Good{
    int w;
    int v;
}goods[101];
int dp[101][1001];
int n,S;//n表示有n个物品,S表示背包的最大容积
for (i = 1; i <= n; i++)
{
    for (j = S; j >= goods[i].w; j--)
    {
         dp[j] = max(dp[j], dp[j - goods[i].w] + goods[i].v);
    }      
}

在这段代码中,之所将j初始化为S,逆序循环更新状态是为了保证在更新dp[j]时,dp[j-goods[i].w]的状态尚未因为本次更新而发生改变,即等价于由

dp[i-1][j-goods[i].w]转移得到dp[i][j]。保证了更新dp[j]时,dp[j-goods[i].w]是没有放入物品i时的数据dp[i-1][j-goods[i].w]。

在解决完全背包问题时,可以借鉴这个思路。在完全背包中,每个物品可以被无限次选择,那么状态dp[i][j]恰好可以由可能已经放入物品i的状态dp[i][j-goods[i].w]转移而来。可以将上面的代码改写如下:

for (i = 1; i <= n; i++)
{
    for (j = goods[i].w; j <= S; j++)
        dp[j] = max(dp[j], dp[j - goods[i].w] + goods[i].v);
}

这样不需要将物品拆分,但是本质上并没什么区别.

3.多层背包

多重背包问题是0-1背包问题和完全背包问题的综合体,可以描述如下:从n种物品向容积为V的背包装入,其中每种物品的体积为w,价值为v,数量为k,问装入的最大价值总和?

我们知道0-1背包问题是背包问题的基础,所以在解决多重背包问题的时候,要将多重背包向0-1背包上进行转换。在多重背包问题中,每种物品有k个,可以将每种物品看作k种,这样就可以使用0-1背包的算法。但是,这样会增加数据的规模。因为该算法的时间复杂度为O(V*∑ni=1ki),所以要降低每种物品的数量ki。

代码和完全背包的代码基本一样,唯一的不一样就是完全背包的每一种物品数量不限.

#include<bits/stdc++.h>
use namespace std;
const int maxn=10005;
int w[maxn],v[maxn]
int dp[maxn];
int main()
{
    int n,maxv;
   cin>>n>>maxv; int index=0; for(int i=0;i<n;i++) { int a,b,c; cin>>a>>b>>c; while(a--) { w[index]=b; v[index++]=c; } } dp[0]=1; for(int i=0;i<index;i++) { for(int j=maxv;j>=w[i];j--) { dp[j]=max(dp[j],dp[j-w[i]]+v[i]); } }    cout<<dp[maxv]<<endl; return 0; }

  

对于完全背包,不推荐使用上面这个方法,下面有一种复杂度更小的代码:

#include<bits/stdc++.h>
using namespace std;
const int maxn=1005;
struct node{
	int w;
	int v;
}p[maxn];
int dp[maxn];

int main(void)
{
	int n,maxv;
	cin>>n>>maxv;
	int index=0;
	for(int i=0;i<n;i++)
	{
		int a,b,c;
		cin>>a>>b>>c;
		while(a>0) {
			if(a&1) {
				p[index].v=c;
				p[index++].w=b;
			}
			a/=2;
			b*=2;
			c*=2;
		}
	}
	dp[0]=0;
	for(int i=0;i<index;i++)
	{
		for(int j=maxv;j>=p[i].w;j--)
		{
			dp[j]=max(dp[j],dp[j-p[i].w+p[i].v]);
		}
	}
	cout<<dp[maxv]<<endl;
    return 0;
}

  这里在划分k个物品的时候用了小技巧,用二进制将所有的情况都表示了出来减少了情况的种类,复杂度相比上一个代码简单了许多,.尤其是在物品数量多的时候.

动态规划之背包问题

标签:可见   max   经典   思路   模板   它的   its   百度百科   程序设计   

原文地址:https://www.cnblogs.com/bqyb/p/12043812.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!