spoj(GSS3) 1716 Can you answer these queries III(线段树求最大子段和)

http://www.spoj.com/problems/GSS3/

题目大意：给一串长度为N(N<=50000)的序列，有M(M<=50000)条操作。

          操作分两种：

          0 x y :将下标为x的数值改为y

          1 x y :询问[x,y]闭区间的最大连续子段和

          用数组lsum保存当前区间左连续最大子段和，rsum保存右连续最大子段和，sum保存区间子段总和。

          左区间合并：右区间左连续最大和加上左区间总和为左连续区间，右区间合并同理。

          中段连续区间合并：右区间左连续最大加上左区间右连续最大，用msum保存

          询问：从左往右依次询问左区间，中段，以及右区间

代码如下：

#include<iostream>
#include<cstring>
#include<queue>
#include<cmath>
#include<cstdio>
using namespace std;
#define N 1250000
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
int sum[N],rsum[N],msum[N],lsum[N],L,R,ret;
int MAX(int a,int b)
{
	return a>b?a:b;
}
void PushUP(int rt)
{
	sum[rt]=sum[rt<<1]+sum[rt<<1|1];
	lsum[rt]=MAX(lsum[rt<<1],lsum[rt<<1|1]+sum[rt<<1]);//左区间合并取大
	rsum[rt]=MAX(rsum[rt<<1|1],rsum[rt<<1]+sum[rt<<1|1]);//右区间合并取大
	msum[rt]=MAX(MAX(msum[rt<<1],msum[rt<<1|1]),lsum[rt<<1|1]+rsum[rt<<1]);//中断区间合并取大
}
void build(int l,int r,int rt)
{
	if(l==r)
	{
		scanf("%d",&sum[rt]);
		lsum[rt]=msum[rt]=rsum[rt]=sum[rt];
		return;
	}
	int m=(l+r)>>1;
	build(lson);
	build(rson);
	PushUP(rt);
}
void query(int a,int b,int l,int r,int rt)
{
	if(a<=l&&r<=b)
	{
		L=MAX(msum[rt],R+lsum[rt]);//上一区间的右连续最大加上当前区间的左连续最大与中段比大，两者都不能并到下一区间
		R=MAX(rsum[rt],R+sum[rt]);//上一区间的右连续最大加上当前区间的总和与当前区间右连续最大取大，并到下一区间
		ret=MAX(ret,L);//贪心思想，每次取大保存
		ret=MAX(ret,R);
		return;
	}
	int m=(l+r)>>1;
	if(m>=a) query(a,b,lson);
	if(m<b) query(a,b,rson);
}
void update(int a,int b,int l,int r,int rt)
{
	if(l==r)
	{
		sum[rt]=msum[rt]=lsum[rt]=rsum[rt]=b;
		return;
	}
	int m=(l+r)>>1;
	if(m>=a) update(a,b,lson);
	else update(a,b,rson);
	PushUP(rt);
}
int main()
{
    int a,b,c,n,m,i;
	while(~scanf("%d",&n))
	{
		build(1,n,1);
		scanf("%d",&m);
		while(m--)
		{
			scanf("%d%d%d",&a,&b,&c);
			if(a==1)
			{
				L=R=0;
				ret=-99999999;
				query(b,c,1,n,1);
				printf("%d\n",ret);
			}
			else
				update(b,c,1,n,1);
		}
	}
        return 0;
}