背包問題

背包問題

背包問題(Knapsack problem)是一種組合最佳化的NP完全問題。問題可以描述為:給定一組物品,每種物品都有自己的重量和價格,在限定的總重量內,我們如何選擇,才能使得物品的總價格最高。問題的名稱來源於如何選擇最合適的物品放置於給定背包中。相似問題經常出現在商業、組合數學,計算複雜性理論、密碼學和套用數學等領域中。也可以將背包問題描述為決定性問題,即在總重量不超過W的前提下,總價值是否能達到V?它是在1978年由Merkel和Hellman提出的。

背包問題已經研究了一個多世紀,早期的作品可追溯到1897年數學家托比亞斯·丹齊格(Tobias Dantzig,1884-1956)的早期作品,並指的是包裝你最有價值或有用的物品而不會超載你的行李的常見問題。

基本介紹

  • 中文名:背包問題
  • 外文名:Knapsack Problem
  • 套用領域:運籌學,套用數學,密碼學等
  • 提出者:Merkel和Hellman
套用,定義,基礎背包,題目,基本思路,空間複雜,示例程式,遞歸實現,程式,測試數據,總結,完全背包,題目,基本思路,簡單有效,轉為問題,實現,總結,多重問題,題目,基本算法,轉為問題,算法,小結,三種背包,問題,背包混合,小結,二維費用,問題,算法,限制,小結,分組背包,問題,算法,小結,依賴問題,簡化問題,算法,一般問題,小結,泛化物品,定義,泛化物品,問題泛化,小結,問法變化,輸出方案,求方案總,最優方案,小結,

套用

1998年的石溪布魯克大學算法庫的研究表明,在75個算法問題中,背包問題是第18個最受歡迎,第4個最需要解決的問題(前三為後kd樹,後綴樹和bin包裝問題)。
背包問題出現在各種領域的現實世界的決策過程中,例如尋找最少浪費的方式來削減原材料,選擇投資和投資組合,選擇資產支持資產證券化,和生成密鑰為Merkle-Hellman和其他背包密碼系統。
背包算法的一個早期套用是在測試的構建和評分中,測試者可以選擇他們回答哪些問題。對於小例子來說,這是一個相當簡單的過程,為測試者提供這樣的選擇。例如,如果考試包含12個問題,每個問題的價值為10分,測試者只需回答10個問題即可獲得100分的最高分。然而,在點值的異質分布的測試 - 即不同的問題值得不同的點值 - 更難以提供選擇。 Feuerman和Weiss提出了一個系統,其中學生被給予一個異質測試,共有125個可能的點。學生被要求儘可能回答所有的問題。在總點數加起來為100的問題的可能子集中,背包算法將確定哪個子集給每個學生最高的可能得分。

定義

我們有n種物品,物品j的重量為wj,價格為pj。
我們假定所有物品的重量和價格都是非負的。背包所能承受的最大重量為W。
如果限定每種物品只能選擇0個或1個,則問題稱為0-1背包問題。
可以用公式表示為:
最大化
受限於
如果限定物品j最多只能選擇bj個,則問題稱為有界背包問題。
可以用公式表示為:
最大化
受限於
如果不限定每種物品的數量,則問題稱為無界背包問題。
各類複雜的背包問題總可以變換為簡單的0-1背包問題進行求解。

基礎背包

題目

有N件物品和一個容量為V的背包。第i件物品的重量是w[i],價值是v[i]。求解將哪些物品裝入背包可使這些物品的重量總和不超過背包容量,且價值總和最大。

基本思路

這是最基礎的背包問題,特點是:每種物品僅有一件,可以選擇放或不放。
用子問題定義狀態:即f[i][v]表示前i件物品恰放入一個容量為v的背包可以獲得的最大價值。則其狀態轉移方程便是:
f[i][v]=max{ f[i-1][v], f[i-1][v-w[i]]+v[i] }。
可以壓縮空間,f[v]=max{f[v],f[v-w[i]]+v[i]}
這個方程非常重要,基本上所有跟背包相關的問題的方程都是由它衍生出來的。所以有必要將它詳細解釋一下:“將前i件物品放入容量為v的背包中”這個子問題,若只考慮第i件物品的策略(放或不放),那么就可以轉化為一個只牽扯前i-1件物品的問題。如果不放第i件物品,那么問題就轉化為“前i-1件物品放入容量為v的背包中”,價值為f[i-1][v];如果放第i件物品,那么問題就轉化為“前i-1件物品放入剩下的容量為v-w[i]的背包中”,此時能獲得的最大價值就是f [i-1][v-w[i]]再加上通過放入第i件物品獲得的價值v[i]。
注意f[v]有意義若且唯若存在一個前i件物品的子集,其費用總和為f[v]。所以按照這個方程遞推完畢後,最終的答案並不一定是f[N] [V],而是f[N][0..V]的最大值。如果將狀態的定義中的“恰”字去掉,在轉移方程中就要再加入一項f[v-1],這樣就可以保證f[N] [V]就是最後的答案。至於為什麼這樣就可以,由你自己來體會了。

