背包问题

问题描述：有不同价值、不同重量的物品n件，求从这n件物品中选取一部分物品的选择方案，使选中物品的总重量不超过指定的限制重量，但选中物品的价值之和最大。

设n个物品的重量和价值分别存储于数组w[ ]和v[ ]中，限制重量为tw。考虑一个n元组（x₀，x₁，…，x_n-1），其中x_i=0 表示第i个物品没有选取，而x_i=1则表示第i个物品被选取。显然这个n元组等价于一个选择方案。用枚举法解决背包问题，需要枚举所有的选取方案，而根据上述方法，我们只要枚举所有的n元组，就可以得到问题的解。

显然，每个分量取值为0或1的n元组的个数共为2ⁿ个。而每个n元组其实对应了一个长度为n的二进制数，且这些二进制数的取值范围为0～2ⁿ-1。因此，如果把0～2ⁿ-1分别转化为相应的二进制数，则可以得到我们所需要的2ⁿ个n元组。

【算法】

maxv=0;

for (i=0;i<2ⁿ;i++)

{     B[0..n-1]=0;

       把i转化为二进制数，存储于数组B中;

       temp_w=0;

       temp_v=0;

       for (j=0;j

       {     if (B[j]==1)

              {     temp_w=temp_w+w[j];

                     temp_v=temp_v+v[j];

              }

              if ((temp_w<=tw)&&(temp_v>maxv))

              {     maxv=temp_v;

                     保存该B数组；

              }

       }

}

 

设n件物品的重量分别为w₀、w₁、…、w_n-1，物品的价值分别为v₀、v₁、…、v_n-1。采用递归寻找物品的选择方案。设前面已有了多种选择的方案，并保留了其中总价值最大的方案于数组option[ ]，该方案的总价值存于变量maxv。当前正在考察新方案，其物品选择情况保存于数组cop[ ]。假定当前方案已考虑了前i-1件物品，现在要考虑第i件物品；当前方案已包含的物品的重量之和为tw；至此，若其余物品都选择是可能的话，本方案能达到的总价值的期望值为tv。算法引入tv是当一旦当前方案的总价值的期望值也小于前面方案的总价值maxv时，继续考察当前方案变成无意义的工作，应终止当前方案，立即去考察下一个方案。因为当方案的总价值不比maxv大时，该方案不会被再考察，这同时保证函数后找到的方案一定会比前面的方案更好。

对于第i件物品的选择考虑有两种可能：

（1）              考虑物品i被选择，这种可能性仅当包含它不会超过方案总重量限制时才是可行的。选中后，继续递归去考虑其余物品的选择。

（2）              考虑物品i不被选择，这种可能性仅当不包含物品i也有可能会找到价值更大的方案的情况。

按以上思想写出递归算法如下：

try(物品i，当前选择已达到的重量和，本方案可能达到的总价值tv)

{    

       if(包含物品i是可以接受的)

       {     将物品i包含在当前方案中；

              if (i

                     try(i+1,tw+物品i的重量,tv);

              else

                    

以当前方案作为临时最佳方案保存;

                     恢复物品i不包含状态；

              }

             

              if (不包含物品i仅是可男考虑的)

                     if (i

                            try(i+1,tw,tv-物品i的价值)；

                     else

                           

以当前方案作为临时最佳方案保存;

       }

       为了理解上述算法，特举以下实例。设有4件物品，它们的重量和价值见表：

物品	0	1	2	3
重量	5	3	2	1
价值	4	4	3	1

 

并设限制重量为7。则按以上算法，下图表示找解过程。由图知，一旦找到一个解，算法就进一步找更好的佳。如能判定某个查找分支不会找到更好的解，算法不会在该分支继续查找，而是立即终止该分支，并去考察下一个分支。

Try(0,0,12)

Try(1,5,12)

Try(1,0,8)

Try(2,5,8)

Try(3,7,8)

Try(2,3,8)

Try(3,5,8)

不能得到更好的解

不能得到更好的解

超重

不能得到更好的解

得到解：(1,0,1,0) maxv=7

得到解：(0,1,1,1) maxv=8

不能得到更好的解

超重

按上述算法编写函数和程序如下：

【程序】

# include

# define   N     100

double     limitW,totV,maxV;

int   option[N],cop[N];

struct      {     double     weight;

                     double     value;

              }a[N];

int   n;

void find(int i,double tw,double tv)

{     int k;

      

       if (tw+a[i].weight<=limitW)

       {     cop[i]=1;

              if (i find(i+1,tw+a[i].weight,tv);

              else

              {     for (k=0;k

                            option[k]=cop[k];

                     maxv=tv;

              }

              cop[i]=0;

}

      

       if (tv-a[i].value>maxV)

              if (i find(i+1,tw,tv-a[i].value);

              else

              {     for (k=0;k

                            option[k]=cop[k];

                     maxv=tv-a[i].value;

              }

}

 

void main()

{     int k;

       double w,v;

       printf(“输入物品种数\n”);

       scanf((“%d”,&n);

       printf(“输入各物品的重量和价值\n”);

       for (totv=0.0,k=0;k

       {     scanf(“% 1f%1f”,&w,&v);

              a[k].weight=w;

              a[k].value=v;

              totV+=V;

       }

       printf(“输入限制重量\n”);

       scanf(“%1f”,&limitV);

       maxv=0.0;

       for (k=0;k  cop[k]=0;

       find(0,0.0,totV);

       for (k=0;k

              if (option[k])  printf(“M”,k+1);

       printf(“\n总价值为%.2f\n”,maxv);

}

       作为对比，下面以同样的解题思想，考虑非递归的程序解。为了提高找解速度，程序不是简单地逐一生成所有候选解，而是从每个物品对候选解的影响来形成值得进一步考虑的候选解，一个候选解是通过依次考察每个物品形成的。对物品i的考察有这样几种情况：当该物品被包含在候选解中依旧满足解的总重量的限制，该物品被包含在候选解中是应该继续考虑的；反之，该物品不应该包括在当前正在形成的候选解中。同样地，仅当物品不被包括在候选解中，还是有可能找到比目前临时最佳解更好的候选解时，才去考虑该物品不被包括在候选解中；反之，该物品不包括在当前候选解中的方案也不应继续考虑。对于任一值得继续考虑的方案，程序就去进一步考虑下一个物品。

【程序】

# include

# define   N     100

double     limitW;

int   cop[N];

struct      ele   {     double     weight;

                            double     value;

                     } a[N];

int   k,n;

struct      {     int          flg;

                     double     tw;

                     double     tv;

              }twv[N];

void next(int i,double tw,double tv)

{     twv[i].flg=1;

       twv[i].tw=tw;

       twv[i].tv=tv;

}

double find(struct ele *a,int n)

{     int i,k,f;

       double maxv,tw,tv,totv;

       maxv=0;

       for (totv=0.0,k=0;k

              totv+=a[k].value;

       next(0,0.0,totv);

       i=0;

       While (i>=0)

       {     f=twv[i].flg;

              tw=twv[i].tw;

              tv=twv[i].tv;

              switch(f)

              {     case 1:     twv[i].flg++;

                                   if (tw+a[i].weight<=limitW)

                                          if (i

                                          {     next(i+1,tw+a[i].weight,tv);

                                                 i++;

                                          }

                                          else

                                          {     maxv=tv;

                                                 for (k=0;k

                                                        cop[k]=twv[k].flg!=0;

                                          }

                                   break;

                     case 0:     i--;

                                   break;

                     default:   twv[i].flg=0;

                                   if (tv-a[i].value>maxv)

                                          if (i

                                          {     next(i+1,tw,tv-a[i].value);

                                                 i++;

                                          }

                                          else

                                          {     maxv=tv-a[i].value;

                                                 for (k=0;k

                                                        cop[k]=twv[k].flg!=0;

                                          }

                                   break;

              }

       }

       return maxv;

}

 

void main()

{     double maxv;

       printf(“输入物品种数\n”);

       scanf((“%d”,&n);

       printf(“输入限制重量\n”);

       scanf(“%1f”,&limitW);

printf(“输入各物品的重量和价值\n”);

       for (k=0;k

              scanf(“%1f%1f”,&a[k].weight,&a[k].value);

       maxv=find(a,n);

       printf(“\n选中的物品为\n”);

for (k=0;k

              if (option[k])  printf(“M”,k+1);

       printf(“\n总价值为%.2f\n”,maxv);

}