bzoj4094 [Usaco2013 Dec]Optimal Milking 线段树

Description

---

Farmer John最近购买了N(1 <= N <= 40000)台挤奶机，编号为1 … N，并排成一行。第i台挤奶机每天能够挤M(i)单位的牛奶 (1 < =M(i) <=100,000)。由于机器间距离太近，使得两台相邻的机器不能在同一天使用。Farmer John可以自由选择不同的机器集合在不同的日子进行挤奶。在D(1 < = D < = 50,000)天中，每天Farmer John对某一台挤奶机进行维护，改变该挤奶机的产量。Farmer John希望设计一个挤奶方案，使得挤奶机能够在D天后获取最多的牛奶。

Solution

---

我非常的菜因此没有秒掉这道题，一开始居然忽略数据范围脑补网络流的搞法

开一棵线段树，每个节点记录这段区间分别只要左端点、只要右端点、两边都要、两边都不要的最大值，然后分类上传即可

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

typedef long long LL;
const int N=40005;

LL s1[N<<2],s2[N<<2],s3[N<<2],s4[N<<2];

int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}

void update(LL &x,LL y) {
    x=std:: max(x,y);
}

void modify(int now,int tl,int tr,int x,int v) {
    if (tl==tr) {
        s3[now]=v;
        return ;
    }
    int mid=(tl+tr)>>1;
    if (x<=mid) modify(now<<1,tl,mid,x,v);
    else modify(now<<1|1,mid+1,tr,x,v);
    s1[now]=0;
    update(s1[now],s1[now<<1]+s1[now<<1|1]);
    update(s1[now],s1[now<<1]+s4[now<<1|1]);
    update(s1[now],s3[now<<1]+s4[now<<1|1]);

    s2[now]=0;
    update(s2[now],s2[now<<1]+s2[now<<1|1]);
    update(s2[now],s4[now<<1]+s2[now<<1|1]);
    update(s2[now],s4[now<<1]+s3[now<<1|1]);

    s3[now]=0;
    update(s3[now],s1[now<<1]+s2[now<<1|1]);
    update(s3[now],s1[now<<1]+s3[now<<1|1]);
    update(s3[now],s3[now<<1]+s2[now<<1|1]);

    s4[now]=0;
    update(s4[now],s2[now<<1]+s4[now<<1|1]);
    update(s4[now],s4[now<<1]+s1[now<<1|1]);
    update(s4[now],s4[now<<1]+s4[now<<1|1]);
}

int main(void) {
    int n=read(),T=read(); LL ans=0;
    rep(i,1,n) {
        int x=read();
        modify(1,1,n,i,x);
    }
    for (;T--;) {
        int x=read(),v=read();
        modify(1,1,n,x,v);
        ans+=std:: max(s1[1],std:: max(s2[1],std:: max(s3[1],s4[1])));
    }
    printf("%lld\n", ans);
    return 0;
}