【多重背包+二进制优化】HDU 2191 珍惜现在，感恩生活

http://acm.hdu.edu.cn/showproblem.php?pid=2191

悼念512汶川大地震遇难同胞——珍惜现在，感恩生活

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 38069    Accepted Submission(s): 16116

Problem 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

二进制优化：

基本思想为：把同种多件物体转换为多种单件物体。

我们已经知道，1、2、4、8 、16 、 32……2^n  可以组成从1到2^（n+1）-1中的任何数，

所以，对于给定的n，我们只要把它分解为从2^x次幂的数，（最后剩下的一个直接补齐，不要求2^x）就可以组成任意件小于等于x的数且不会超过

对于分分好的这些组，每组当成一个新的物品，都具有新的树形，就成功的把同种多件物体转换成了多种单件物体，就可以用01背包来写，在物品件数较多时优化很大！

本题思路

典型的 多重背包问题，背包模型中的背包容量 为题目中的经费金额， 背包问题中的物品数量，为题目中的大米数量。

使用二进制优化+多重背包

AC CODE

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int nmax=100+10;
int v[nmax];//物品价值 (每袋重量）
int w[nmax];//物品重量 (每袋价格）
int c[nmax];//件数
int dp[nmax];
//最大化v,购买大米的最多重量
 
int value[1111];//分解后的物品价值 
int weight[1111];//分解后的物品重量 


int main(int argc, char** argv) {
	int t;
	int m;//经费金额 (背包容量） 
	int n;//大米种类 (物品数量) 
	while(cin>>t){
		while(t--){
			int cnt=0;
			cin>>m>>n;
			for(int i=0;i<n;i++){
				cin>>w[i]>>v[i]>>c[i];
				
				//对每一种类的c[i]件物品进行二进制分解
				for(int j=1;j<=c[i];j<<=1){//右移=*2 
					value[cnt]=j*v[i];
					weight[cnt]=j*w[i];
					cnt++;
					c[i]-=j; 
				} 
				if(c[i]>0){
					value[cnt]=c[i]*v[i];
					weight[cnt]=c[i]*w[i];
					cnt++;
				}
			}
			//经过上述分解，value为分解后物品价值，weight为分解后物品质量，cnt为n
			//01背包
	        memset(dp,0,sizeof(dp));
			for(int i=0;i<cnt;i++){
				for(int j=m;j>=weight[i];j--){
					dp[j]=max(dp[j],dp[j-weight[i]]+value[i]);
				}
			} 
			cout<<dp[m]<<endl;
		}
	}
	return 0;
}

另一种写法：参自背包九讲

#include <iostream> 
#include <cstring>  
using namespace std;  
  
int weight[110],Value[110],num[110];  
int f[1100];  
int limit;  
  
inline void ZeroOnePack(int w,int v)  
{  
    int j;  
    for(j=limit;j>=w;j--)  
    {  
        if(f[j-w]+v > f[j])  
            f[j]=f[j-w]+v;  
    }  
}  
  
inline void CompletePack(int w,int v)  
{  
    int j;  
    for(j=w;j<=limit;j++)  
    {  
        if(f[j-w]+v > f[j])  
            f[j]=f[j-w]+v;  
    }  
}  
  
inline void MultiplePack(int w,int v,int amount)  
{  
    if(amount * w >= limit)  
    {  
        CompletePack(w,v);  
        return ;  
    }  
    for(int k=1;k<amount;k<<=1)  
    {  
        ZeroOnePack(k*w,k*v);  
        amount -= k;  
    }  
    ZeroOnePack(amount*w,amount*v);  
}  
  
int main()  
{  
    int T,n;  
    cin>>T;  
    while(T--)  
    {  
        cin>>limit>>n;  
          
        for(int i=0;i<n;i++)  
            cin>>weight[i]>>Value[i]>>num[i];  
          
        memset(f,0,sizeof(f));  
          
        for(int i=0;i<n;i++)  
            MultiplePack(weight[i],Value[i],num[i]);  
          
        cout<<f[limit]<<endl;  
    }  
    return 0;  
}