3163: [Heoi2013]Eden的新背包问题 多重背包

开始还在想离线把询问分组什么的…后来发现完全不用。
正反两遍多重背包，每次查询找最值就行了。。。
第一次写多重背包，原来一直以为是$O(n^3)$的，原来是$O(n^2logn)$的啊

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<set>
#define N 1005
using namespace std;
int n,Q;
int f[N][N],g[N][N];
struct node
{
    int a,b,c;
}term[N];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
int main()
{
    n=read();
    for (int i=1;i<=n;i++)
        term[i].a=read(),term[i].b=read(),term[i].c=read();
    for (int i=1;i<=n;i++)
    {
        for (int j=0;j<=1000;j++)
            f[i][j]=f[i-1][j];
        int tmp=term[i].c;
        for (int j=1;tmp;j=min(j<<1,tmp))
        {
            for (int k=1000;k>=j*term[i].a;k--)
                f[i][k]=max(f[i][k],f[i][k-j*term[i].a]+j*term[i].b);
            tmp-=j;
        }
    }
    for (int i=n;i;i--)
    {
        for (int j=0;j<=1000;j++)
            g[i][j]=g[i+1][j];
        int tmp=term[i].c;
        for (int j=1;tmp;j=min(j<<1,tmp))
        {
            for (int k=1000;k>=j*term[i].a;k--)
                g[i][k]=max(g[i][k],g[i][k-j*term[i].a]+j*term[i].b);
            tmp-=j;
        }
    }
    Q=read();
    while (Q--)
    {
        int x=read()+1,y=read();
        int ans=0;
        for (int i=0;i<=y;i++)
            ans=max(ans,f[x-1][i]+g[x+1][y-i]);
        printf("%d\n",ans);
    }
    return 0;
}