BZOJ 4094 [Usaco2013 Dec]Optimal Milking 线段树

Description

Farmer John最近购买了N(1 <= N <= 40000)台挤奶机，编号为1 ... N，并排成一行。第i台挤奶机每天能够挤M(i

)单位的牛奶 (1 < =M(i) <=100,000)。由于机器间距离太近，使得两台相邻的机器不能在同一天使用。Farmer Jo

hn可以自由选择不同的机器集合在不同的日子进行挤奶。在D(1 < = D < = 50,000)天中，每天Farmer John对某一

台挤奶机进行维护，改变该挤奶机的产量。Farmer John希望设计一个挤奶方案，使得挤奶机能够在D天后获取最多

的牛奶。

Input

第1行：两个整数N和D

第2..N+1行：每台挤奶机的M(i)

第N+2..N+D+1行：两个整数i和m，表示每天对机器i进行维护，机器i的产量为m。

Output

最大产量

Sample Input

5 3
1
2
3
4
5
5 2
2 7
1 10

Sample Output

32
【样例解释】
第1天，最优方案为2+4=6 ( 方案1+3+2一样)
第2天，最优方案为7+4=11
第3天，最优方案为10+3+2=15

HINT

传送门

一开始以为间隔地取肯定最大了……事实上很容易举出反例啊= =

然后瞬间不会了2333

其实观察到每个区间，对其它区间的影响只有它的左右端点有没有取。

那么只要弄一棵线段树，维护区间取左、取右、取左右、不取左右的最大值即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int
    N=40005;
int n,D;ll w[N];
struct Segment{ll lrx,lx,rx,x;}seg[N<<2];
void up(int u){
    int l=u<<1,r=u<<1|1;
    seg[u].lrx=max(seg[l].lx+max(seg[r].lrx,seg[r].rx),seg[r].rx+max(seg[l].lrx,seg[l].lx));
    seg[u].lx=max(seg[l].lrx+seg[r].x,seg[l].lx+max(seg[r].x,seg[r].lx));
    seg[u].rx=max(seg[r].lrx+seg[l].x,seg[r].rx+max(seg[l].x,seg[l].rx));
    seg[u].x=max(seg[l].rx+seg[r].x,seg[l].x+seg[r].lx);
}
void build(int id,int l,int r){
    if (l==r){
        seg[id].x=seg[id].lx=seg[id].rx=0;
        seg[id].lrx=w[l];
        return;
    }
    int mid=(l+r)>>1;
    build(id<<1,l,mid);
    build(id<<1|1,mid+1,r);
    up(id);
}
void update(int id,int l,int r,int g,ll num){
    if (l==r){
        seg[id].x=seg[id].lx=seg[id].rx=0;
        seg[id].lrx=num;
        return;
    }
    int mid=(l+r)>>1;
    if (g<=mid) update(id<<1,l,mid,g,num);
        else update(id<<1|1,mid+1,r,g,num);
    up(id);
}
int main(){
    scanf("%d%d",&n,&D);
    for (int i=1;i<=n;i++) scanf("%lld",&w[i]);
    build(1,1,n);
    int x;
    ll y,ans=0;
    while (D--){
        scanf("%d%lld",&x,&y);
        update(1,1,n,x,y);
        ll t=0;
        t=max(seg[1].lrx,seg[1].lx);
        t=max(t,seg[1].rx),t=max(t,seg[1].x);
        ans+=t;
    }
    printf("%lld\n",ans);
    return 0;
}