1. 程式人生 > >01揹包問題(動態規劃)

01揹包問題(動態規劃)

揹包問題

最近剛學了01揹包問題,但是聽老師講再加上以前自己看書看的,發現有很多地方很容易搞混,原理就是劃分找動態轉移方程,但是寫程式時會遇到困難,趁著今天有空,就特意整理一下01揹包問題。

動態規劃的基本思想: 動態規劃演算法通常用於求解具有某種最優性質的問題。即我們平常所說的最優子結構性質。動態規劃演算法與分治法類似,其基本思想也是將待求解問題分解成若干個子問題,先求解子問題,然後從這些子問題的解得到原問題的解。與分治法最大的區別是,適合於用動態規劃求解的問題,經分解得到子問題往往不是互相獨立的,即下一個子階段的求解是建立在上一個子階段的解的基礎上,進行進一步的求解。若用分治法來解這類問題,則分解得到的子問題數目太多,有些子問題被重複計算了很多次。如果我們能夠儲存已解決的子問題的答案,而在需要時再找出已求得的答案,這樣就可以避免大量的重複計算,節省時間。我們可以用一個表來記錄所有已解的子問題的答案。不管該子問題以後是否被用到,只要它被計算過,就將其結果填入表中。

問題描述:

給定N中物品和一個揹包。物品i的重量是Wi,其價值位Vi ,揹包的容量為C。問應該如何選擇裝入揹包的物品,使得轉入揹包的物品的總價值為最大??

在選擇物品的時候,對每種物品i只有兩種選擇,即裝入揹包或不裝入揹包。不能講物品i裝入多次,也不能只裝入物品的一部分。因此,該問題被稱為0-1揹包問題。

問題分析:令V(i,j)表示在前i(1<=i<=n)個物品中能夠裝入容量為就j(1<=j<=C)的揹包中的物品的最大價值,則可以得到如下的動態規劃函式:

(1)   V(i,0)=V(0,j)=0 

(2)   (a)  V(i,j)=V(i-1,j)    j<wi
(b) V(i,j)=max{V(i-1,j) ,V(i-1,j-wi)+vi) } j>wi

(1)式表明:如果第i個物品的重量大於揹包的容量,則裝人前i個物品得到的最大價值和裝入前i-1個物品得到的最大價是相同的,即物品i不能裝入揹包。 (2)式表明:如果第i個物品的重量小於揹包的容量,則會有一下兩種情況:(a)如果第i個物品沒有裝入揹包,則揹包中物品價值就等於把前i-1個物品裝入容量為j的揹包中所取得的價值。(b)如果把第i個物品裝入揹包,則揹包物品的價值等於第i-1個物品裝入容量位j-wi 的揹包中的價值加上第i個物品的價值vi; 顯然,取二者中價值最大的作為把前i個物品裝入容量為j的揹包中的最優解。

#include <iostream>
#include <algorithm>

using namespace std;

int V[200][200];//前i個物品裝入容量為j的揹包中獲得的最大價值

int KnapSack(int n,int w[],int v[],int x[],int C)
{
    int i,j;
    for(i=0; i<=n; i++)
        V[i][0]=0;
    for(j=0; j<=C; j++)
        V[0][j]=0;
    for(i=0; i<=n-1; i++)
        for(j=0; j<=C; j++)
            if(j<w[i])
                V[i][j]=V[i-1][j];
            else
                V[i][j]=max(V[i-1][j],V[i-1][j-w[i]]+v[i]);
    j=C;
    for(i=n-1; i>=0; i--)
    {
        if(V[i][j]>V[i-1][j])
        {
            x[i]=1;
            j=j-w[i];
        }
        else
            x[i]=0;
    }
    cout << "選中的物品時:" << endl;
    for(i=0; i<n; i++)
        cout << x[i];
    cout << endl;
    return V[n-1][C];

}

int main()
{
    int s;///獲得的最大價值
    int w[15];///物品的重量
    int v[15];///物品的價值
    int x[15];///物品的選取狀態
    int n,i;
    int C;///揹包最大容量
    cout << "請輸入揹包的最大容量:" << endl;
    cin >> C;

    cout << "請輸入物品數:" << endl;
    cin >> n;
    cout << "請分別輸入物品的重量:" << endl;
    for(i=0; i<n; i++)
        cin >> w[i];

    cout << "請分別輸入物品的價值:" << endl;
    for(i=0; i<n; i++)
        cin >> v[i];

    s=KnapSack(n,w,v,x,C);

    cout << "物品價值最大為:" << endl << s << endl;
}