[bzoj3163][HEOI2013]Eden的新背包问题

题目大意

现在有n个物品,第i个物品有c[i]个,每购买第i个物品一个需要a[i]元,可获b[i]代价。
现有m个询问,每次询问形如:第x个物品禁止购买,你有y元的话,你能获得的最大价值是多少?询问之间互相独立。
n<=1000,m<=3*10^5。

多重背包DP

先不考虑禁止购买,这是经典的多重背包问题。
设f[i,j]表示做到第i个物品,目前已经用了j元。
转移显然。
现在我们考虑禁止购买的情况。
我们可以处理前缀DP与后缀DP,然后询问时转化为背包合并问题。容量最大为1000,所以复杂度为1000*m。DP可以使用单调队列优化,但本题没有必要。

参考程序

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<deque>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
int pre[1000+10][1000+10],suf[1000+10][1000+10],a[1000+10],b[1000+10],c[1000+10];
deque<int> dl[100+10];
int i,j,k,l,t,n,m,ans;
int main(){
    scanf("%d",&n);
    fo(i,1,n) scanf("%d%d%d",&a[i],&b[i],&c[i]);
    fo(i,0,c[1]){
        if (a[1]*i>1000) break;
        pre[1][a[1]*i]=b[1]*i;
    }
    fo(i,2,n){
        fo(j,0,100) dl[j].clear();
        fo(j,0,1000){
            k=j%a[i];
            while (!dl[k].empty()){
                l=dl[k].back();
                if (pre[i-1][l]+(j-l)/a[i]*b[i]>pre[i-1][j]) break;
                dl[k].pop_back();
            }
            dl[k].push_back(j);
            while (1){
                l=dl[k].front();
                if ((j-l)/a[i]<=c[i]) break;
                dl[k].pop_front();
            }
            l=dl[k].front();
            pre[i][j]=pre[i-1][l]+(j-l)/a[i]*b[i];
            if (j) pre[i][j]=max(pre[i][j],pre[i][j-1]);
        }
    }
    fo(i,0,c[n]){
        if (a[n]*i>1000) break;
        suf[n][a[n]*i]=b[n]*i;
    }
    fd(i,n-1,1){
        fo(j,0,100) dl[j].clear();
        fo(j,0,1000){
            k=j%a[i];
            while (!dl[k].empty()){
                l=dl[k].back();
                if (suf[i+1][l]+(j-l)/a[i]*b[i]>suf[i+1][j]) break;
                dl[k].pop_back();
            }
            dl[k].push_back(j);
            while (1){
                l=dl[k].front();
                if ((j-l)/a[i]<=c[i]) break;
                dl[k].pop_front();
            }
            l=dl[k].front();
            suf[i][j]=suf[i+1][l]+(j-l)/a[i]*b[i];
            if (j) suf[i][j]=max(suf[i][j],suf[i][j-1]);
        }
    }
    scanf("%d",&m);
    while(m--){
        scanf("%d%d",&j,&k);
        j++;
        ans=0;
        fo(i,0,k) ans=max(ans,pre[j-1][i]+suf[j+1][k-i]);
        printf("%d\n",ans);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值