codevs5429(单调队列优化多重背包)

最新推荐文章于 2023-08-16 20:13:19 发布

qkoqhh

最新推荐文章于 2023-08-16 20:13:19 发布

阅读量336

点赞数

分类专栏： DP

本文链接：https://blog.csdn.net/qkoqhh/article/details/82822255

版权

DP 专栏收录该内容

116 篇文章 1 订阅

订阅专栏

其实这个知识点并不是很难。。只是觉得并不会这么极限卡log就没去管他。。直到今天差点被卡TAT

稍微推导一下

设d[i][j]为前i种物品体积为j的最大价值

$d[i][j]=max\{d[i-1][j-kv_i]+kw_i\}$

此时 $\small 0\leq k\leq min\{j/v_i,c[i]\}$

令 j=pv_i+q ， k=p-k ，有

$\small d[i][j]=max\{d[i-1][q+(p-k)v_i]+kw_i\}=max\{d[i-1][q+kv_i]-kw_i\}+pw_i$

此时 $\small p-min\{j/v_i,c[i]\}\leq k\leq p$ 即 $\small max\{0,p-c[i]\}\leq k\leq p$

然后考虑到枚举k的时候k一定为整数，所以直接 $\small p-c[i]\leq k\leq p$ 就行。。然后发现这就是划窗。。

所以枚举q后枚举k，然后直接单调队列维护最值即可。。

然后时间复杂度降成了O(nm)而且好像并不比二进制压缩复杂多少。。所以以后就写单调队列了。。

/**
 *　　　　　　　　┏┓　　 　┏┓
 * 　　　　　　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　　　┃　　　　　　　┃ 　
 * 　　　　　　　┃　　　━　　 　┃
 * 　　　　　　　┃　＞　　　＜　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┃...　⌒　... 　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┗━┓　　　┏━┛
 * 　　　　　　　　　┃　　　┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃  　　　　　　
 * 　　　　　　　　　┃　　　┃
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┗━━━┓
 * 　　　　　　　　　┃　　　　　　　┣┓
 * 　　　　　　　　　┃　　　　　　　┏┛
 * 　　　　　　　　　┗┓┓┏━┳┓┏┛
 * 　　　　　　　　　　┃┫┫　┃┫┫
 * 　　　　　　　　　　┗┻┛　┗┻┛
 */ 
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#include<stdlib.h>
#include<assert.h>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,1,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 7005 
#define nm 7005
#define pi 1.1415926535897931
const int inf=1e9+7;
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=0;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
 



int n,m,c[NM],w[NM],b[NM],qh,qt,q[NM];
int d[NM],f[NM];



int main(){
    n=read();m=read();
    inc(i,1,n)c[i]=read(),w[i]=read(),b[i]=read();
    inc(i,1,n){
	inc(j,0,c[i]-1){
	    qh=1;qt=0;
	    for(int k=0;j+k*c[i]<=m;k++){
		while(qh<=qt&&q[qh]<k-b[i])qh++;
		while(qh<=qt&&f[j+q[qt]*c[i]]-q[qt]*w[i]<=f[j+k*c[i]]-k*w[i])qt--;
		q[++qt]=k;
		d[j+k*c[i]]=f[j+q[qh]*c[i]]-q[qh]*w[i]+k*w[i];
	    }
	}
	inc(i,1,m)f[i]=d[i];
    }
    return 0*printf("%d\n",d[m]);
}