多重背包问题

还是一个体积为m的背包，共有n种物品，其中每种物品有固定的体积，权重和数量 v[i], w[i], s[i], 要求从这i种物品中做出选择，使得在总体积不超过m的条件下，装下的物品有最大的权值。

首先，一种暴力方法就是枚举选择0~s个第i件物品的所有方法，代码如下

#include <iostream>

using namespace std;

const int N = 110; 

int n, m;
int v[N], w[N], s[N];
int f[N][N];

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 j = 0; j <= m; j++)
            for (int k = 0; k <= s[i] && k * v[i] <= j; k++)
                f[i][j] = max(f[i][j], f[i-1][j - k * v[i]] + k * w[i]);
                
    cout << f[n][m] << endl;
    return 0;
}

但是这个时间复杂度是n * m * s， 只能对比较小的数据进行处理，很容易久会tle，因此，我们需要对它进行优化。

但是这个与完全背包不同，不能以同样的方式进行优化（但是完全背包可以用多重背包的方式优化），我们另辟蹊径，于是就有了一种经典的二进制优化方式。

首先我们可以证明，任给整数n，对于最大的整数k, 使得 2^k <= n && 2 ^ (k + 1 ) > n 的条件成立的时候，总会存在常数c（c < 2^(k+1)）使得n 可以被 1, 2, 4, 8, 16, ..., 2^k, c 这列数中的某些元素组合表示。于是，用同样的思路，我们可以将第i种物品分成logs 组，再做一次01背包问题即可。代码如下：

 

#include <iostream>

using namespace std;

const int N = 12000, M = 2010;
// N = 物体种类数量 * logs = 1000 * lg(2000) = 1000 * (1 + lg1000) = 12000
// M = 总体积

int n, m;
int v[N], w[N], s[N];
int f[M];

int main()
{
    cin >> n >> m;
    
    int cnt = 0;//计数变量，用来表示重新分类之后物体的标号，每一类物品视为一个新的整体
    for (int i = 0; i <= n; i++)
    {
        int a, b, s;
        cin >> a >> b >> s;
        int k = 1;
        while(k <= s)
        {
            cnt ++;
            v[cnt] = k * a; // 第cnt类的物体（相当于是从整体来看）的总体积
            w[cnt] = k * b;
            s -= k;
            k *= 2;
        }
        if(s > 0) // 现在的s就等于是常数c
        {
            cnt ++;
            v[cnt] = s * a;
            w[cnt] = s * b;
        }
    }
    n = cnt;一共有cnt个新的物品（每一个物品是之前固定数量物品的总和）
        
    for (int i = 1; i <= n; i++)//简化成01背包问题
        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;
}

再附一道多重背包的小习题：

HDU 2844

这道题相当于是把体积以及权重融在了一起，代码如下：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define N 700000
#define M 110
using namespace std;

int n, m;
int value[M], number[M];
int v[N], f[N];

int main()
{
	while(scanf("%d%d", &n, &m), !(n == 0 && m == 0))
	{
		memset(f, 0, sizeof f);
		for (int i = 1; i <= n; i++)
			scanf("%d", &value[i]);
		for (int i = 1; i <= n; i++)
			scanf("%d", &number[i]);

		int cnt = 0;
		for (int i = 1; i <= n; i++)
		{
			int k = 1;
			while(k <= number[i])
			{
				cnt++;
				v[cnt] = k * value[i];
				number[i] -= k;
				k *= 2;
			}

			if(number[i] > 0)
			{
				cnt++;
				v[cnt] = number[i] * value[i];
			}
		}

		n = cnt;

		for (int i = 1; i <= n; i++)
			for (int j = m; j >= v[i]; j--)
				f[j] = max(f[j], f[j - v[i]] + v[i]);

		int res = 0;
		for (int i = 1; i <= m; i++)
			if(f[i] == i) res ++;
		cout << res << endl;
	}
}