貪心算法 - 0/1背包問題


1、問題描述

     給定n種物品和一背包。物品i的重量是wi,其價值為vi,背包的容量為C。問:應如何選擇裝入背包的物品,使得裝入背包中物品的總價值最大?

     形式化描述:給定c >0, wi >0, vi >0 , 1≤i≤n.要求找一n元向量(x1,x2,…,xn,), xi∈{0,1}, ∋ ∑ wi xi≤c,且∑ vi xi達最大.即一個特殊的整數規划問題。

       2、最優性原理

     設(y1,y2,…,yn)是 (3.4.1)的一個最優解.則(y2,…,yn)是下面相應子問題的一個最優解:

     證明:使用反證法。若不然,設(z2,z3,…,zn)是上述子問題的一個最優解,而(y2,y3,…,yn)不是它的最優解。顯然有
                                    ∑vizi > ∑viyi   (i=2,…,n)
     且                           w1y1+ ∑wizi<= c
     因此                       v1y1+ ∑vizi (i=2,…,n) > ∑ viyi, (i=1,…,n) 
     說明(y1,z2, z3,…,zn)是(3.4.1)0-1背包問題的一個更優解,導出(y1,y2,…,yn)不是背包問題的最優解,矛盾。

       3、遞推關系

    設所給0-1背包問題的子問題

     

     的最優值為m(i,j),即m(i,j)是背包容量為j,可選擇物品為i,i+1,…,n時0-1背包問題的最優值。由0-1背包問題的最優子結構性質,可以建立計算m(i,j)的遞歸式:

     注:(3.4.3)式此時背包容量為j,可選擇物品為i。此時在對xi作出決策之后,問題處於兩種狀態之一:
    (1)背包剩余容量是j,沒產生任何效益;
    (2)剩余容量j-wi,效益值增長了vi ;
     使用遞歸C++代碼如下:

  

復制代碼
#include<iostream>
using namespace std; const int N=3; const int W=50; int weights[N+1]={0,10,20,30}; int values[N+1]={0,60,100,120}; int V[N+1][W+1]={0}; int knapsack(int i,int j) { int value; if(V[i][j]<0) { if(j<weights[i]) { value=knapsack(i-1,j); } else { value=max(knapsack(i-1,j),values[i]+knapsack(i-1,j-weights[i])); } V[i][j]=value; } return V[i][j]; } int main() { int i,j; for(i=1;i<=N;i++) for(j=1;j<=W;j++) V[i][j]=-1; cout<<knapsack(3,50)<<endl; cout<<endl; }
復制代碼

不使用遞歸的C++代碼:簡單一點的修改http://www.cppblog.com/Geek/archive/2009/12/02/102393.html

復制代碼
//3d10-1 動態規划 背包問題 #include <iostream> using namespace std; const int N = 4; void Knapsack(int v[],int w[],int c,int n,int m[][10]); void Traceback(int m[][10],int w[],int c,int n,int x[]); int main() { int c=8; int v[]={0,2,1,4,3},w[]={0,1,4,2,3};//下標從1開始 int x[N+1]; int m[10][10]; cout<<"待裝物品重量分別為:"<<endl; for(int i=1; i<=N; i++) { cout<<w[i]<<" "; } cout<<endl; cout<<"待裝物品價值分別為:"<<endl; for(int i=1; i<=N; i++) { cout<<v[i]<<" "; } cout<<endl; Knapsack(v,w,c,N,m); cout<<"背包能裝的最大價值為:"<<m[1][c]<<endl; Traceback(m,w,c,N,x); cout<<"背包裝下的物品編號為:"<<endl; for(int i=1; i<=N; i++) { if(x[i]==1) { cout<<i<<" "; } } cout<<endl; return 0; } void Knapsack(int v[],int w[],int c,int n,int m[][10]) { int jMax = min(w[n]-1,c);//背包剩余容量上限 范圍[0~w[n]-1] for(int j=0; j<=jMax;j++) { m[n][j]=0; } for(int j=w[n]; j<=c; j++)//限制范圍[w[n]~c]  { m[n][j] = v[n]; } for(int i=n-1; i>1; i--) { jMax = min(w[i]-1,c); for(int j=0; j<=jMax; j++)//背包不同剩余容量j<=jMax<c  { m[i][j] = m[i+1][j];//沒產生任何效益  } for(int j=w[i]; j<=c; j++) //背包不同剩余容量j-wi >c  { m[i][j] = max(m[i+1][j],m[i+1][j-w[i]]+v[i]);//效益值增長vi  } } m[1][c] = m[2][c]; if(c>=w[1]) { m[1][c] = max(m[1][c],m[2][c-w[1]]+v[1]); } } //x[]數組存儲對應物品0-1向量,0不裝入背包,1表示裝入背包 void Traceback(int m[][10],int w[],int c,int n,int x[]) { for(int i=1; i<n; i++) { if(m[i][c] == m[i+1][c]) { x[i]=0; } else { x[i]=1; c-=w[i]; } } x[n]=(m[n][c])?1:0; }
復制代碼

運行結果:

 

算法執行過程對m[][]填表及Traceback回溯過程如圖所示:

      從m(i,j)的遞歸式容易看出,算法Knapsack需要O(nc)計算時間; Traceback需O(n)計算時間;算法總體需要O(nc)計算時間。當背包容量c很大時,算法需要的計算時間較多。例如,當c>2^n時,算法需要Ω(n2^n)計算時間。


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM