从零开始的DP学习之肆
当DP方程中的一部分具有某种单调性时可以用数据结构或者预处理维护来降低复杂度
一开始没有看懂题，尴尬，后来发现题目可以简化成这个样子： 将一个序列划分为若干段，每段长度不超过$L$，求每段中最大值之和的最小值
看起来可以直接二分，然而大概并不可行，不过我们有一个明显的$O(n^2)$的DP思路，和我最近做的一道题思路一样。设$dp[i]$表示以$i$为结尾的最小花费，枚举结尾$i$，在$1->i$中找到所有与$i$间距不超过$L$的$j$来转移，代价就是$(j,i)$中的最大值，可以ST表预处理。即$dp[i]=min(dp[i],dp[j]+max(j,i))(j<=i\&\&\sum\limits_{k=j}^i width[k]<=L)$
考虑如何优化，发现$dp$数组在顺序下是单调不下降的，而当右端点$i$确定时，所有的$j$的最大值在顺序下是单调不上升的，所以我们可以考虑用线段树维护这个最小花费，也就是维护$dp[j]+max(j,i)$的最小值。我们枚举右端点时候相当于移动一个长度不超过$L$的滑动窗口，然后我们对于每个新加入的高度就维护一下它左边这段的最大高度，求出$dp[i]$后再单点修改一下下一个位置就可以了
具体来说是用线段树维护四个值+一个标记：最小值，最大值，dp数组，最小花费。其中最小值是为了在维护最大值的时候二分用的，当找到一个要修改的区间的时候我们二分出要修改的那一块，具体来说就是这块中的最小值小于我们要修改成的值再进去修改。然后说说这只蒟蒻都WA了什么鬼畜的错误，第一次是pushup用串了(雾)；第二次是因为调试的时候写的和平时不太一样了，然后线段树release了叶子节点=。=；最后一次发现是区间没卡准(之前都是怎么过的2333)
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=;
long long dp[*N],mini[*N],maxx[*N];
long long laz[*N],val[*N];
long long w[N],h[N];
long long n,m,f,tot,ans;
void pushup(int nde)
{
int ls=*nde,rs=*nde+;
maxx[nde]=max(maxx[ls],maxx[rs]);
mini[nde]=min(mini[ls],mini[rs]);
val[nde]=min(val[ls],val[rs]);
}
void release(int nde)
{
if(laz[nde])
{
int ls=*nde,rs=*nde+;
laz[ls]=mini[ls]=maxx[ls]=laz[nde];
laz[rs]=mini[rs]=maxx[rs]=laz[nde];
val[ls]=dp[ls]+maxx[ls],val[rs]=dp[rs]+maxx[rs]; laz[nde]=;
}
}
void change1(int nde,int l,int r,int nl,int nr,long long task)
{
if(l>nr||r<nl) return ;
if(l!=r) release(nde); int mid=(l+r)/,ls=*nde,rs=*nde+;
if(l>=nl&&r<=nr)
{
if(task<=maxx[nde])
{
if(task>mini[ls]) change1(ls,l,mid,nl,nr,task);
if(task>mini[rs]) change1(rs,mid+,r,nl,nr,task);
}
else
{
maxx[nde]=mini[nde]=laz[nde]=task;
val[nde]=dp[nde]+maxx[nde]; return ;
}
}
else
change1(ls,l,mid,nl,nr,task),change1(rs,mid+,r,nl,nr,task);
pushup(nde);
}
void change2(int nde,int l,int r,int pos,long long task)
{
if(l==pos&&r==pos) {dp[nde]=task; return ;}
release(nde); int mid=(l+r)/,ls=*nde,rs=*nde+;
if(pos<=mid) change2(ls,l,mid,pos,task);
else change2(rs,mid+,r,pos,task);
dp[nde]=min(dp[ls],dp[rs]);
}
long long query(int nde,int l,int r,int nl,int nr)
{
if(l>nr||r<nl)
return 1e18;
else if(l>=nl&&r<=nr)
return val[nde];
else
{
int mid=(l+r)/,ls=*nde,rs=*nde+; release(nde);
return min(query(ls,l,mid,nl,nr),query(rs,mid+,r,nl,nr));
}
}
int main ()
{
scanf("%lld%lld",&n,&m),f=;
for(int i=;i<=n;i++)
scanf("%lld%lld",&h[i],&w[i]);
for(int i=;i<=n;i++)
{
tot+=w[i];
while(tot>m) tot-=w[f++];
change1(,,n,,i,h[i]);
ans=query(,,n,f,i);
if(i!=n) change2(,,n,i+,ans);
}
printf("%lld",ans);
return ;
}