多重背包问题 以及优化

有 N 种物品和一个容量是 V 的背包。

第 i 种物品最多有 si 件，每件体积是 vi，价值是 wi。

求解将哪些物品装入背包，可使物品体积总和不超过背包容量，且价值总和最大。
输出最大价值。

输入格式
第一行两个整数，N，V，用空格隔开，分别表示物品种数和背包容积。

接下来有 N 行，每行三个整数 vi,wi,si，用空格隔开，分别表示第 i 种物品的体积、价值和数量。

输出格式
输出一个整数，表示最大价值。

数据范围
0<N,V≤100
0<vi,wi,si≤100

数据范围很小可以用三重循环：
时间复杂度O(n^2V)

#include<iostream>
#include<algorithm>
using namespace std;

const int N = 110;

int f[N];

int main()
{
    int n, m;
    cin >> n >> m;
    
    for(int i = 0; i < n; i ++)
    {
        int v, w, s;
        cin >> v >> w >> s;
        
        for(int j = m; j >= 0; j--)
        {
            for(int k = 1; k <= s && k * v <= j; k ++)
                f[j] = max(f[j], f[j - v * k] + k * w);
        }
    }
    
    cout << f[m] << endl;
}

数据范围
0<N≤1000
0<V≤2000
0<vi,wi,si≤2000
这样的数据范围就不能用三重循环，可以用二进制优化或者单调队列优化；

二进制优化代码：
时间复杂度：O(n^2logS)

#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;

const int N = 2010;

int f[N];
struct x
{
    int v, w;
};

vector<x> ve;

int main()
{
    int n, m;
    cin >> n >> m;
    
    for(int i = 0; i < n; ++i)
    {
        int v, w, s;
        cin >> v >> w >> s;
        
        for(int j = 1; j <= s; j *= 2)
        {
            s -= j;
            ve.push_back({v * j, w * j});
        }
        if(s)
            ve.push_back({v * s, w * s});
    }
    
    for(auto i : ve)
    {
        for(int j = m; j >= i.v; j --)
            f[j] = max(f[j], f[j - i.v] + i.w);
    }
    
    cout << f[m] << endl;
    
    return 0;
}

单调队列优化代码：
数据范围
0<N≤1000
0<V≤2000
0<vi,wi,si≤2000
时间复杂度o(VN)

#include <iostream>
#include <cstring>

using namespace std;

const int N = 20010;

int dp[N], pre[N], q[N];
int n, m;

int main() {
    cin >> n >> m;
    for (int i = 0; i < n; ++i) {
        memcpy(pre, dp, sizeof(dp));
        int v, w, s;
        cin >> v >> w >> s;
        for (int j = 0; j < v; ++j) {
            int head = 0, tail = -1;
            for (int k = j; k <= m; k += v) {

                if (head <= tail && k - s*v > q[head])
                    ++head;

                while (head <= tail && pre[q[tail]] - (q[tail] - j)/v * w <= pre[k] - (k - j)/v * w)
                    --tail;

                if (head <= tail)
                    dp[k] = max(dp[k], pre[q[head]] + (k - q[head])/v * w);

                q[++tail] = k;
            }
        }
    }
    cout << dp[m] << endl;
    return 0;
}