acwing算法提高之动态规划--背包模型（四）_体积恰好为j的完全背包-CSDN博客

本文链接：https://blog.csdn.net/YMWM_/article/details/134960897

文章介绍了完全背包模型在货币系统中的多个问题求解，如求方案数、最简面值表示、混合背包和有依赖的背包，以及背包问题的最优方案数计算，展示了C++代码实现。

摘要由CSDN通过智能技术生成

1 基础知识

暂无。。。

2 模板

暂无。。。

3 工程化

题目1：货币系统。

解题思路：完全背包模型求方案数。

状态定义f[i][j]：从前i个物品中选体积恰好为j的方案数。
状态转移f[i][j]，以下情况的累加和，

不选择第i个物品，即f[i-1][j]。
选择第i个物品1次，即f[i-1][j-v[i]]。
选择第i个物品2次，即f[i-1][j-2*v[i]]。
……
选择第i个物品max次，即f[i-1][j-max*v[i]]。

状态转移可以简化成，

状态转移f[i][j]，以下情况的累加和，

不选择第i个物品，即f[i-1][j]。
选择第i个物品1次~max次的累加和，为f[i][j-v[i]]。

初始化，f[0][0] = 1。

最终答案为f[n][m]。

同样的，也可以将状态降维，C++代码如下，

#include <iostream>

using namespace std;

const int M = 3010;
int n, m;
long long f[M];

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

题目2：货币系统，求最简的面值表示。

解题思路：完全背包求方案数的应用。

C++代码如下，

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 110, M = 25010;
int n, m;
int v[N];
int f[M];

int main() {
    int T;
    cin >> T;
    
    while (T--) {
        cin >> n;
        for (int i = 0; i < n; ++i) cin >> v[i];
        sort(v, v + n);
        
        int res = 0;
        m = v[n - 1];
        memset(f, 0, sizeof f);
        f[0] = 1;
        for (int i = 0; i < n; ++i) {
            if (f[v[i]] == 0) res += 1;
            for (int j = v[i]; j <= m; ++j) {
                f[j] += f[j - v[i]];
            }
        }
        cout << res << endl;
    }
    
    return 0;
}

题目3：混合背包问题。

解题思路：完全背包归为一类，体积从小到大枚举；01背包和多重背包归为一类，体积从大到小枚举。

C++代码如下，

#include <iostream>

using namespace std;

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

int main() {
    cin >> n >> m;
    for (int i = 0; i < n; ++i) {
        int v, w, s;
        cin >> v >> w >> s;
        
        if (s == 0) {
            //完全背包问题
            for (int j = v; j <= m; ++j) {
                f[j] = max(f[j], f[j - v] + w);
            }
        } else {
            if (s == -1) s = 1;
            for (int k = 1; k <= s; k <<= 1) {
                for (int j = m; j >= k * v; --j) {
                    f[j] = max(f[j], f[j - k * v] + k * w);
                }
                s -= k;
            }
            if (s) {
                for (int j = m; j >= s * v; --j) {
                    f[j] = max(f[j], f[j - s * v] + s * w);
                }
            }
        }
    }
    
    cout << f[m] << endl;
       
    return 0;
}

题目4：有依赖的背包问题。

解题思路：关键在于树结构的dfs。

C++代码如下，

#include <iostream>
#include <cstring>
#include <algorithm>
#include <vector>

using namespace std;

const int N = 110;
int n, m;
int v[N], w[N];
int root;
vector<vector<int>> g; //g[i]表示i号结点的子结点
int f[N][N]; //f[i][j]表示以i号结点为根结点的树且体积不超过j的case下，能获得到的最大价值

void dfs(int u) {
    for (auto son : g[u]) {
        //结点son为结点u的子结点
        dfs(son);
        
        for (int j = m - v[u]; j >= 0; --j) {
            for (int k = 0; k <= j; ++k) {
                f[u][j] = max(f[u][j], f[u][j - k] + f[son][k]);
            }
        }
    }
    
    //将物品u加进去
    for (int i = m; i >= v[u]; --i) f[u][i] = f[u][i - v[u]] + w[u];
    for (int i = 0; i < v[u]; ++i) f[u][i] = 0;
    return;
}

int main() {
    cin >> n >> m;
    g.resize(n + 1);
    for (int i = 1; i <= n; ++i) {
        int p;
        cin >> v[i] >> w[i] >> p;
        //第i号结点的父结点是p
        if (p == -1) {
            root = i;
        } else {
            g[p].emplace_back(i);
        }
    }
    
    dfs(root);
    
    cout << f[root][m] << endl;
    
    return 0;
}

题目5：背包问题求最优方案数。

解题思路：多定义一个g数组，用来存储最优方案数，即g[j]表示体积恰好为j时的最优方案数。而f[j]表示体积恰好为j时的最大价值。初始化很关键，memset(f, -0x3f, sizeof f)，f[0] = 0，g[0] = 1。

C++代码如下，

#include <iostream>
#include <cstring>

using namespace std;

const int N = 1010, mod = 1e9 + 7;
int n, m;
int f[N]; //f[i]表示体积恰好为j的最大价值
int g[N]; //g[i]表示体积恰好为j的最大价值的方案数

int main() {
    cin >> n >> m;
    
    memset(f, -0x3f, sizeof f);
    f[0] = 0;
    g[0] = 1;
    
    for (int i = 0; i < n; ++i) {
        int v, w;
        cin >> v >> w;
        
        for (int j = m; j >= v; --j) {
            int maxv = max(f[j], f[j - v] + w);
            int cnt = 0;
            if (maxv == f[j]) cnt += g[j];
            if (maxv == f[j - v] + w) cnt += g[j - v];
            g[j] = cnt % mod;
            f[j] = maxv;
        }
    }
    
    int res = 0;
    for (int j = 0; j <= m; ++j) res = max(res, f[j]);
    
    int cnt = 0;
    for (int j = 0; j <= m; ++j) {
        if (f[j] == res) {
            cnt = (cnt + g[j]) % mod;
        }
    }
    
    cout << cnt << endl;
    
    return 0;
}

题目6：能量石。

解题思路：第一步贪心，先把顺序确定好。然后就是DP，这里j表示体积恰好等于j，所以状态初始化时需要memset(f, -0x3f, sizeof f)，f[0] = 0。

C++代码如下，

#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>

using namespace std;

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

struct STONE {
    STONE (int a, int b, int c) : s(a), e(b), l(c) {}
    int s;
    int e;
    int l;
};

int main() {
    int T;
    cin >> T;
    
    for (int C = 1; C <= T; ++C) {
        cin >> n;
        vector<STONE> stones;
        m = 0;
        for (int i = 1; i <= n; ++i) {
            int s, e, l;
            cin >> s >> e >> l;
            STONE stone(s, e, l);
            stones.emplace_back(stone);
            m += s;
        }
        
        sort(stones.begin(), stones.end(), [](const STONE stone1, const STONE stone2) {
            return stone1.s * stone2.l < stone2.s * stone1.l;
        });
        
        memset(f, -0x3f, sizeof f);
        f[0] = 0;
        
        for (int i = 0; i < n; ++i) {
            int s = stones[i].s;
            int e = stones[i].e;
            int l = stones[i].l;
            for (int j = m; j >= s; --j) {
                f[j] = max(f[j], f[j - s] + e - (j - s) * l);
            }
        }
        
        int res = 0;
        for (int j = 0; j <= m; ++j) res = max(res, f[j]);
        //Case #1: 105
        cout << "Case #" << C << ": " << res << endl;
    }
    return 0;
}