【NOIP2013模拟联考5】军训(training)

Description

给出一个序列，序列的每个位置有两个值，H和G。
现在要把这个序列分成几份。每一段必须是连续的，并且要求每一段的max(Hi)的和<=lim,且使每一段的∑Gi的最大值最小。
求这个值。
n<=20000

Solution

看到双最值，直接二分就好~(≧▽≦)/~啦啦啦
关键是如何判断。
我们可以写一个n^2的dp，设fi表示以这一段i结尾的最小的max(Hi)的和，那么转移方程很好写。

$$Fi=min_{j符合要求}(Fj+max_{k=j+1}^{i}(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);
}