中医药暑假训练赛三 c题 题解 （多重背包问题）

原题：

Description

  主角kirito是使用世界首款完全潜行游戏“刀剑神域（Sword Art Online）”的玩家。曾经很幸运的参与过封闭测试，并买下正式版的kirito，正准备体验游戏的第一次正式营运。但在登入后不久，kirito发现“登出”指令竟然消失，而与此同时自称是SAO游戏设计者“茅场晶彦”的人说：“无法完成攻略就无法离开游戏，只有打倒位于“艾恩葛朗特”顶楼，第100层的头目－达成“完全攻略”才是离开这个世界唯一的方法。并且，在游戏内GAME OVER或是尝试脱下NERvGear，玩家会立刻被NERvGear发出的高频率微波破坏脑部而死亡。”唯有接受这个矛盾事实的人，才能够存活下去。
  自己也被卷入其中的kirito，在游戏的舞台——巨大浮游城堡“艾恩葛朗特”里，以不与人组队的独行剑士身份，逐渐崭露头角，并获得“黑色剑士”的称号。kirito以完全攻略的条件——到达城堡最上层为目标，持续进行严酷且漫长的冒险，在这期间他邂逅了女性细剑使——“闪光”亚丝娜，以及公会“血盟骑士团团长”希兹克利夫，他的命运也一步步产生了巨大的变化。kirito能否从游戏里全身而退……

由于kirito是封弊者，kirito有一个二刀流技能，可以使用星曝气流斩，斩杀了强大的守关BOSS。

但是星曝气流斩需要很庞大的法力值。
现在商店有N个药品，kirito的物品栏有W的容量。
第i个药品有重量w_i,可以恢复法力值v_i，有数量c_i个。
现在请你帮助kirito计算他可以恢复的最大法力值。

Input

第一行两个整数N,W(1 <= N <= 300,1 <= W <= 500000 )
接下来N行，每行三个整数w_i,v_i,c_i(1 <= w_i <= 10000,1 <= v_i <= 10000, 1 <= c_i <= 500)

Output

输出一个整数

Sample Input

3 6

2 2 5

3 3 8

1 4 1

Sample Output

9

 

题目大意：

就是个多重背包问题， 有N种物品和一个容量为V的背包。第i种物品最多有n[i]件可用，每件费用是c[i]，价值是w[i]。求解将哪些物品装入背包可使这些物品的费用总和不超过背包容量，且价值总和最大。

思路：

  1.对于解决多重背包问题，得先知道01背包（每个物品就一件）和完全背包（每个物品无限件）

   https://blog.csdn.net/qq_34374664/article/details/52230368 这里面讲的很详细，很通俗易懂。

  实在看不懂，可以记住套用下面的核心模板。

 01背包模板:     状态转移方程：dp[j] = max(dp[i-1][j],dp[i-1][j - size[i]] + Value[i])

		for (i = 0; i<n; i++)//n表示有几种物品

			for (j = Limit; j >= size[i]; j--)//关键是逆序

                //Limit表示最大容量，size[i]表示第i种物品的重量

					dp[j] = max(dp[j],dp[j - size[i]] + Value[i]);
				
                   //dp[j]表示j容量下的存放的最大价值，Value[i]表示第i种物品的价值。
			

 

完全背包模板：    状态转移方程：f[i][v]=max{f[i-1][v-k*c[i]]+k*w[i]|0<=k*c[i]<=v}

f[0] = 0;//f[0]=0, 关键步骤

		for (int i = 1; i <= n; i++)

			for (int V = c[i]; V <= v; V++)//顺序循环，因为比较的是当前状态

				f[v]=max(f[v],f[v-c[i]]+w[i]);

				

 

2.在了解完以上两个问题后，在对多重背包进行解决。有两种方法。

 一种是把所有物品当做不同的个体，物品总数就由n变成了n+c[1]+c[2]+...+c[i](c[i]为第i种物品的数量），这样就变成了

 01背包问题，可以套用01背包的模板解决。

 但有个问题如果物品数量太多可能会导致运行时间超时，很不幸，这道题目就会超时。

这时需要换个分法，全分成1，1，1，1，1，1......太多了。另一个分法的名字叫做二进制优化，将c[i]个物品分成，1,2,4......份，

剩下不能被2的N次分的算成一整份

二进制优化代码

for (j = 1; j <= c[i]; j <<= 1)
			{
				//<<右移1位，相当于乘2
				Value[count1] = j*v[i];
				size1[count1++] = j*w[i];
				c[i] -= j;
			}
			if (c[i] > 0)
			{
				Value[count1] = c[i] * v[i];
				size1[count1++] = c[i] * w[i];
			}

 

另一种方法就是把多重背包问题当成完全背包问题解决，不过k的范围由（0<=k<=v/w[i])变为(0<=k<=c[i])，对此需要多加一层k的内循环。可能也会导致超时。

AC代码：

#include<cstdio>
#include <iostream>
#include<cstring>
using namespace std;

int v[50008], w[50008], c[50008], dp[500008];
int count1, Value[50008], size1[50008];
int main()
{
	int nCase, Limit, nKind, i, j, k;
	//v[]存价值，w[]存尺寸，c[]存件数
	//在本题中，价值是米的重量，尺寸是米的价格

	//count1存储分解完后的物品总数
	//Value存储分解完后每件物品的价值
	//size1存储分解完后每件物品的尺寸


	count1 = 0;
	while (scanf("%d %d", &nKind, &Limit) != EOF) {
		for (i = 0; i<nKind; i++)
		{
			cin >> w[i] >> v[i] >> c[i];

			//对该种类的c[i]件物品进行二进制分解
			for (j = 1; j <= c[i]; j <<= 1)
			{
				//<<右移1位，相当于乘2
				Value[count1] = j*v[i];
				size1[count1++] = j*w[i];
				c[i] -= j;
			}
			if (c[i] > 0)
			{
				Value[count1] = c[i] * v[i];
				size1[count1++] = c[i] * w[i];
			}
		}

		//经过上面对每一种物品的分解，
		//现在Value[]存的就是分解后的物品价值
		//size1[]存的就是分解后的物品尺寸
		//count就相当于原来的n

		//下面就直接用01背包算法来解
		memset(dp, 0, sizeof(dp));

		for (i = 0; i<count1; i++)
		{

			for (j = Limit; j >= size1[i]; j--)
			{
               //状态转移方程dp[j]=max(dp[j],dp[j-size1[i]]+Value[i])
				if (dp[j] < dp[j - size1[i]] + Value[i])
				{
					dp[j] = dp[j - size1[i]] + Value[i];
				}
				//   cout<<dp[j]<<" ";
			}
			// cout<<endl;
		}

		cout << dp[Limit] << endl;
	}
	return 0;
}