dp学习笔记 4 多重背包(单调队列优化)

题目
题意: 多重背包.
思路: 单调队列优化，O(n*m).

  原来我们常用的写法是对于某个体积j往前找，枚举用多少个物品i, j-vi,j-2vi.
  我们也可以直接往后推，因为对于每个体积j，j%v < v, 枚举0 - v-1为初始体积再枚举物品i的使用次数，也可以推导出所有可放入物品i的相应状态。
这个有什么用呢，对于每个状态 j + xv 可以由于 j + (x-1)v+w, j + (x-2)v+2w,j+(x-k)v+kw 等等推出.
f[j+xv] = max(f[j+kv]+(x-k)w)
f[j+xv] = max(f[j+kv] - kw) + xw
如果我们把kw提出来，这个max只和k相关了。
就变成了一个单调队列问题，把0-v-1为初始值各自进行一次单调队列。
放入单调队列的权值自然是f[j+kv] - kw，持续时间是k + t，因为一个物品最多使用t次。
时间复杂度: O(nm)
代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int q[N][2]; //0:消失时间 1:最大值
int f[N];
void solve()
{
    cin>>n>>m;
    for(int i=1;i<=n;++i)
    {
    	int v,w,t; cin>>v>>w>>t;
    	for(int j=0;j<v;++j) //取模后为0-v-1
    	{
    		int hh = 1,tt = 0;
    		for(int k=j,x=1;k<=m;k+=v,x++)
    		{
    			int a = x+t; //一个物品最多使用t次
    			int b = f[k]-x*w;
    			for(;hh<=tt&&q[tt][1]<=b;--tt);
    			q[++tt][0] = a; q[tt][1] = b;
    			f[k] = q[hh][1] + x*w;
    			for(;hh<=tt&&q[hh][0]==x;++hh);
    		}
    	}
    }
    cout<<f[m];
}
signed main(void)
{  
   T = 1;
   // OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}