空間複雜

以上方法的時間和空間複雜度均為O(N*V),其中時間複雜度基本已經不能再最佳化了,但空間複雜度卻可以最佳化到O(V)。
先考慮上面講的基本思路如何實現,肯定是有一個主循環i=1..N,每次算出來二維數組f[i][0..V]的所有值。那么,如果只用一個數組f [0..V],能不能保證第i次循環結束後f[v]中表示的就是我們定義的狀態f[i][v]呢?
f[i][v]是由f[i-1][v]和f [i-1][v-w[i]]兩個子問題遞推而來,能否保證在推f[v]時(也即在第i次主循環中推f[v]時)能夠得到f[v]和f[v -w[i]]的值呢?事實上,這要求在每次主循環中我們以v=V..0的順序推f[v],這樣才能保證推f[v]時f[v-w[i]]保存的是狀態f[i-1][v-c[i]]的值。偽代碼如下:
for i=1..N
for v=V..0
f[v]=max{f[v],f[v-w[i]]+v[i]};
其中的f[v]=max{f[v],f[v-w[i]]}一句恰就相當於我們的轉移方程f[i][v]=max{f[i-1][v],f[i-1][v-w[i]]},因為的
f[v-w[i]]就相當於原來的f[i-1][v-w[i]]。如果將v的循環順序從上面的逆序改成順序的話,那么則成了f[i][v]由f[i][v-w[i]]推知,與本題意不符,但它卻是另一個重要的背包問題P02最簡捷的解決方案,故學習只用一維數組解01背包問題是十分必要的。

示例程式

(完全背包問題求解)
#include<iostream>using namespace std;int findM(int N,int K,int G[],int W[]){    int *M=new int[N+1],i,j,k;         for(i=0;i<N+2;i++)          M[i]=0;      for(i=0;i<K;i++)      {                  for(j=N;j>=G[i];j--)               {                               for(k=1;j-k*G[i]>=0;k++)                       {                          M[j]=M[j]>k*W[i]+M[j-k*G[i]]?M[j]:k*W[i]+M[j-k*G[i]];                           }                   }        }             return M[N];}int main(){        int N,K,i;        while(cin>>N>>K)        {        int *G=new int[K];                int *W=new int[K];                for(i=0;i<K;i++)               cin>>G[i]>>W[i];                cout<<findM(N,K,G,W)<<endl;                delete []G;               delete []W;            }         return 0;}

遞歸實現

//現在設A[i][v]表示在剩餘空間為v時選取當前物品i的最大值,B[i][v]表示不選取當前物品i的最大值,所以總的最大值必然是max(A[n][v],B[n][v]),詳細程式見如下:
#include<fstream>#include<iostream>using namespace std;#define MAXSIZE 1000int A[MAXSIZE+1][MAXSIZE+1],B[MAXSIZE+1][MAXSIZE+1];int c[MAXSIZE+1],w[MAXSIZE+1];int F(int n ,int v){if(n==0)return 0;if(!A[n][v]&&v>=c[n])A[n][v]=F(n-1,v-c[n])+w[n];if(!B[n][v])B[n][v]=F(n-1,v);return A[n][v]>B[n][v]?A[n][v]:B[n][v];}int main(int argc,char*argv[]){int n,v;memset(A,0,sizeof(A));memset(B,0,sizeof(B));ifstreamin("in.txt");ofstreamout("out.txt");cin>>n>>v;for(int i=1;i<=n;i++)cin>>c[i]>>w[i];cout<<F(n,v);return 0;}

程式

程式一:
vari,j,v,n:longint;f,c,w:array[0..100] of longint;function max(a,b:longint):longint;beginif a>b then exit(a) else exit(b);end;beginread(n,v);fillchar(f,sizeof(f),0);for i:=1 to n doread(c[i],w[i]);for i:=1 to n dofor j:=v downto c[i] dof[j]:=max(f[j],f[j-c[i]]+w[i]);writeln(f[v]);end.
程式二(順推法):
var m,n,x,i:integer;c,w:array[1..30] of integer;f:array[0..30,0..300] of integer;function max(x,y:integer):integer;beginif x>y then max:=x else max:=y;end;beginreadln(n,m);for i:=1 to n doreadln(c[i],w[i]);for i:=1 to n dofor x:=1 to m doif x>=c[i] then f[i,x]:=max(f[i-1,x-c[i]]+w[i],f[i-1,x])else f[i,x]:=f[i-1,x];writeln(f[n,m]);end.

