根据lyd 的模板,我们只需要在每次更新信息的时候 同时维护
t[p].sum=t[p*2].sum+t[p*2+1].sum;
t[p].lx=max(t[p*2].lx, t[p*2].sum+t[p*2+1].lx);
t[p].rx=max(t[p*2+1].rx, t[p*2].rx+t[p*2+1].sum);
t[p].dat=Max(t[p*2].dat, t[p*2+1].dat, t[p*2].rx+t[p*2+1].lx);
p节点所在区间的和,p节点所在区间的紧靠左边的最大字段和, 紧靠右边的最大字段和,区间内的最大字段和
每次change单点更新节点的时候 从下往上更新这些值
buil 建树的时候也需要从下往上更新这些值
问题来了区间查询ask (lyd没讲清楚,wa了一下午,最后还是去查了题解,巨菜ε=(´ο`*)))唉)
对于区间l~r 我们如果要查询l~r 的最大连续字段和, 那么肯定是分三种情况,l~r在左子区间,l~r在右子区间
l~r跨越两个区间,前两种情况我们只要递归查询左右子区间就行了。
第三种情况,这时候就要用到我们之前一直维护的那些信息了。
首先我们递归找l~mid 左子区间 让他返回上述那些 信息,然后找mid~r 右子区间 同时也让他返回上面那些信息,
要同时返回这么多信息,肯定是不能直接返回某个值,于是我们用返回结构体(a和b)来保存这些信息,然后
就根据上面的递推式得到
c.sum=a.sum+b.sum;
c.lx=max(a.lx, a.sum+b.lx);
c.rx=max(b.rx, a.rx+b.sum);
c.dat=Max(a.dat, b.dat, a.rx+b.lx);
返回这个c 就是我们想要得到的
刚开始思路 以为可以直接得到最大连续字段和,而第三种情况必须根据上面的信息推导出来,所以每次递归都需要
返回这些信息,最后根据这些信息得出结果。
#include<cstdio>
#include<iostream>
#include<cmath>
#include<queue>
#include<cstring>
#define ll long long
using namespace std;
const int SIZE=500005;
ll a[SIZE];
int n,m;
ll Max(ll a, ll b, ll c)
{
return max(max(a,b) , c);
}
struct segment_tree
{
int l,r;
ll sum;
ll lx;
ll rx;
ll dat;
}t[SIZE*4];
void build(int p, int l, int r)
{
t[p].l=l;
t[p].r=r;
if(l==r)
{
t[p].sum=a[l];
t[p].lx=a[l];
t[p].rx=a[l];
t[p].dat=a[l];
return;
}
int mid=(l+r)>>1;
build(p*2,l,mid);
build(p*2+1,mid+1,r);
t[p].sum=t[p*2].sum+t[p*2+1].sum;
t[p].lx=max(t[p*2].lx, t[p*2].sum+t[p*2+1].lx);
t[p].rx=max(t[p*2+1].rx, t[p*2].rx+t[p*2+1].sum);
t[p].dat=Max(t[p*2].dat, t[p*2+1].dat, t[p*2].rx+t[p*2+1].lx);
}
void change(int p, int x, int v)
{
if(t[p].l==t[p].r)
{
t[p].sum=v;
t[p].lx=v;
t[p].rx=v;
t[p].dat=v;
//cout<<p<<"ass "<<t[p].dat<<endl;
return ;
}
int mid=(t[p].l+t[p].r)>>1;
if(x<=mid)change(p*2, x, v);
else change(p*2+1,x,v);
t[p].sum=t[p*2].sum+t[p*2+1].sum;
t[p].lx=max(t[p*2].lx, t[p*2].sum+t[p*2+1].lx);
t[p].rx=max(t[p*2+1].rx, t[p*2].rx+t[p*2+1].sum);
t[p].dat=Max(t[p*2].dat, t[p*2+1].dat, t[p*2].rx+t[p*2+1].lx);
//cout<<p<<"dd "<<t[p].dat<<endl;
}
segment_tree ask(int p, int l, int r)
{
if(l<=t[p].l&&r>=t[p].r)return t[p];
int mid=(t[p].l+t[p].r)>>1;
if(r<=mid)return ask(p*2,l,r);
else if(l>mid)return ask(p*2+1,l,r);
else
{
segment_tree a=ask(p*2,l,r);
segment_tree b=ask(p*2+1,l,r);
segment_tree c;
c.sum=a.sum+b.sum;
c.lx=max(a.lx, a.sum+b.lx);
c.rx=max(b.rx, a.rx+b.sum);
c.dat=Max(a.dat, b.dat, a.rx+b.lx);
return c;
}
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
build(1,1,n);
int k;
int x,y;
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&k,&x,&y);
if(k==2)
{
change(1,x,y);
}
if(k==1)
{
int temp1=min(x,y);
int temp2=max(x,y);
x=temp1;
y=temp2;
segment_tree ans=ask(1,x,y);
printf("%lld\n",ans.dat);
}
}
return 0;
}
写线段树一定要心平气和!!!理解每次递归所返回的信息能不能对所求的东西有所帮助。