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

题目链接

题目大意：现在有n个物品，第i个物品有c[i]个，每购买第i个物品一个需要a[i]元，可获b[i]代价。
有m个询问，每次询问形如：第x个物品禁止购买，你有y元的话，你能获得的最大价值是多少？询问之间互相独立。

题解：CDQ分治
Solve(l,r)表示，当前维护的dp数组，记录的答案是除去[l,r]外的物品的答案
Solve(l,mid)时，用[mid+1,r]内的物品转移dp数组。当l=r时，将所有ban=l=r的询问一起（查dp数组求答案）

我的收获：233

#include <bits/stdc++.h>
using namespace std;

const int N=611111,V=2111;

int n,m;
int t,head[N];
int c[N],limit[N],w[N];
int dp[15][V],ans[N];

struct que{
    int pos,val;
    que(int a=0,int b=0){
        pos=a;val=b;
    }
}q[N];

struct edge{int pos,money,nex;}e[N];

void add(int u,int v,int w){e[t].pos=v,e[t].money=w,e[t].nex=head[u],head[u]=t++;}

void Dynamic_Pro(int *f,int x){
    for(int i=0;i<c[x];i++){
        int l=1,r=1;q[1]=que(0,f[i]);
        for(int j=1;j*c[x]<=(1000-i);j++){
            while(l<=r&&q[r].val<=f[j*c[x]+i]-j*w[x])r--;
            q[++r]=que(j,f[j*c[x]+i]-j*w[x]);
            while(l<=r&&q[l].pos<j-limit[x])l++;
            f[j*c[x]+i]=max(f[j*c[x]+i],q[l].val+j*w[x]);
        }
    }
}

void solve(int l,int r,int d){
    if(l==r){
        for(int i=head[l];~i;i=e[i].nex) ans[e[i].pos]=dp[d-1][e[i].money];
        return ;
    }
    int mid=(l+r)>>1;
    for(int i=0;i<V;i++) dp[d][i]=dp[d-1][i];
    for(int i=mid+1;i<=r;i++) Dynamic_Pro(dp[d],i);
    solve(l,mid,d+1);
    for(int i=0;i<V;i++) dp[d][i]=dp[d-1][i];
    for(int i=l;i<=mid;i++) Dynamic_Pro(dp[d],i);
    solve(mid+1,r,d+1);
}

void work()
{
    scanf("%d",&m);
    for(int i=1;i<=m;i++){
        int b,c;scanf("%d%d",&b,&c);
        add(b+1,i,c);
    }
    solve(1,n,1);
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
}

void init()
{
    memset(head,-1,sizeof(head));
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d%d%d",&c[i],&w[i],&limit[i]);
}

int main()
{
    init();
    work();
    return 0;
}