0-1背包问题

版权声明:本文为博主原创文章,转载请注明出处,谢谢! https://blog.csdn.net/xp731574722/article/details/70766804

0-1 背包问题:给定 n 种物品和一个容量为 C 的背包,物品 i 的重量是 wi,其价值为 vi 。

问:应该如何选择装入背包的物品,使得装入背包中的物品的总价值最大?


分析一波,面对每个物品,我们只有选择拿取或者不拿两种选择,不能选择装入某物品的一部分,也不能装入同一物品多次。


解决办法:声明一个 大小为  m[n][c] 的二维数组,m[ i ][ j ] 表示 在面对第 i 件物品,且背包容量为  j 时所能获得的最大价值 ,那么我们可以很容易分析得出 m[i][j] 的计算方法,

(1). j < w[i] 的情况,这时候背包容量不足以放下第 i 件物品,只能选择不拿

m[ i ][ j ] = m[ i-1 ][ j ]

(2). j>=w[i] 的情况,这时背包容量可以放下第 i 件物品,我们就要考虑拿这件物品是否能获取更大的价值。

如果拿取,m[ i ][ j ]=m[ i-1 ][ j-w[ i ] ] + v[ i ]。 这里的m[ i-1 ][ j-w[ i ] ]指的就是考虑了i-1件物品,背包容量为j-w[i]时的最大价值,也是相当于为第i件物品腾出了w[i]的空间。

如果不拿,m[ i ][ j ] = m[ i-1 ][ j ] , 同(1)

究竟是拿还是不拿,自然是比较这两种情况那种价值最大。


由此可以得到状态转移方程:


 
 
  1. if(j>=w[i])
  2. m[i][j]=max(m[i -1][j],m[i -1][j-w[i]]+v[i]);
  3. else
  4. m[i][j]=m[i -1][j];


例:0-1背包问题。在使用动态规划算法求解0-1背包问题时,使用二维数组m[i][j]存储背包剩余容量为j,可选物品为i、i+1、……、n时0-1背包问题的最优值。绘制

价值数组v = {8, 10, 6, 3, 7, 2},

重量数组w = {4, 6, 2, 2, 5, 1},

背包容量C = 12时对应的m[i][j]数组。

0123456789101112
1000888888888
20008810101010181818
30668814141616181824
40669914141717191924
50669914141717192124
626891114161719192124
(第一行和第一列为序号,其数值为0)
如m[2][6],在面对第二件物品,背包容量为6时我们可以选择不拿,那么获得价值仅为第一件物品的价值8,如果拿,就要把第一件物品拿出来,放第二件物品,价值10,那我们当然是选择拿。m[2][6]=m[1][0]+10=0+10=10;依次类推,得到m[6][12]就是考虑所有物品,背包容量为C时的最大价值。


 
 
  1. #include <iostream>
  2. #include <cstring>
  3. using namespace std;
  4. const int N= 15;
  5. int main()
  6. {
  7.     int v[N]={ 0, 8, 10, 6, 3, 7, 2};
  8.     int w[N]={ 0, 4, 6, 2, 2, 5, 1};
  9.     int m[N][N];
  10.     int n= 6,c= 12;
  11.     memset(m, 0, sizeof(m));
  12.     for( int i= 1;i<=n;i++)
  13.     {
  14.         for( int j= 1;j<=c;j++)
  15.         {
  16.             if(j>=w[i])
  17.                 m[i][j]=max(m[i -1][j],m[i -1][j-w[i]]+v[i]);
  18.             else
  19.                 m[i][j]=m[i -1][j];
  20.         }
  21.     }
  22.     for( int i= 1;i<=n;i++)
  23.     {
  24.         for( int j= 1;j<=c;j++)
  25.         {
  26.             cout<<m[i][j]<< ' ';
  27.         }
  28.         cout<< endl;
  29.     }
  30.     return 0;
  31. }

到这一步,可以确定的是可能获得的最大价值,但是我们并不清楚具体选择哪几样物品能获得最大价值。

另起一个 x[ ] 数组,x[i]=0表示不拿,x[i]=1表示拿。

m[n][c]为最优值,如果m[n][c]=m[n-1][c] ,说明有没有第n件物品都一样,则x[n]=0 ; 否则 x[n]=1。当x[n]=0时,由x[n-1][c]继续构造最优解;当x[n]=1时,则由x[n-1][c-w[i]]继续构造最优解。以此类推,可构造出所有的最优解。(这段全抄算法书,实在不知道咋解释啊。。)


 
 
  1. void traceback()
  2. {
  3. for( int i=n;i> 1;i--)
  4. {
  5. if(m[i][c]==m[i -1][c])
  6. x[i]= 0;
  7. else
  8. {
  9. x[i]= 1;
  10. c-=w[i];
  11. }
  12. }
  13. x[ 1]=(m[ 1][c]> 0)? 1: 0;
  14. }

例:

某工厂预计明年有A、B、C、D四个新建项目,每个项目的投资额Wk及其投资后的收益Vk如下表所示,投资总额为30万元,如何选择项目才能使总收益最大?

Project

Wk

Vk

A

15

12

B

10

8

C

12

9

D

8

5

结合前面两段代码

 
 
  1. #include <iostream>
  2. #include <cstring>
  3. using namespace std;
  4. const int N= 150;
  5. int v[N]={ 0, 12, 8, 9, 5};
  6. int w[N]={ 0, 15, 10, 12, 8};
  7. int x[N];
  8. int m[N][N];
  9. int c= 30;
  10. int n= 4;
  11. void traceback()
  12. {
  13. for( int i=n;i> 1;i--)
  14. {
  15. if(m[i][c]==m[i -1][c])
  16. x[i]= 0;
  17. else
  18. {
  19. x[i]= 1;
  20. c-=w[i];
  21. }
  22. }
  23. x[ 1]=(m[ 1][c]> 0)? 1: 0;
  24. }
  25. int main()
  26. {
  27. memset(m, 0, sizeof(m));
  28. for( int i= 1;i<=n;i++)
  29. {
  30. for( int j= 1;j<=c;j++)
  31. {
  32. if(j>=w[i])
  33. m[i][j]=max(m[i -1][j],m[i -1][j-w[i]]+v[i]);
  34. else
  35. m[i][j]=m[i -1][j];
  36. }
  37. } /*
  38. for(int i=1;i<=6;i++)
  39. {
  40. for(int j=1;j<=c;j++)
  41. {
  42. cout<<m[i][j]<<' ';
  43. }
  44. cout<<endl;
  45. }
  46. */
  47. traceback();
  48. for( int i= 1;i<=n;i++)
  49. cout<<x[i];
  50. return 0;
  51. }

输出x[i]数组:0111,输出m[4][30]:22。

得出结论:选择BCD三个项目总收益最大,为22万元。


不过这种算法只能得到一种最优解,并不能得出所有的最优解。


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值