acwing算法提高之动态规划--背包模型（一）

1 基础知识

背包问题分类：

1. 01背包问题：N个物品有体积和价值属性，每个物品只有1个，背包容量V，求背包能装下的最大价值。
2. 完全背包问题：每个物品有无限个。
3. 多重背包问题：每个物品有$x_i$个。
   3.1 二进制优化。
   3.2 单调队列优化（即滑动窗口求最值）。
4. 分组背包问题：每组只能选择一类物品。
5. 混合背包问题。

上述问题是求最大价值，还有下列问题，

1. 求取最大价值时的具体方案。
2. 总价值为$W$时的方案数。
3. 二维费用的背包问题。
4. 有依赖的背包问题。

2 模板

暂无。。。

3 工程化

题目1：多重背包问题III

解题思路：使用单调队列来优化。

C++代码如下，

#include <iostream>
#include <deque>

using namespace std;

const int N = 1010, V = 20010;
int n, m;
int v[N], w[N], s[N];
int f[N][V];
int q[V];

int main() {
    cin >> n >> m;
    for (int i = 1; i <= n; ++i) cin >> v[i] >> w[i] >> s[i];
    
    for (int i = 1; i <= n; ++i) {
        for (int r = 0; r < v[i]; ++r) {
            int hh = 0, tt = -1;
            for (int j = r; j <= m; j += v[i]) {
                if (hh <= tt && q[hh] < j - s[i] * v[i]) ++hh; //队头不在滑动区间内，弹出队头
                f[i][j] = f[i-1][j]; //不选择第i个物品
                if (hh <= tt) f[i][j] = max(f[i][j], f[i-1][q[hh]] + (j - q[hh]) / v[i] * w[i]); //选择第i个物品
                while (hh <= tt && f[i-1][q[tt]] + (j - q[tt]) / v[i] * w[i] <= f[i-1][j]) --tt; //弹出队尾
                q[++tt] = j; //往队尾插入元素
            }
        }
        
        //(2)用deque也会超时
        // for (int r = 0; r < v[i]; ++r) {
        //     deque<int> dq;
        //     for (int j = r; j <= m; j += v[i]) {
        //         if (!dq.empty() && dq.front() < j - s[i] * v[i]) dq.pop_front(); //如果队头不在滑动区间内，则弹出
        //         f[i][j] = f[i-1][j]; //不放第i件物品
        //         if (!dq.empty()) f[i][j] = max(f[i][j], f[i-1][dq.front()] + (j - dq.front()) / v[i] * w[i]); //放第i件物品
        //         while (!dq.empty() && f[i-1][dq.back()] + (j - dq.back()) / v[i] * w[i] <= f[i-1][j]) dq.pop_back();
        //         dq.push_back(j);
        //     }
        // }
        
        //(1)暴力也会超时
        // for (int j = 0; j <= m; ++j) {
        //     f[i][j] = f[i-1][j];
        //     for (int k = 0; k <= s[i]; ++k) {
        //         if (k * v[i] <= j) {
        //             f[i][j] = max(f[i][j], f[i-1][j - k * v[i]] + k * w[i]);
        //         }
        //     }
        // }
    }
    
    cout << f[n][m] << endl;
    
    return 0;
}

同样的，也可以使用二进制表示来优化，但也会超时，C++代码如下，

#include <iostream>

using namespace std;

const int V = 20010, N = 10 * V + 10;
int n, m;
int v[N], w[N], s[N];
int f[V];
int cnt;

int main() {
    cin >> n >> m;
    for (int i = 1; i <= n; ++i) {
        int a, b, c;//a表示体积，b表示价值，c表示数量
        cin >> a >> b >> c;
        
        int k = 1;
        while (c >= k) {
            cnt ++;
            
            v[cnt] = k * a;
            w[cnt] = k * b;
            c -= k;
            
            k <<= 1;
        }
        
        if (c > 0) {
            cnt ++;
            v[cnt] = c * a;
            w[cnt] = c * b;
        }
    }
    
    //转换成01背包问题
    for (int i = 1; i <= cnt; ++i) {
        for (int j = m; j >= v[i]; --j) {
            f[j] = max(f[j], f[j - v[i]] + w[i]);
        }
    }
    
    cout << f[m] << endl;
    
    return 0;
}

