POJ2754 Similarity of necklaces 2 多重背包

　　题目链接：http://poj.org/problem?id=2754

　　当时居然没有想到是多重背包，敲了个分组背包，TLE。。。

　　后来才知道可以转化为多重背包，即把每个下界L[i]放为0，上限为U[i]-L[i]，就是物品的个数，价值为P[i]，花费为M[i]。本来题目要求Σ(M[i]*Table[i])=0时,Σ(P[i]*Table[i])的最大值。现在Table[i]的每个值都加了-L[i]，即Σ(M[i]*Table[i] + M[i]*L[i]) 时，Σ(P[i]*Table[i] + P[i]*L[i])的最值，由Σ(M[i]*Table[i])=0可知，背包的容量为Σ(M[i]*L[i])，典型的多重背包问题了。最后再把累加的值减去就行了。

//STATUS:C++_AC_375MS_580KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
#define LL __int64
#define pii pair<int,int>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int N=210,INF=0x3f3f3f3f,MOD=100000000;
const double DNF=100000000000;

int f[N*25*20],pa[N],m[N],l[N],up[N];
int n,rp,rm;

void zo(int v,int w)
{
    int j;
    for(j=rm;j>=v;j--){
        if(f[j-v]>=0)
        f[j]=Max(f[j],f[j-v]+w);
    }
}

int main()
{
 //   freopen("in.txt","r",stdin);
    int i,j,k;
    while(~scanf("%d",&n))
    {
        mem(f,-1);
        for(i=rm=rp=0;i<n;i++){
            scanf("%d%d%d%d",&pa[i],&m[i],&l[i],&up[i]);
            rp+=l[i]*pa[i];
            rm-=l[i]*m[i];
            up[i]-=l[i];
        }
        f[0]=0;
        for(i=0;i<n;i++){
            for(k=1;up[i]>k;up[i]-=k,k<<=1){
                zo(k*m[i],k*pa[i]);
            }
            if(up[i])zo(up[i]*m[i],up[i]*pa[i]);
        }

        printf("%d\n",f[rm]+rp);
    }
    return 0;
}