背包问题

总结一下遇到的三种背包问题：0/1背包，完全背包，多重背包。

背包问题(Knapsack problem)是一种组合优化的NP完全问题。问题可以描述为：给定一组物品，每种物品都有自己的重量和价格，在限定的总重量内，我们如何选择，才能使得物品的总价格最高。问题的名称来源于如何选择最合适的物品放置于给定背包中。相似问题经常出现在商业、组合数学，计算复杂性理论、密码学和应用数学等领域中。也可以将背包问题描述为决定性问题，即在总重量不超过W的前提下，总价值是否能达到V？它是在1978年由Merkel和Hellman提出的。

0/1背包

这是最基础的背包问题，特点是：每种物品仅有一件，可以选择放或不放（0不放，1为放，所以称为0/1背包）。
用子问题定义状态：即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]}
对应题目：蒜头君的购物袋2

#include <iostream>
using namespace std;

int items[110];
int values[110];
int dp[110][1010];
int main()
{
    int V, n;
    cin >> V >> n;
    for (int i = 1; i <= n; i++) {
        cin >> items[i];
        cin >> values[i];
    }
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= V; j++) {
            if (j < items[i]) {
                dp[i][j] = dp[i-1][j];
            }else{
                dp[i][j] = max(dp[i-1][j], dp[i-1][j-items[i]] + values[i]);
            }
        }
    }
    cout << dp[n][V];
    return 0;
}

完全背包

有N种物品和一个容量为V的背包，每种物品都有无限件可用。第i种物品的体积是c，价值是w。将哪些物品装入背包可使这些物品的体积总和不超过背包容量，且价值总和最大。
完全的地方体现在，每种物品都有无限件，可以从0/1背包的问题中转移过来。

for (int i = 1; i <= n; ++ i)
    for (int j = c[i]; j <= v; ++ j)
        dp[j] = max(dp[j - c[i]] + w[i], dp[j]);

对应题目：蒜头君的购物袋3

#include<iostream>
using namespace std;
int dp[10001];
int main(){
    int n,V;
    cin>>n>>V;

    int v[1001],p[1001];
    for(int i=1;i<=n;++i)   cin>>v[i]>>p[i];

    for(int i=1;i<=n;++i)
        for(int j=v[i];j<=V;++j)
            dp[j]=max(dp[j-v[i]]+p[i],dp[j]);

    cout<<dp[V];
}

多重背包

有N种物品和一个容量为V的背包。第i种物品最多有n件可用，每件体积是c，价值是w。求解将哪些物品装入背包可使这些物品的体积总和不超过背包容量，且价值总和最大。
这个问题可以转化为0/1背包问题，即将多个相同的物体拆开来看待。但是如果考察多重背包，这样直接拆分一般是会让程序超时，为了避免超时，我们可以采用二进制拆分来优化。具体的意思是，例如有一个物品有13个，可以拆为1,2,4,6，即2^0,2^1,2^2，而2^3>13-2^0+2^1+2^2，故将剩余的6作为一组。这样以来，就可以用这四个值表示从1-13当中任意一个数。
对应题目：平分娃娃

#include <iostream>
using namespace std;

int nums[10];
int dp[1000000];
int vv[500];
int main()
{
    int v = 0;
    for (int i = 1; i <= 6; i++) {
        cin >> nums[i];
        v += nums[i] * i;
    }
    if (v % 2 == 1) {
        cout << "Can't be divided.";
        return 0;
    }else{
        v = v/2;
        int num = 1;
        for (int i = 1; i <= 6; i++) {
            int temp = nums[i], k = 1;
            for (k = 1; k <= temp; k=k*2) {
                vv[num++] = i * k;
                temp -= k;
            }
            if (temp > 0) {
                vv[num++] = temp * i;
            }
        }
        for (int i = 1; i <= num - 1; i++) {
            for (int j = v; j >= vv[i]; j--) {
                dp[j] = max(dp[j], dp[j-vv[i]]+vv[i]);
            }
        }
        if (dp[v] == v) {
            cout << "Can be divided.";
        }else{
            cout << "Can't be divided.";
        }
    }


    return 0;
}