測試數據

//in.txt:
5 100
77 92
22 22
29 87
50 46
99 90
//out.txt
133
//in.txt:
8 200
79 83
58 14
86 54
11 79
28 72
62 52
15 48
68 62
//out.txt
334
C動態規划算法的實現(完整代碼)
#include<stdio.h>#include<malloc.h>typedefstruct{intobject;intweight;intvalue;}KnapSack;KnapSack*knapsack;//背包數組,用malloc或new動態創建intnum;//物體的個數intcontainer;//背包的最大容量int**array=NULL;//用來存放子問題的結果//動態創建背包voidCreate_KnapSack(){charc;printf("inputthenumberofobjects\n");scanf("%d",&num);knapsack=newKnapSack[num+1];printf("inputweightandvalueof%dobjects,like1:410\n",num);for(inti=1;i<=num;i++){scanf("%d%c%d%c%d",&knapsack[i].object,&c,&knapsack[i].weight,&c,&knapsack[i].value);getchar();//為了獲取空格或其他輸入,聲明下scanf挺噁心}intk=knapsack[num].value;printf("%d",k);printf("inputthevolumeoftheknapsack:\n");scanf("%d",&container);}//確定最優子問題voidResolve_KnapSack(){intk=knapsack[num].value;printf("%d",k);//創建動態二維數組m[num][container]array=(int**)malloc((num+1)*sizeof(int*));for(inti=0;i<=num;i++)array[i]=(int*)malloc((container+1)*sizeof(int));//for(intj=0;j<=container;j++)array[num][j]=(j>=knapsack[num].weight)?knapsack[num].value:0;//子問題的最優結果for(intm=num-1;m>0;m--)for(intn=0;n<=container;n++)if(n>knapsack[m].weight&&array[m+1][n]<=array[m+1][n-knapsack[m].weight]+knapsack[m].value)array[m][n]=array[m+1][n-knapsack[m].weight]+knapsack[m].value;//else包括兩種情況,共同點是該物體沒有被使用elsearray[m][n]=array[m+1][n];}//往回找,確定某個物體i是否被使用bool*Trace_back(){intc=container;bool*used;used=(bool*)malloc(sizeof(bool)*(num+1));for(inti=1;i<num;i++)if(array[i][c]==array[i+1][c])used[i]=0;else{used[i]=1;c-=knapsack[i].weight;}used[num]=(c>=knapsack[num].weight)?1:0;returnused;}//用來輸出被使用的物體及其相關值voidPrint_KnapSack(bool*used){printf("theobjectsusedasfollows:\n");for(inti=1;i<=num;i++)if(used[i])printf("%d:%d%d\n",knapsack[i].object,knapsack[i].weight,knapsack[i].value);}voidmain(){bool*used;Create_KnapSack();Resolve_KnapSack();used=Trace_back();Print_KnapSack(used);}

總結

0/1背包問題是最基本的背包問題,它包含了背包問題中設計狀態、方程的最基本思想,另外,別的類型的背包問題往往也可以轉換成0/1背包問題求解。故一定要仔細體會上面基本思路的得出方法,狀態轉移方程的意義,以及最後怎樣最佳化的空間複雜度。

完全背包

題目

有N種物品和一個容量為V的背包,每種物品都有無限件可用。第i種物品的體積是c,價值是w。將哪些物品裝入背包可使這些物品的體積總和不超過背包容量,且價值總和最大。

基本思路

這個問題非常類似於01背包問題,所不同的是每種物品有無限件。也就是從每種物品的角度考慮,與它相關的策略已並非取或不取兩種,而是有取0件、取1件、取2件……等很多種。如果仍然按照解01背包時的思路,令f[i,v]表示前i種物品恰放入一個容量為v的背包的最大權值。仍然可以按照每種物品不同的策略寫出狀態轉移方程,像這樣:f[i,v]=max{f[i,v-vi]+wi,f[i-1,v]}。這跟01背包問題一樣有O(N*V)個狀態需要求解,但求解每個狀態的時間則不是常數了,求解狀態f[v]的時間是O(v/c),總的複雜度是超過O(VN)的。
將01背包問題的基本思路加以改進,得到了這樣一個清晰的方法。這說明01背包問題的方程的確是很重要,可以推及其它類型的背包問題。但我們還是試圖改進這個複雜度

簡單有效