题目2：采药。

解题思路：直接应用01背包问题求解即可。C++代码如下，

#include <iostream>

using namespace std;

const int M = 1010;
int f[M];
int m, n;

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

题目3：装箱问题。

解题思路：01背包问题的应用。

#include <iostream>

using namespace std;

const int M = 20010;
int n, m;
int f[M];

int main() {
    cin >> m >> n;
    for (int i = 0; i < n; ++i) {
        int v;
        cin >> v;
        for (int j = m; j >= v; --j) {
            f[j] = max(f[j], f[j - v] + v);
        }
    }
    cout << m - f[m] << endl;
    
    return 0;
}

题目4：宠物小精灵之收服。

解题思路：二维费用的背包问题。

状态定义f[i][j][k]：从前i个物品中选择，总体积1不超过j，总体积2不超过k的选法。其值表示该选法能够获得的最大价值。

状态转移，求以下情况的最大值，

1. 不选择第i个物品，f[i-1][j][k]。
2. 选择第i个物品，f[i-1][j-v1[i]][k-v2[i]] + w[i]。

C++代码如下，

#include <iostream>

using namespace std;

const int N = 110;
const int V1_MAX = 1010, V2_MAX = 510;
int V1, V2, n;
int v1[N], v2[N];
int f[N][V1_MAX][V2_MAX];

int main() {
    cin >> V1 >> V2 >> n;
    //V1表示小智的精灵球数量
    //V2表示皮卡丘的初始体力值
    //n表示野生小精灵的数量
    
    for (int i = 1; i <= n; ++i) {
        cin >> v1[i] >> v2[i];
        //v1[i]表示收服该小精灵需要的精灵球数量
        //v2[i]表示收服过程中对皮卡丘造成的伤害
    }
    
    for (int i = 1; i <= n; ++i) {
        for (int j = 0; j <= V1; ++j) {
            for (int k = 0; k <= V2 - 1; ++k) {
                f[i][j][k] = f[i-1][j][k];
                if (j >= v1[i] && k >= v2[i]) {
                    f[i][j][k] = max(f[i][j][k], f[i-1][j - v1[i]][k - v2[i]] + 1);
                }
            }
        }
    }
    
    cout << f[n][V1][V2-1] << " ";
    int k = V2-1;
    while (k >= 0 && f[n][V1][k] == f[n][V1][V2-1]) k--;
    k += 1;
    cout << V2 - k << endl;
        
    return 0;
}

想想咋优化？可以类比01背包问题，可以将三维状态降阶至二维状态，将“前i个物品”给隐去。

C++代码如下，

#include <iostream>

using namespace std;

const int N = 110;
const int V1_MAX = 1010, V2_MAX = 510;
int V1, V2, n;
int v1[N], v2[N];
int f[V1_MAX][V2_MAX];

int main() {
    cin >> V1 >> V2 >> n;
    //V1表示小智的精灵球数量
    //V2表示皮卡丘的初始体力值
    //n表示野生小精灵的数量
    
    for (int i = 1; i <= n; ++i) {
        cin >> v1[i] >> v2[i];
        //v1[i]表示收服该小精灵需要的精灵球数量
        //v2[i]表示收服过程中对皮卡丘造成的伤害
    }
    
    for (int i = 1; i <= n; ++i) {
        for (int j = V1; j >= v1[i]; --j) {
            for (int k = V2 - 1; k >= v2[i]; --k) { //由于皮卡丘的体力值不能为0，故消耗的总体力值为V2-1
                f[j][k] = max(f[j][k], f[j - v1[i]][k - v2[i]] + 1);
            }
        }
    }
    
    cout << f[V1][V2-1] << " ";
    int k = V2-1;
    while (k >= 0 && f[V1][k] == f[V1][V2-1]) k--;
    k += 1;
    cout << V2 - k << endl;
        
    return 0;
}