[bzoj4094][Usaco2013 Dec]Optimal Milking 线段树

4094: [Usaco2013 Dec]Optimal Milking

Time Limit: 10 Sec   Memory Limit: 128 MB
[ Submit][ Status][ Discuss]

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

ta写成tn

多么痛的领悟

#include <bits/stdc++.h>
using namespace std;
const int N = 40000 + 5; long long ans;
int tl[N<<2],tr[N<<2],ta[N<<2],tn[N<<2],n,d,a[N];
void update( int k ){
	tl[k] = max( tl[k<<1] + tn[k<<1|1], max( ta[k<<1] + tn[k<<1|1], tl[k<<1] + tl[k<<1|1] ) );
	tr[k] = max( tn[k<<1] + tr[k<<1|1], max( tn[k<<1] + ta[k<<1|1], tr[k<<1] + tr[k<<1|1] ) );
	tn[k] = max( tn[k<<1] + tn[k<<1|1], max( tn[k<<1] + tl[k<<1|1], tr[k<<1] + tn[k<<1|1] ) );
	ta[k] = max( ta[k<<1] + tr[k<<1|1], max( tl[k<<1] + tr[k<<1|1], tl[k<<1] + ta[k<<1|1] ) );
}
void build( int k, int l, int r ){
	if( l == r ){
		ta[k] = a[l]; return;
	}
	int mid = l + r >> 1;
	build( k<<1, l, mid );
	build( k<<1|1, mid+1, r );
	update( k );
}
void change( int k, int l, int r, int x, int val ){
	if( l == r ){ ta[k] = val; return; }
	int mid = l + r >> 1;
	if( x <= mid ) change( k<<1, l, mid, x, val );
	if( x >  mid ) change( k<<1|1, mid+1, r, x, val );
	update( k );
}
int main(){
	scanf( "%d%d", &n, &d );
	for( int i = 1; i <= n; i++ ) scanf( "%d", &a[i] );
	build( 1, 1, n );
	while( d-- ){
		int i,m;
		scanf( "%d%d", &i, &m );
		change( 1, 1, n, i, m );
		ans += max( max( tl[1], tr[1] ), max( tn[1], ta[1] ) );
	}
	printf( "%lld\n", ans );
	return 0;
}