完全背包問題有一個很簡單有效的最佳化,是這樣的:若兩件物品i、j滿足c<=c[j]且w>=w[j],則將物品j去掉,不用考慮。這個最佳化的正確性顯然:任何情況下都可將價值小體積高的j換成物美價廉的i,得到至少不會更差的方案。對於隨機生成的數據,這個方法往往會大大減少物品的件數,從而加快速度。然而這個並不能改善最壞情況的複雜度,因為有可能特別設計的數據可以一件物品也去不掉。
完全背包還有另一種最佳化,代碼如下:
var p,t:array[1..10000] of integer;m,n:integer;min:longint;procedure init;var i:integer;beginreadln(m,n);min:=maxlongint;for i:=1 to n dobeginreadln(p[i],t[i]);if t[i]<min then min:=t[i];end;end;procedure qsort(l,r:integer);var i,j,x,temp:longint;begini:=l;j:=r;x:=t[(l+r)div2];while i<jdobeginwhile(i<j)and(t[i]<x) do inc(i);while(i<j)and(x<t[j]) do dec(j);if i<=j thenbegintemp:=t[i];t[i]:=t[j];t[j]:=temp;temp:=p[i];p[i]:=p[j];p[j]:=temp;inc(i);dec(j);end;end;if i<r then qsort(i,r);if l<j then qsort(l,j);end;function max(a,b:longint):longint;beginif a>b then max:=a else max:=b;end;procedure work;var f:array[0..10000] of longint;i,j:longint;beginfillchar(f,sizeof(f),0);for i:=min to m dobeginf[i]:=f[i-1];for j:=1 to n doif i-t[j]>=0 then f[i]:=max(f[i],f[i-t[j]]+p[j]) else break;end;writeln(f[m]);end;begininit;qsort(1,n);work;end.

轉為問題

既然01背包問題是最基本的背包問題,那么我們可以考慮把完全背包問題轉化為01背包問題來解。最簡單的想法是,考慮到第i種物品最多選V/c 件,於是可以把第i種物品轉化為V/c件體積及價值均不變的物品,然後求解這個01背包問題。這樣完全沒有改進基本思路的時間複雜度,但這畢竟給了我們將完全背包問題轉化為01背包問題的思路:將一種物品拆成多件物品。
更高效的轉化方法是:把第i種物品拆成體積為c*2^k、價值為w*2^k的若干件物品,其中k滿足c*2^k<V。這是二進制的思想,因為不管最優策略選幾件第i種物品,總可以表示成若干個2^k件物品的和。這樣把每種物品拆成O(log(V/c))件物品,是一個很大的改進。但我們有更優的O(VN)的算法。* O(VN)的算法 這個算法使用一維數組,先看偽代碼:<pre class"example"> for i=1..N for v=0..V f[v]=max{f[v],f[v-c]+w};
你會發現,這個偽代碼與P01的偽代碼只有v的循環次序不同而已。為什麼這樣一改就可行呢?首先想想為什麼P01中要按照v=V..0的逆序來循環。這是因為要保證第i次循環中的狀態f[v]是由狀態f[v-c]遞推而來。換句話說,這正是為了保證每件物品只選一次,保證在考慮“選入第i件物品”這件策略時,依據的是一個絕無已經選入第i件物品的子結果f[v-c]。而完全背包的特點恰是每種物品可選無限件,所以在考慮“加選一件第i種物品”這種策略時,卻正需要一個可能已選入第i種物品的子結果f[v-c],所以就可以並且必須採用v= 0..V的順序循環。這就是這個簡單的程式為何成立的道理。
這個算法也可以以另外的思路得出。例如,基本思路中的狀態轉移方程可以等價地變形成這種形式:f[v]=max{f[v],f[v-c]+w},將這個方程用一維數組實現,便得到了上面的偽代碼。

實現

vari,j,v,n:longint;f,c,w:array[0..100] of longint;functionmax(a,b:longint):longint;beginifa>bthenexit(a)elseexit(b);end;beginread(n,v);fillchar(f,sizeof(f),0);fori:=1tondoread(c[i],w[i]);fori:=1tondoforj:=c[i]tovdof[j]:=max(f[j],f[j-c[i]]+w[i]);writeln(f[v]);end.

總結

完全背包問題也是一個相當基礎的背包問題,它有兩個狀態轉移方程,分別在“基本思路”以及“O(VN)的算法“的小節中給出。希望你能夠對這兩個狀態轉移方程都仔細地體會,不僅記住,也要弄明白它們是怎么得出來的,最好能夠自己想一種得到這些方程的方法。事實上,對每一道動態規劃題目都思考其方程的意義以及如何得來,是加深對動態規劃的理解、提高動態規劃功力的好方法。

多重問題

題目

有N種物品和一個容量為V的背包。第i種物品最多有n件可用,每件體積是c,價值是w。求解將哪些物品裝入背包可使這些物品的體積總和不超過背包容量,且價值總和最大。

基本算法

