Description
给出一个序列,序列的每个位置有两个值,H和G。
现在要把这个序列分成几份。每一段必须是连续的,并且要求每一段的max(Hi)的和<=lim,且使每一段的∑Gi的最大值最小。
求这个值。
n<=20000
Solution
看到双最值,直接二分就好~(≧▽≦)/~啦啦啦
关键是如何判断。
我们可以写一个n^2的dp,设fi表示以这一段i结尾的最小的max(Hi)的和,那么转移方程很好写。
Fi=minj符合要求(Fj+maxik=j+1(Hk))
那么这个j的左端点是可以预处理出来的。
我们考虑优化。把这个东西分成两段来用线段树维护。(鄙人蒟蒻一只,不会用stl)
维护每一个区间的min(Fi),和使得min(Fi+max(Hk))最小的Fi和max(Hk)
max还有个常数更小的方法来维护,因为max影响的只是i到它左边第一个比它大的数,预处理就好了。
然后就开打喽。
Code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define N 20005
using namespace std;
struct note{int f,mx,mi;}tr[N*5];
int n,lim,l,r,h[N],g[N],f[N],d[N],pre[N],lazy[N*5];
void mx(int x,int z) {
tr[x].mx=z;lazy[x]=z;tr[x].mi=tr[x].f;
}
void down(int x) {
if (lazy[x]) {
mx(x*2,lazy[x]);
mx(x*2+1,lazy[x]);
lazy[x]=0;
}
}
note merge(note y,note z) {
note x;
x.f=min(y.f,z.f);
if (y.mi+y.mx<z.mi+z.mx) x.mi=y.mi,x.mx=y.mx;
else x.mi=z.mi,x.mx=z.mx;
return x;
}
void change(int v,int l,int r,int x,int y,int z) {
if (l==x&&r==y) {mx(v,z);return;}
int m=(l+r)/2;down(v);
if (y<=m) change(v*2,l,m,x,y,z);
else if (x>m) change(v*2+1,m+1,r,x,y,z);
else change(v*2,l,m,x,m,z),change(v*2+1,m+1,r,m+1,y,z);
tr[v]=merge(tr[v*2],tr[v*2+1]);
}
void modify(int v,int l,int r,int x,int y) {
if (l==r) {tr[v].f=tr[v].mi=y;return;}
int m=(l+r)/2;down(v);
if (x<=m) modify(v*2,l,m,x,y);
else modify(v*2+1,m+1,r,x,y);
tr[v]=merge(tr[v*2],tr[v*2+1]);
}
int find(int v,int l,int r,int x,int y) {
if (l==x&&r==y) return tr[v].mx+tr[v].mi;
int m=(l+r)/2;down(v);
if (y<=m) return find(v*2,l,m,x,y);
else if (x>m) return find(v*2+1,m+1,r,x,y);
else return min(find(v*2,l,m,x,m),find(v*2+1,m+1,r,m+1,y));
}
bool check(int x) {
memset(tr,0,sizeof(tr));int j=1;
fo(i,1,n) {
if (g[i]-g[i-1]>x) return 0;
while (g[i]-g[j-1]>x) j++;
change(1,0,n,pre[i],i,h[i]);
f[i]=find(1,0,n,j-1,i-1);
modify(1,0,n,i,f[i]);
}
if (f[n]<=lim) return 1;else return 0;
}
int main() {
freopen("training5.in","r",stdin);
scanf("%d%d",&n,&lim);l=1;
fo(i,1,n) scanf("%d%d",&h[i],&g[i]),g[i]+=g[i-1];
fo(i,1,n) {
while (d[0]&&h[i]>h[d[d[0]]]) d[0]--;
pre[i]=d[d[0]];d[++d[0]]=i;
}r=g[n];
while (l<r) {
int mid=(l+r)/2;
if (check(mid)) r=mid;else l=mid+1;
}
printf("%d",l);
}