HDU-2191（简单背包问题的转化）

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

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2191
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.基本算法：
这题目和完全背包问题很类似。基本的方程只需将完全背包问题的方程略微一改即可，因为对于第i种物品有n[i]+1种策略：取0件，取1件……取 n[i]件。令f[i][v]表示前i种物品恰放入一个容量为v的背包的最大权值，则：f[i][v]=max{f[i-1][v-k*c[i]]+ k*w[i]|0<=k<=n[i]}。复杂度是O(V*∑n[i])。
代码如下：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
using namespace std;
int c[1111],dp[1111][1111];
int w[1111],num[1112];
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        int n,m;
        scanf("%d%d",&m,&n);
        for(int i=1;i<=n;i++)
            scanf("%d%d%d",&w[i],&c[i],&num[i]);
        for(int i=1;i<=n;i++)
        {
            for(int j=m;j>=0;j--)
            {
                dp[i][j]=dp[i-1][j];//进入下个循环之前就要把上一层的值赋给现在的dp[i][j]
                for(int k=0;k<=num[i];k++)
                {
                    if(j>=k*w[i])
                    dp[i][j]=max(dp[i][j],dp[i-1][j-k*w[i]]+k*c[i]);
                }
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}

2.转化为01背包问题

另一种好想好写的基本方法是转化为01背包求解：把第i种物品换成n[i]件01背包中的物品，则得到了物品数为∑n[i]的01背包问题，直接求解，复杂度仍然是O(V*∑n[i])。
代码如下：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
using namespace std;
int c[1111],dp[1111];
int w[1111],num[1112];
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        int n,m;
        scanf("%d%d",&m,&n);
        for(int i=1;i<=n;i++)
            scanf("%d%d%d",&w[i],&c[i],&num[i]);//也可以在这里把物品给分成一件一件的
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=num[i];j++)//此处num数组相当于把一类东西拆分成一件一件的与01背包相同
            {
                for(int k=m;k>=w[i];k--)
                {
                    dp[k]=max(dp[k],dp[k-w[i]]+c[i]);
                }
            }
        }
        printf("%d\n",dp[m]);
    }
    return 0;
}

3.二进制优化的多重背包
我们期望将它转化为01背包问题之后能够像完全背包一样降低复杂度。仍然考虑二进制的思想，我们考虑把第i种物品换成若干件物品，使得原问题中第i种物品可取的每种策略——取0..n[i]件——均能等价于取若干件代换以后的物品。另外，取超过n[i]件的策略必不能出现。

方法是：将第i种物品分成若干件物品，其中每件物品有一个系数，这件物品的费用和价值均是原来的费用和价值乘以这个系数。使这些系数分别为 1,2,4,…,2^(k-1),n[i]-2^k+1，且k是满足n[i]-2^k+1>0的最大整数。例如，如果n[i]为13，就将这种物品分成系数分别为1,2,4,6的四件物品。

分成的这几件物品的系数和为n[i]，表明不可能取多于n[i]件的第i种物品。另外这种方法也能保证对于0..n[i]间的每一个整数，均可以用若干个系数的和表示，这个证明可以分0..2^k-1和2^k..n[i]两段来分别讨论得出，并不难，希望你自己思考尝试一下。

这样就将第i种物品分成了O(log n[i])种物品，将原问题转化为了复杂度为O(V*∑log n[i])的01背包问题，是很大的改进。
代码如下：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
using namespace std;
int c[2111],dp[11111];// 这里的dp数组不知道为什么开的大一点反而用的时间比较短
int w[2111];
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        int n,m;
        int a,b,x,num=0;
        scanf("%d%d",&m,&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d%d",&a,&b,&x);
            int r=1;
            while(r<x)  //这个循环将x件物品通过二进制组成更少的件物品
            {
                w[num]=a*r;
                c[num]=b*r;
                num++;
                x=x-r; //已经取出r件，则要从x中减去
                r=r*2; //二进制进位
            }
            w[num]=a*x; //最后r>x的话直接跳出
            c[num]=b*x;//就像13可以拆成1，2，4，6.其中的6比8小直接跳出
            num++;
        }
        //普通的01背包
        for(int i=0;i<num;i++)
        {
            for(int j=m;j>=w[i];j--)
            {
                    dp[j]=max(dp[j],dp[j-w[i]]+c[i]);
            }
        }
        printf("%d\n",dp[m]);
    }
    return 0;
}