HDU 2844 Coins

题目分析

这道题很明显是一道完全背包的题，但是对于只学过01背包的菜菜来说很明心要超时，于是开始学习背包九讲，这里我先摘取一段优化的代码：

另一种好想好写的基本方法是转化为01背包求解：把第i种物品换成n[i]件01背包中的物品，则得到了物品数为Σn[i]的01背包问题，直接求解，复杂度仍然是O(V*Σ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的四件物品。

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

下面给出O(log amount)时间处理一件多重背包中物品的过程，其中amount表示物品的数量。

> procedure MultiplePack(cost,weight,amount)
if cost*amount>=V
    CompletePack(cost,weight)
    return
integer k=1
while k<amount
    ZeroOnePack(k*cost,k*weight)
    amount=amount-k
    k=k*2
ZeroOnePack(amount*cost,amount*weight)

希望你仔细体会这个伪代码，如果不太理解的话，不妨翻译成程序代码以后，单步执行几次，或者头脑加纸笔模拟一下，也许就会慢慢理解了。

初始化问题

如果是第一种问法，要求恰好装满背包，那么在初始化时除了f[0]为0其它f[1..V]均设为-∞，这样就可以保证最终得到的f[N]是一种恰好装满背包的最优解。

为什么呢？可以这样理解：初始化的f数组事实上就是在没有任何物品可以放入背包时的合法状态。如果要求背包恰好装满，那么此时只有容量为0的背包可能被价值为0的nothing“恰好装满”，其它容量的背包均没有合法的解，属于未定义的状态，它们的值就都应该是-∞了。如果背包并非必须被装满，那么 任何容量的背包都有一个合法解“什么都不装”，这个解的价值为0，所以初始时状态的值也就全部为0了。

对于这道题进行判断的时候我采用了另一种只属于本题的方法，就是判断f[i] == i，很明显，只有恰好装满这个背包才是可行的。其余即使能得到结果但是不等于i也是不可行的。当然你也可以采用初始化时除了f[0]其余全部初始化为负无穷的方法，这样只需要判断大于0即可。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int INF = 0x3f3f3f3f;
int f[100050],a[105],c[105];

int MultiplePack(int n,int m)
{
    for(int i = 1; i <= n; i++)
    {
        if(a[i]*c[i] >= m)  //这里用完全背包处理
        {
            for(int j = a[i]; j <= m; j++)
                f[j] = max(f[j], f[j-a[i]]+a[i]);
            continue;
        }
        for(int k = 1; k <= c[i]; k *= 2)  //多重背包优化log n;
        {
            for(int j = m; j >= k*a[i]; j--)
                f[j] = max(f[j], f[j - k*a[i]] + k*a[i]);
            c[i] -= k;
        }
        if(c[i] > 0)
            for(int j = m; j >= c[i]*a[i]; j--)
                f[j] = max(f[j], f[j - c[i]*a[i]] + c[i]*a[i]);
    }
    int ans = 0;
    for(int i = 1; i <= m; i++)
        if(f[i] == i) ans++;
    return ans;
}

int main()
{
    int n,m;
    while(scanf("%d%d", &n, &m) != EOF)
    {
        if(n == 0 && m == 0) break;
        for(int i = 1; i <= n; i++)
            scanf("%d", &a[i]);
        for(int i = 1; i <= n; i++)
            scanf("%d", &c[i]);
        memset(f,0,sizeof(f));
        printf("%d\n", MultiplePack(n, m));
    }
    return 0;
}