多重背包

多重背包

一、问题类型

题目描述

有 N 种物品和一个容量是 V 的背包。
第 i 种物品最多有 si 件，每件体积是 vi，价值是 wi。
求解将哪些物品装入背包，可使物品体积总和不超过背包容量，且价值总和最大。 输出最大价值。

输入格式

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

输出格式

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

二、分析

法1

思路：

用三层循环：物品件数，背包容量，选用件数；
在完全背包基础上进行修改，将最后一层循环范围改为：0 <= k <= s[i];
但还要确保背包容量大于物品的体积，所以在枚举中要满足：k * w[i] <= j;

划分：

前i件物品，背包容量为j时，背包所能装下的最大价值；

阶段；

dp[i][j]:前i件物品，背包容量为j时，背包所能装下的最大价值；

转移：

i:物品编号,1~n;
j:背包容量,1~m;
k:选i号物品的件数,0~s[i];
且一定满足k * w[i] <= j;
状态转移方程：
$ dp[i][j] = \max{dp[i][j], dp[i - 1][j - k * w[i]] + k * c[i]} $

法2

思路：

对法1进行滚动数组优化；

划分：

背包容量为j时，背包最大价值；

阶段：

dp[j]:背包容量为j时，背包最大价值；

转移：

:物品编号,1~n;
j:背包容量,1~m;
k:选i号物品的件数,0~s[i];
且一定满足k * w[i] <= j;
状态转移方程：
$ dp[j] = \max{dp[j], dp[j - k * w[i]] + k * c[i]} $

法3

思路:

多重背包本身是01背包，但由于有多个价值相同的物品，所以有更加优化的解法,就是用二进制优化，将相同价值的物品分为 $ 2_0， 2_1，2_2 … $ 个，最后再有一部分不能继续向二进制高阶划分的，将这些被划分的物品分别合成一个物品，那么就得到了价值为1倍的、2倍的、4倍的物品，然后再按照01背包的做法，通过二进制数的加和可以实现物品的所有放入情况。

三、代码

法1

#include <cstdio>
#include <algorithm>
using namespace std;
int n, m, w[105], c[105], s[105], dp[105][1005];
int main() {
	scanf("%d %d", &n, &m);
	for (int i = 1; i <= n; i++) {
		scanf("%d %d %d", &w[i], &c[i], &s[i]);
	}
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			for (int k = 0; k <= s[i] && k * w[i] <= j; k++) {
				dp[i][j] = max(dp[i][j], dp[i - 1][j - k * w[i]] + k * c[i]);
			}
		}
	}
	printf("%d", dp[n][m]); 
	return 0;
}

法2

#include <cstdio>
#include <algorithm>
using namespace std;
int n, m, w[105], c[105], s[105], dp[1005];
int main() {
	scanf("%d %d", &n, &m);
	for (int i = 1; i <= n; i++) {
		scanf("%d %d %d", &w[i], &c[i], &s[i]);
	}
	for (int i = 1; i <= n; i++) {
		for (int j = m; j >= 1; j--) {
			for (int k = 0; k <= s[i] && k * w[i] <= j; k++) {
				dp[j] = max(dp[j], dp[j - k * w[i]] + k * c[i]);
			}
		}
	}
	printf("%d", dp[m]); 
	return 0;
}

法3

#include <cstdio>
#include <algorithm>
using namespace std;
int n, m, w[400005], c[400005], dp[400005], a, b, x, d, n1;
int main() {
	scanf("%d %d", &n, &m);
	for (int i = 1; i <= n; i++) {
		d = 1;
		scanf("%d %d %d", &a, &b, &x);
		while (x >= d) {
			w[++n1] = a * d;
			c[n1] = b * d;
			x -= d;
			d *= 2;
		}
		w[++n1] = a * x;
		c[n1] = b * x;
	}
	for (int i = 1; i <= n1; i++) {
		for (int j = m; j >= w[i]; j--) {
			dp[j] = max(dp[j], dp[j - w[i]] + c[i]);
		}
	}
	printf("%d", dp[m]);
	return 0;
}