【DP】【单调队列】多重背包

传送门

Description

急！灾区的食物依然短缺！
为了挽救灾区同胞的生命，心系灾区同胞的你准备自己采购一些粮食支援灾区，现在假设你一共有资金n元，而市场有m种大米，每种大米都是袋装产品，其价格不等，并且只能整袋购买。
请问：你用有限的资金最多能采购多少公斤粮食呢？

后记：
人生是一个充满了变数的生命过程，天灾、人祸、病痛是我们生命历程中不可预知的威胁。
月有阴晴圆缺，人有旦夕祸福，未来对于我们而言是一个未知数。那么，我们要做的就应该是珍惜现在，感恩生活——
感谢父母，他们给予我们生命，抚养我们成人；
感谢老师，他们授给我们知识，教我们做人
感谢朋友，他们让我们感受到世界的温暖；
感谢对手，他们令我们不断进取、努力。
同样，我们也要感谢痛苦与艰辛带给我们的财富～

Input
输入数据首先包含一个正整数C，表示有C组测试用例，每组测试用例的第一行是两个整数n和m(1<=n<=100, 1<=m<=100),分别表示经费的金额和大米的种类，然后是m行数据，每行包含3个数p，h和c(1<=p<=20,1<=h<=200,1<=c<=20)，分别表示每袋的价格、每袋的重量以及对应种类大米的袋数。

Output
对于每组测试数据，请输出能够购买大米的最多重量，你可以假设经费买不光所有的大米，并且经费你可以不用完。每个实例的输出占一行。

Sample Input

1
8 2
2 100 4
4 100 2

Sample Output

400

Author
lcy

Source
2008-06-18《 ACM程序设计》期末考试——四川加油！中国加油！

Recommend
lcy | We have carefully selected several similar problems for you: 1114 2602 2159 1171 1059

---

题目大意

你有一个容量为M的背包，和N种物品。
每种物品都有三个属性，vi，wi，与ci，分别表示这种物品的体积、价值和件数。
你的任务是，从这些所给物品中，选出若干件，其体积之和不能超过背包容量，并且使所选物品的权值的和最大。

---

解题思路

先回顾最赤裸裸的多重背包做法。

for (int i=1;i<=n;i++)
    for (int j=m;j>=0;j--)
       for (int k=1;k<=s[i]&&k*w[i]<=j;k++)
          F[j]=min(f[j],f[j-k*w[i]]+v[i]*k);

接下来开始变式

w[i]=物品重量，v[i]=物品价值，num[i]=题目给出的物品数量，c[i]=物品正确数量

单调队列优化的原理是啥呢？朴素状态转移方程：
令 $c[i]=min(num[i],j/w[i]);$
$f[i][j]=max(f[i-1][j-k\ast w[i]]+k\ast v[i])$ (1 <= k <= c[i])
这里的 k 是指取第 i 种物品 k 件。

用物品剩余空间（余数）（就是接下来的$b$） 分层（众所周知，背包的计算是分层的）

–*
友情提示
按照接下来的思路想求出程序中的变式很难，真的转换到疯掉qwq，所以程序中用了另一种方式来理解
如果选择简单理解法，请不要食用下面的解题思路（和程序中的理解方式完全不一样）
–*

如果令 $a=j/w[i],b=j$%$w[i]$那么 $j=a\ast w[i]+b$.
这里用 $k$ 表示的意义改变，$k$ 表示取第 $i$ 种物品的件数比 $a$ 少几件，等于第$i$种物品放弃多少件物品。
那么

f[i][j] = max(f[i-1][b+k*w[i]] - k*v[i]) + a*v[i]
              (a-c[i]<=k<=a)

$b$是其他物品占的空间，$k\ast w[i]$是$i$物品放弃的空间 也是其他物品占的空间
因为后面默认给$i$物品所有价值（$a\ast v[i]$），所以要减去$i$放弃空间的价值（$k\ast v[i]$）

可以发现，$f[i-1][b+k\ast w[i]]-k\ast v[i]$ 在同一层只与 $k$ 有关，而这个 $k$ 是一段连续的。我们要做的就是求出
$f[i-1][b+k\ast w[i]]-k\ast v[i]$ 在 $k$ 取可行区间内时的最大值，这就可以使用单调队列优化。

注意
这是基础思想，程序中的写法是运用了这个思想而已

---

//下面用的是正向思维法，和正确思维不一样，但也说得通
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int T,n,m,w,v,c,head,tail,f[2000],que[2000],s[2000];
int main(){
	scanf("%d",&T);
	while(T--){
    	scanf("%d%d",&m,&n);
    	memset(f,0,sizeof(f));
    	for(int i=1;i<=n;i++){
    		scanf("%d%d%d",&w,&v,&c);
    		if(m/w<c)c=m/w;
    		//求出正确的物品数量
    		for(int j=0;j<w;j++){
    		//枚举余数==分层
    			head=1,tail=1;
    			for(int k=0;k<=(m-j)/w;k++){//枚举拿了的件数
    				int tmp=f[j+k*w]-k*v;//d所贡献的权值
     				while(head<tail&&s[tail-1]<=tmp)tail--;//递减，比当前小，对后面的答案没有贡献
    			    s[tail]=tmp,que[tail++]=k;
     			    while(head<tail&&k-que[head]>c)head++;//超出能拿的范围
     			    if(s[head]+k*v>f[j+k*w])//不知道为什么网站判我max和min有问题，被迫手打
    			       f[j+k*w]=s[head]+k*v;//状态转移
    			}
    		}
    	}
    	printf("%d\n",f[m]);
	}
}