SPOJ HKNAP 部分贪心

题意

有S个袋子，袋子容量为Y，合并两个袋子花费为C，雕像有N种。接下来N行每行有一个W和一个V，分别代表雕像的重量和体积。问所带的雕像重量-袋子花费最多可以是多少。

题解

这道题是出现在部分贪心这篇OI论文中的。尽管主题是部分贪心，但实际上就是背包的空间优化。不过不是从算法上进行优化，而是从题目的信息上进行优化。
首先的话，根据题目体积只有1-18这18种情况，那么我们就只保留对于每个体积重量最大的元素就可以了。然后我们求出来性价比最高的元素，作为主要填充物。
我们假设主要填充物的W为wi，V为vi。我们可以分析一下，可以发现其他填充物的数量一定是小于vi的。因为如果其他填充物的数量>=vi，那么一定可以抽出来v[i]*vi的空间，这块空间显而易见要填充性价比最高的元素，因此可以得到其他每种填充物数量一定 < vi。这样的话，我们可以得到填充其他填充物的容量上限lim=vi*18（为了方便，我们把这个上限设置的大一点，允许装vi个主要填充物）。
上述分析是第一步，我们可以根据这一步的结果求出来不填充主要填充物的时候的最大收益。其实就是一个简单的可重复选择元素的背包问题。
求出来这个以后，我们要考虑一下合并的问题。合并的背包数量不超过lim。我们可以假设背包的体积为1（最小值）。如果我们合并之后装的体积>lim，那么就代表一定有一部分空间充满了主要填充物（应为其他物品最多只能装到lim）。既然这一部分装的是主要填充物，那么这一部分和剩余部分是没有必要合并的。
可能有人会有疑问，有没有可能存在不能拆分的情况，比如2个容量3的背包，3个体积2的雕像，如果不合并，就只能装2个，只有合并了才能装3个。我们可以再分析一下，就会发现这个担心是多余的。因为我们组成的背包容量为y*(lim+x)，而lim=vi*18。这样的话，就会发现y*(vi*18+x)，我们很容易就能分离出一个vi*xx的背包，使得总容量为y*(vi*(18-xx)+x) < y*lim。
这样的话，我们就可以计算出合并1-lim种背包的最大收益。这时候需要再进行一次背包计算。这个子问题是，我们有s个背包，我们有lim种背包，问我们能获得的最大收益是多少。进行一下背包计算，输出dp[s]即可。

代码

#include <iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<string>
#include<set>
#include<map>
#include<bitset>
#include<stack>
#include<string>
#define UP(i,l,h) for(int i=l;i<h;i++)
#define DOWN(i,h,l) for(int i=h-1;i>=l;i--)
#define W(a) while(a)
#define MEM(a,b) memset(a,b,sizeof(a))
#define LL long long
#define INF 0x3f3f3f3f3f3f3f3f
#define MAXN 1050
#define MOD 1000000007
#define EPS 1e-3
#define int LL
using namespace std;

int w[20];
int dp[MAXN],ans[MAXN];
int lim,n,s,y,c,wi,vi;

int getW(int x){
    int num=ceil((double)(x*y-lim)/vi);
    return num*wi+dp[x*y-vi*num];
}

main() {
    int t;
    scanf("%lld",&t);
    W(t--) {
        MEM(w,0);
        MEM(dp,0);
        MEM(ans,0);
        scanf("%lld%lld%lld%lld",&n,&s,&y,&c);
        int allv=0;
        UP(i,0,n) {
            int a,b;
            scanf("%lld%lld",&a,&b);
            w[b]=max(w[b],a);
        }
        wi=w[1],vi=1;
        UP(i,2,19) if(w[i]*vi>i*wi) wi=w[i],vi=i;
        lim=vi*18;
        UP(i,1,19) UP(j,0,lim+1) if(j-i>=0) dp[j]=max(dp[j],dp[j-i]+w[i]);
        UP(i,0,lim+1) dp[i]=max(dp[i],dp[i-1]);
        UP(i,1,lim+1) ans[i]=getW(i)-c*(i-1);
        MEM(dp,0);
        UP(i,0,lim+1) UP(j,0,s+1) if(j-i>=0) dp[j]=max(dp[j],dp[j-i]+ans[i]);
        printf("%lld\n",dp[s]);
    }
}