這題目和完全背包問題很類似。基本的方程只需將完全背包問題的方程略微一改即可,因為對於第i種物品有n+1種策略:取0件,取1件……取 n件。令f[v]表示前i種物品恰放入一個容量為v的背包的最大權值,則:f[v]=max{f[v-k*c]+ k*w|0<=k<=n}。複雜度是O(V*∑n)。

轉為問題

另一種好想好寫的基該方法是轉化為01背包求解:把第i種物品換成n件01背包中的物品,則得到了物品數為∑n的01背包問題,直接求解,複雜度仍然是O(V*∑n)。
但是我們期望將它轉化為01背包問題之後能夠像完全背包一樣降低複雜度。仍然考慮二進制的思想,我們考慮把第i種物品換成若干件物品,使得原問題中第i種物品可取的每種策略——取0..n件——均能等價於取若干件代換以後的物品。另外,取超過n件的策略必不能出現。
方法是:將第i種物品分成若干件物品,其中每件物品有一個係數,這件物品的費用和價值均是原來的費用和價值乘以這個係數。使這些係數分別為 1,2,4,...,2^(k-1),n-2^k+1,且k是滿足n-2^k+1>0的最大整數。例如,如果n為13,就將這種物品分成係數分別為1,2,4,6的四件物品。
分成的這幾件物品的係數和為n,表明不可能取多於n件的第i種物品。另外這種方法也能保證對於0..n間的每一個整數,均可以用若干個係數的和表示,這個證明可以分0..2^k-1和2^k..n兩段來分別討論得出,並不難,希望你自己思考嘗試一下。
這樣就將第i種物品分成了O(log n)種物品,將原問題轉化為了複雜度為O(V*∑log n)的01背包問題,是很大的改進。

算法

多重背包問題同樣有O(VN)的算法。這個算法基於基本算法的狀態轉移方程,但套用單調佇列的方法使每個狀態的值可以以均攤O⑴的時間求解。由於用單調佇列最佳化的DP已超出了NOIP的範圍,故本文不再展開講解。我最初了解到這個方法是在樓天成的“男人八題”幻燈片上。

小結

這裡我們看到了將一個算法的複雜度由O(V*∑n)改進到O(V*∑log n)的過程,還知道了存在套用超出NOIP範圍的知識的O(VN)算法。希望你特別注意“拆分物品”的思想和方法,自己證明一下它的正確性,並用儘量簡潔的程式來實現。

三種背包

問題

如果將P01、P02、P03混合起來。也就是說,有的物品只可以取一次(01背包),有的物品可以取無限次(完全背包),有的物品可以取的次數有一個上限(多重背包)。應該怎么求解呢?

背包混合

考慮到在P01和P02中最後給出的偽代碼只有一處不同,故如果只有兩類物品:一類物品只能取一次,另一類物品可以取無限次,那么只需在對每個物品套用轉移方程時,根據物品的類別選用順序或逆序的循環即可,複雜度是O(VN)。偽代碼如下:
for i=1..N
if 第i件物品是01背包
for v=V..0
f[v]=max{f[v],f[v-c]+w};
else if 第i件物品是完全背包
for v=0..V
f[v]=max{f[v],f[v-c]+w};
再加上多重背包
如果再加上有的物品最多可以取有限次,那么原則上也可以給出O(VN)的解法:遇到多重背包類型的物品用單調佇列解即可。但如果不考慮超過NOIP範圍的算法的話,用P03中將每個這類物品分成O(log n)個01背包的物品的方法也已經很優了。

小結

有人說,困難的題目都是由簡單的題目疊加而來的。這句話是否公理暫且存之不論,但它在本講中已經得到了充分的體現。本來01背包、完全背包、多重背包都不是什麼難題,但將它們簡單地組合起來以後就得到了這樣一道一定能嚇倒不少人的題目。但只要基礎紮實,領會三種基本背包問題的思想,就可以做到把困難的題目拆分成簡單的題目來解決。

二維費用

問題

二維費用的背包問題是指:對於每件物品,具有兩種不同的費用;選擇這件物品必須同時付出這兩種代價;對於每種代價都有一個可付出的最大值(背包容量)。問怎樣選擇物品可以得到最大的價值。設這兩種代價分別為代價1和代價2,第i件物品所需的兩種代價分別為a和b。兩種代價可付出的最大值(兩種背包容量)分別為V和U。物品的價值為w。

算法

費用加了一維,只需狀態也加一維即可。設f[v]表示前i件物品付出兩種代價分別為v和u時可獲得的最大價值。狀態轉移方程就是:fi[i][v][u]=max{f[i-1][v][u],f[i-1][v-a[i]][u-b[i]]+w[i]}。如前述方法,可以只使用二維的數組:當每件物品只可以取一次時變數v和u採用逆序的循環,當物品有如完全背包問題時採用順序的循環。當物品有如多重背包問題時拆分物品。

