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

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

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

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

18 22 100 44 100 2

Sample Output

400

Author

lcy

Source

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

Recommend

lcy

   多重背包问题。研究了一晚上的背包问题。顺便写个模板方便自己以后参考。

#include<algorithm>#include<iostream>#include<string.h>#include<sstream>#include<stdio.h>#include<math.h>#include<vector>#include<string>#include<queue>#include<set>#include<map>//#pragma comment(linker,"/STACK:1024000000,1024000000")using namespace std;const int INF=0x3f3f3f3f;const double eps=1e-8;const double PI=acos(-1.0);const int maxn=100010;typedef __int64 ll;int n,lim,f[100010],p[110],h[110],c[110];void zero(int wei,int val)//01背包求解。wei为费.花val为价值下同.{    for(int i=lim; i>=wei; i--)        f[i]=max(f[i],f[i-wei]+val);}void complete(int wei,int val)//完全背包求解{    for(int i=wei; i<=lim; i++)        f[i]=max(f[i],f[i-wei]+val);}void mul(int wei,int val,int mou)//混合背包求解{    if(wei*mou>=lim)//如果数目充足转化为完全背包    {        complete(wei,val);        return ;    }//二进制优化。多重背包    for(int k=1; k<mou; k<<=1)//必须这样能用lowbit优化。因为这样可以组成任意数。而lowbit不行    {        zero(k*wei,k*val);        mou-=k;    }    zero(mou*wei,mou*val);}int main(){    int i,m,cas;    scanf("%d",&cas);    while(cas--)    {        scanf("%d%d",&n,&m);        lim=n;        for(i=0; i<m; i++)            scanf("%d%d%d",&p[i],&h[i],&c[i]);        memset(f,0,sizeof f);        for(i=0; i<m; i++)//遍历每种选择。          mul(p[i],h[i],c[i]);        printf("%d\n",f[n]);    }    return 0;}

补上O(V*N)的代码：

#include<algorithm>#include<iostream>#include<string.h>#include<sstream>#include<stdio.h>#include<math.h>#include<vector>#include<string>#include<queue>#include<set>#include<map>//#pragma comment(linker,"/STACK:1024000000,1024000000")using namespace std;const int INF=0x3f3f3f3f;const double eps=1e-8;const double PI=acos(-1.0);const int maxn=110;typedef __int64 ll;int n,lim,f[maxn],p[110],h[110],c[110];int q[maxn],id[maxn],head,tail;void zero(int wei,int val)//01背包求解。wei为费.花val为价值下同.{    for(int i=lim; i>=wei; i--)        f[i]=max(f[i],f[i-wei]+val);}void complete(int wei,int val)//完全背包求解{    for(int i=wei; i<=lim; i++)        f[i]=max(f[i],f[i-wei]+val);}void mul(int wei,int val,int mou)//混合背包求解{    if(wei*mou>=lim)//如果数目充足转化为完全背包    {        complete(wei,val);        return ;    }    for(int r=0;r<wei;r++)    {        head=tail=0;        for(int w=r,k=0;w<=lim;w+=wei,k++)        {            int tp=f[w]-k*val;            while(tail>head&&tp>q[tail-1])                tail--;            q[tail]=tp;            id[tail++]=k;            while(k-id[head]>mou)                head++;            f[w]=q[head]+k*val;        }    }}int main(){    int i,m,cas;    scanf("%d",&cas);    while(cas--)    {        scanf("%d%d",&n,&m);        lim=n;        for(i=0; i<m; i++)            scanf("%d%d%d",&p[i],&h[i],&c[i]);        memset(f,0,sizeof f);        for(i=0; i<m; i++)//遍历每种选择。          mul(p[i],h[i],c[i]);        printf("%d\n",f[n]);    }    return 0;}

还是有点差距的。