【算法设计】背包问题

研究生课程系列文章参见索引《在信科的那些课》

题目

一个旅行者准备随身携带一个背包,可以放入背包的物品有n种,每种物品的重量和价值分别为wj, vj . 如果背包的最大重量限制是b, 怎样选择放入背包的物品以使得背包的价值最大?
目标函数:

约束条件:


算法设计


设Fk(y) 表示只允许装前k 种物品,背包总重不超过y 时背包的最大价值。Fk(y)有两种情况:不装第k件物品或至少装1件第k种物品。
如果不装第k件物品,那么只能用前k-1件物品装入背包,背包的限制重量仍为y,所以最大价值是Fk-1(y);
如果装1件第k件物品,那么装入的第k件物品价值为vk,重量为wk,剩下的物品仍要在前k件里选择(因为每件物品可以装多件,如果只能装1件就是在前k-1件里选择)。于是问题规约为背包限制重量y-wk的情况下前k件物品取得最大价值,即Fk(y-wk)+vk。 递推方程与边界条件:



上式初值比较多,F0(y)是不装物品的最大价值;F1(y)是只能装第一件物品时,最多装|_y/w1_|件;递推式Fk(y-wk),y-wk有可能得到负值,即不能再装物品,所以设置最小数以保证在优化问题中淘汰这种情况。

算法实现

标记函数:实现是需要一个ik(y)记录优化函数Fk(y)用到物品的最大标号。计算Fk(y)时,如果Fk-1(y)>Fk(y-wk)+vk,即没有加入第k件物品,ik(y)即为Fk-1(y)的物品最大标号;反正,加入第k件物品,ik(y)记为k。标记函数递归关系:


代码如下:
void Knapsack(int v[N],int w[N],int F[][B+1],int tagi[][B+1]){
	for(int k=0;k<=N;k++){
		F[k][0]=0;
		tagi[k][0]=0;
	}
	for(int y=0;y<=B;y++){
		F[0][y]=0;
		F[1][y]=(int)(y/w[0])*v[0];//只能装第一件物品时
		tagi[0][y]=0;
	}

	for(int k=1;k<=N;k++){
		for(int y=1;y<=B;y++){
			if(y-w[k-1]<0){
				F[k][y]=F[k-1][y];
				tagi[k][y]=tagi[k-1][y]; 
			}
			else{
				//允许装入k件物品,价值的两种情况:
				//不装第k件物品或至少装1件第k件物品
				F[k][y]=F[k-1][y]>F[k][y-w[k-1]]+v[k-1] ? F[k-1][y]:(F[k][y-w[k-1]]+v[k-1]);
				tagi[k][y]=F[k-1][y]>F[k][y-w[k-1]]+v[k-1]?tagi[k-1][y]:k;
			}
		}
	}
}

实例及解的追踪

试验下面的例子:
v1=1,v2=3,v3=5,v4=9,w1=2,w2=3,w3=4,w4=7,b=10
运行结果如下图:


我们需要在标记函数ik(y)中把实际解,及每个物品分别装入多少件追踪出来。
由最后i4(10)开始,i4(10)=4,表示此时第4件物品至少装入1件,占用重量w4=7,于是背包剩余重量为10-7=3;继续查询i4(3),由i4(3)=2,表示剩余物品最大标号为2,第2件物品至少装入1件。剩余重量为0,即不能再装入物品。用公式表示 追踪解的过程:


根据实例,可以理出追踪解的思路,代码如下:
void TrackSolution(int v[N],int w[N],int tagi[][B+1]){
	//x[i-1]标记第i件物品的件数
	int x[N];
	for(int i=0;i
运行结果:


其他题目

背包问题是很经典的动态规划问题,很多问题都是背包的变种,比如下面两个题目:
  • 设P是一台Internet上的Web服务器。T={1,2,...,n}是n个下载请求集合,ai为正整数,表示下载请求i所申请的带宽,已知服务器的最大带宽是正整数K。我们的目标是使带宽得到最大限度的利用,即确定T的一个子集S,使得,且达到最小。设计一个算法求服务器下载问题。
  • 设有n项任务,加工时间分别表示为正整数t1,t2,...,tn。现有2台同样的机器,从0时刻可以安排对这些任务的加工,知道T时刻所有任务完成,总加工时间为T。设计算法使得总加工时间T最小的调度方案。
第一个题目其实就是0-1背包问题,即看做价值和重量相等(都为ai)的物品装入背包,每件物品最多选1件,总重不能超过K,总价值最大的问题。设Fk(y)表示只允许前k个下载请求,最大带宽不超过y时利用最大限度的带宽数。 递推关系和边界条件如下:



注意0-1背包和背包问题的递推关系主要区别是:当选择第k件物品时,Fk(y)表示为Fk-1(y-wk)+vk,而非Fk(y-wk)+vk,即只能在前k-1件物品里继续选择。另外F1(y)的边界函数也不同。
至于第二个题目,其实就是使得一条加工线上的加工时间不超过T/2时加工时间尽可能大的问题,和第一个问题是一样的。


代码下载: http://download.csdn.net/detail/xiaowei_cqu/4775367
参考资料:屈婉玲 刘田等 《算法设计与分析》

(转载请注明作者和出处:http://blog.csdn.net/xiaowei_cqu未经允许请勿用于商业用途)



你可能感兴趣的:(【算法设计】背包问题)