限制

有時,“二維費用”的條件是以這樣一種隱含的方式給出的:最多只能取M件物品。這事實上相當於每件物品多了一種“件數”的費用,每個物品的件數費用均為1,可以付出的最大件數費用為M。換句話說,設f[v][m]表示付出費用v、最多選m件時可得到的最大價值,則根據物品的類型(01、完全、多重)用不同的方法循環更新,最後在f[0..V][0..M]範圍內尋找答案。
另外,如果要求“恰取M件物品”,則在f[0..V][M]範圍內尋找答案。

小結

事實上,當發現由熟悉的動態規劃題目變形得來的題目時,在原來的狀態中加一維以滿足新的限制是一種比較通用的方法。希望你能從本講中初步體會到這種方法。

分組背包

問題

有N件物品和一個容量為V的背包。第i件物品的費用是c,價值是w。這些物品被劃分為若干組,每組中的物品互相衝突,最多選一件。求解將哪些物品裝入背包可使這些物品的費用總和不超過背包容量,且價值總和最大。

算法

這個問題變成了每組物品有若干種策略:是選擇本組的某一件,還是一件都不選。也就是說設f[k][v]表示前k組物品花費費用v能取得的最大權值,則有f[k][v]=max{f[k-1][v],f[k-1][v-c]+w|物品i屬於第k組}。
使用一維數組的偽代碼如下:
for 所有的組k
for v=V..0
for 所有的i屬於組k (我覺得循環順序應改成這樣,大家可以看一下以前的版本自己判斷)
f[v]=max{f[v],f[v-c]+w}
另外,顯然可以對每組中的物品套用P02中“一個簡單有效的最佳化”。

小結

分組的背包問題將彼此互斥的若干物品稱為一個組,這建立了一個很好的模型。不少背包問題的變形都可以轉化為分組的背包問題(例如P07),由分組的背包問題進一步可定義“泛化物品”的概念,十分有利於解題。

依賴問題

簡化問題

這種背包問題的物品間存在某種“依賴”的關係。也就是說,i依賴於j,表示若選物品i,則必須選物品j。為了簡化起見,我們先設沒有某個物品既依賴於別的物品,又被別的物品所依賴;另外,沒有某件物品同時依賴多件物品。

算法

這個問題由NOIP2006金明的預算方案一題擴展而來。遵從該題的提法,將不依賴於別的物品的物品稱為“主件”,依賴於某主件的物品稱為“附屬檔案”。由這個問題的簡化條件可知所有的物品由若干主件和依賴於每個主件的一個附屬檔案集合組成。
按照背包問題的一般思路,僅考慮一個主件和它的附屬檔案集合。可是,可用的策略非常多,包括:一個也不選,僅選擇主件,選擇主件後再選擇一個附屬檔案,選擇主件後再選擇兩個附屬檔案……無法用狀態轉移方程來表示如此多的策略。(事實上,設有n個附屬檔案,則策略有2^n+1個,為指數級。)
考慮到所有這些策略都是互斥的(也就是說,你只能選擇一種策略),所以一個主件和它的附屬檔案集合實際上對應於P06中的一個物品組,每個選擇了主件又選擇了若干個附屬檔案的策略對應於這個物品組中的一個物品,其費用和價值都是這個策略中的物品的值的和。但僅僅是這一步轉化並不能給出一個好的算法,因為物品組中的物品還是像原問題的策略一樣多。
再考慮P06中的一句話:可以對每組中的物品套用P02中“一個簡單有效的最佳化”。這提示我們,對於一個物品組中的物品,所有費用相同的物品只留一個價值最大的,不影響結果。所以,我們可以對主件i的“附屬檔案集合”先進行一次01背包,得到費用依次為0..V-c所有這些值時相應的最大價值f'[0..V-c]。那么這個主件及它的附屬檔案集合相當於V-c+1個物品的物品組,其中費用為c+k的物品的價值為f'[k]+w。也就是說原來指數級的策略中有很多策略都是冗餘的,通過一次01背包後,將主件i轉化為 V-c+1個物品的物品組,就可以直接套用P06的算法解決問題了。

一般問題

更一般的問題是:依賴關係以圖論中“森林”的形式給出(森林即多叉樹的集合),也就是說,主件的附屬檔案仍然可以具有自己的附屬檔案集合,限制只是每個物品最多只依賴於一個物品(只有一個主件)且不出現循環依賴。
解決這個問題仍然可以用將每個主件及其附屬檔案集合轉化為物品組的方式。唯一不同的是,由於附屬檔案可能還有附屬檔案,就不能將每個附屬檔案都看作一個一般的01 背包中的物品了。若這個附屬檔案也有附屬檔案集合,則它必定要被先轉化為物品組,然後用分組的背包問題解出主件及其附屬檔案集合所對應的附屬檔案組中各個費用的附屬檔案所對應的價值。
事實上,這是一種樹形DP,其特點是每個父節點都需要對它的各個兒子的屬性進行一次DP以求得自己的相關屬性。這已經觸及到了“泛化物品”的思想。看完P08後,你會發現這個“依賴關係樹”每一個子樹都等價於一件泛化物品,求某節點為根的子樹對應的泛化物品相當於求其所有兒子的對應的泛化物品之和。

小結

NOIP2006的那道背包問題我做得很失敗,寫了上百行的代碼,卻一分未得。後來我通過思考發現通過引入“物品組”和“依賴”的概念可以加深對這題的理解,還可以解決它的推廣問題。用物品組的思想考慮那題中極其特殊的依賴關係:物品不能既作主件又作附屬檔案,每個主件最多有兩個附屬檔案,可以發現一個主件和它的兩個附屬檔案等價於一個由四個物品組成的物品組,這便揭示了問題的某種本質。
我想說:失敗不是什麼丟人的事情,從失敗中全無收穫才是。

泛化物品

定義

考慮這樣一種物品,它並沒有固定的費用和價值,而是它的價值隨著你分配給它的費用而變化。這就是泛化物品的概念。
更嚴格的定義之。在背包容量為V的背包問題中,泛化物品是一個定義域為0..V中的整數的函式h,當分配給它的費用為v時,能得到的價值就是h(v)。
這個定義有一點點抽象,另一種理解是一個泛化物品就是一個數組h[0..V],給它費用v,可得到價值h[V]。
一個費用為c價值為w的物品,如果它是01背包中的物品,那么把它看成泛化物品,它就是除了h(c)=w其它函式值都為0的一個函式。如果它是完全背包中的物品,那么它可以看成這樣一個函式,僅當v被c整除時有h(v)=v/c*w,其它函式值均為0。如果它是多重背包中重複次數最多為n的物品,那么它對應的泛化物品的函式有h(v)=v/c*w僅當v被c整除且v/c<=n,其它情況函式值均為0。
一個物品組可以看作一個泛化物品h。對於一個0..V中的v,若物品組中不存在費用為v的的物品,則h(v)=0,否則h(v)為所有費用為v的物品的最大價值。P07中每個主件及其附屬檔案集合等價於一個物品組,自然也可看作一個泛化物品。

泛化物品

如果面對兩個泛化物品h和l,要用給定的費用從這兩個泛化物品中得到最大的價值,怎么求呢?事實上,對於一個給定的費用v,只需枚舉將這個費用如何分配給兩個泛化物品就可以了。同樣的,對於0..V的每一個整數v,可以求得費用v分配到h和l中的最大價值f(v)。也即f(v)=max{h(k) +l(v-k)|0<=k<=v}。可以看到,f也是一個由泛化物品h和l決定的定義域為0..V的函式,也就是說,f是一個由泛化物品h和 l決定的泛化物品。
由此可以定義泛化物品的和:h、l都是泛化物品,若泛化物品f滿足f(v)=max{h(k)+l(v-k)|0<=k<=v},則稱f是h與l的和,即f=h+l。這個運算的時間複雜度是O(V^2)。
泛化物品的定義表明:在一個背包問題中,若將兩個泛化物品代以它們的和,不影響問題的答案。事實上,對於其中的物品都是泛化物品的背包問題,求它的答案的過程也就是求所有這些泛化物品之和的過程。設此和為s,則答案就是s[0..V]中的最大值。

問題泛化

一個背包問題中,可能會給出很多條件,包括每種物品的費用、價值等屬性,物品之間的分組、依賴等關係等。但肯定能將問題對應於某個泛化物品。也就是說,給定了所有條件以後,就可以對每個非負整數v求得:若背包容量為v,將物品裝入背包可得到的最大價值是多少,這可以認為是定義在非負整數集上的一件泛化物品。這個泛化物品——或者說問題所對應的一個定義域為非負整數的函式——包含了關於問題本身的高度濃縮的信息。一般而言,求得這個泛化物品的一個子域(例如0..V)的值之後,就可以根據這個函式的取值得到背包問題的最終答案。
綜上所述,一般而言,求解背包問題,即求解這個問題所對應的一個函式,即該問題的泛化物品。而求解某個泛化物品的一種方法就是將它表示為若干泛化物品的和然後求之。

小結

本講可以說都是我自己的原創思想。具體來說,是我在學習函式式編程的 Scheme 語言時,用函式編程的眼光審視各類背包問題得出的理論。這一講真的很抽象,也許在“模型的抽象程度”這一方面已經超出了NOIP的要求,所以暫且看不懂也沒關係。相信隨著你的OI之路逐漸延伸,有一天你會理解的。
我想說:“思考”是一個OIer最重要的品質。簡單的問題,深入思考以後,也能發現更多。

問法變化

以上涉及的各種背包問題都是要求在背包容量(費用)的限制下求可以取到的最大價值,但背包問題還有很多種靈活的問法。
例如,求解最多可以放多少件物品或者最多可以裝滿多少背包的空間。這都可以根據具體問題利用前面的方程求出所有狀態的值(f數組)之後得到。
還有,如果要求的是“總價值最小”“總件數最小”,只需簡單的將上面的狀態轉移方程中的max改成min即可。
下面說一些變化更大的問法。

輸出方案

一般而言,背包問題是要求一個最優值,如果要求輸出這個最優值的方案,可以參照一般動態規劃問題輸出方案的方法:記錄下每個狀態的最優值是由狀態轉移方程的哪一項推出來的,換句話說,記錄下它是由哪一個策略推出來的。便可根據這條策略找到上一個狀態,從上一個狀態接著向前推即可。
還是以01背包為例,方程為f[v]=max{f[v],f[v-c]+w}。再用一個數組g [v],設g[v]=0表示推出f[v]的值時是採用了方程的前一項(也即f[v]=f[v]),g[v]表示採用了方程的後一項。注意這兩項分別表示了兩種策略:未選第i個物品及選了第i個物品。那么輸出方案的偽代碼可以這樣寫(設最終狀態為f[N][V]):
i=N
v=V
while(i>0)
if(g[v]==0)
print "未選第i項物品"
else if(g[v]==1)
print "選了第i項物品"
v=v-c
另外,採用方程的前一項或後一項也可以在輸出方案的過程中根據f[v]的值實時地求出來,也即不須紀錄g數組,將上述代碼中的g [v]==0改成f[v]==f[v],g[v]==1改成f[v]==f[v-c]+w也可。
輸出字典序最小的最優方案
這裡“字典序最小”的意思是1..N號物品的選擇方案排列出來以後字典序最小。以輸出01背包最小字典序的方案為例。
一般而言,求一個字典序最小的最優方案,只需要在轉移時注意策略。首先,子問題的定義要略改一些。我們注意到,如果存在一個選了物品1的最優方案,那么答案一定包含物品1,原問題轉化為一個背包容量為v-c[1],物品為2..N的子問題。反之,如果答案不包含物品1,則轉化成背包容量仍為V,物品為2..N的子問題。不管答案怎樣,子問題的物品都是以i..N而非前所述的1..i的形式來定義的,所以狀態的定義和轉移方程都需要改一下。但也許更簡易的方法是先把物品逆序排列一下,以下按物品已被逆序排列來敘述。
在這種情況下,可以按照前面經典的狀態轉移方程來求值,只是輸出方案的時候要注意:從N到1輸入時,如果f[v]==f及f[v]==f[f-c]+w同時成立,應該按照後者(即選擇了物品i)來輸出方案。

求方案總

對於一個給定了背包容量、物品費用、物品間相互關係(分組、依賴等)的背包問題,除了再給定每個物品的價值後求可得到的最大價值外,還可以得到裝滿背包或將背包裝至某一指定容量的方案總數。
對於這類改變問法的問題,一般只需將狀態轉移方程中的max改成sum即可。例如若每件物品均是01背包中的物品,轉移方程即為f[v]=sum{f[v],f[v-c]+w},初始條件f[0][0]=1。
事實上,這樣做可行的原因在於狀態轉移方程已經考察了所有可能的背包組成方案。

最優方案

這裡的最優方案是指物品總價值最大的方案。還是以01背包為例。
結合求最大總價值和方案總數兩個問題的思路,最優方案的總數可以這樣求:f[v]意義同前述,g[v]表示這個子問題的最優方案的總數,則在求f[v]的同時求g[v]的偽代碼如下:
for i=1..N
for v=0..V
f[v]=max{f[v],f[v-c]+w}
g[v]=0
if(f[v]==f[v])
inc(g[v],g[v]
if(f[v]==f[v-c]+w)
inc(g[v],g[v-c])
如果你是第一次看到這樣的問題,請仔細體會上面的偽代碼。

小結

顯然,這裡不可能窮盡背包類動態規劃問題所有的問法。甚至還存在一類將背包類動態規劃問題與其它領域(例如數論、圖論)結合起來的問題,在這篇論背包問題的專文中也不會論及。但只要深刻領會前述所有類別的背包問題的思路和狀態轉移方程,遇到其它的變形問法,只要題目難度還屬於NOIP,應該也不難想出算法。
觸類旁通、舉一反三,應該也是一個OIer應有的品質吧。

相關詞條

熱門詞條

聯絡我們