题意分析
求一段区间子序列和的最大值,要求在线,操作有:
- 单点值修改
- 区间查询子序列和最大值
如何解决呢
维护四个东西,左边连续最大值,右边连续最大值,区间连续最大值,区间和。
更新操作
- 左边连续最大值,取【左区间左连续最大值,左区间和+右区间左最大值,左区间和+右区间和】中最大值。
- 右边连续最大值,取【右区间右连续最大值,右区间和+左区间右最大值,右区间和+左区间和】中最大值。
- 区间连续最大值,取【左区间连续最大值,右区间连续最大值,左区间右连续最大值+右区间左连续最大值】。
- 区间和左右和求和即可。
查询操作
查询操作有点复杂,首先分类讨论一下。
一个区间的连续最大值无非就以下三种情况,其值等于
- 左区间连续最大值
- 右区间连续最大值
- 左区间右连续最大值+右区间左连续最大值
其实也就是考虑在左边,在右边,和在中间的情况。
在实现的时候,就是按照这个思想,要另外考虑的是,查询区间和维护区间的问题,因为要查询的区间可能不完全包含在节点维护的区间中,所以这里还要分类讨论一下。
- 如果要查询的区间右端点小于等于当前节点的mid,这说明要查询的区间在左区间,向左边查询即可。
- 如果要查询的区间左端点大于当前节点的mid,这说明要查询的区间在右区间,向右边查询即可。
- 剩下的情况就就属于要查询的区间一部分在当前节点的左区间,一部分在右区间,这样的话就要分开查询了。分别向左右查询,然后汇总即可。
对于汇总的答案,直接返回一个结点,该节点就是查询的结果,然后求一下左右查询结果的合并,合并的方法和更新方法类似。
代码总览
#include<bits/stdc++.h>
using namespace std;
const int nmax = 505000;
const int INF = 0x3f3f3f3f;
typedef struct{
int l, r;
int lmax,rmax,tmax,sum;
int mid() {return (l+r) >> 1;}
} Tree;
Tree tree[nmax<<2];
int n,m,num[nmax];
void pushup(int rt){
tree[rt].sum = tree[rt<<1].sum + tree[rt<<1|1].sum;
tree[rt].lmax = max(tree[rt<<1].lmax, max(tree[rt<<1].sum + tree[rt<<1|1].lmax,tree[rt<<1].sum + tree[rt<<1|1].sum));
tree[rt].rmax = max(tree[rt<<1|1].rmax, max(tree[rt<<1|1].sum + tree[rt<<1].rmax,tree[rt<<1].sum + tree[rt<<1|1].sum));
tree[rt].tmax = max(tree[rt<<1].rmax + tree[rt<<1|1].lmax , max(tree[rt<<1].tmax , tree[rt<<1|1].tmax));
}
void build(int l,int r, int rt){
tree[rt].l = l, tree[rt].r = r;
tree[rt].lmax = tree[rt].rmax = tree[rt].sum = tree[rt].tmax = 0;
if(tree[rt].l == tree[rt].r){
tree[rt].lmax = tree[rt].rmax = tree[rt].sum = tree[rt].tmax = num[l];
return;
}
build(l,tree[rt].mid(),rt<<1);
build(tree[rt].mid() + 1, r, rt<<1|1);
pushup(rt);
}
void change(int val ,int pos ,int rt){
if(tree[rt].l == tree[rt].r){
tree[rt].lmax = tree[rt].rmax = tree[rt].sum = tree[rt].tmax = val;
return;
}
if(pos <= tree[rt].mid()) change(val,pos,rt<<1);
else change(val,pos,rt<<1|1);
pushup(rt);
}
Tree query(int l, int r, int rt){
if(tree[rt].l == l && tree[rt].r == r) return tree[rt];
if(r <= tree[rt].mid()) return query(l,r,rt<<1);
else if(l > tree[rt].mid() )return query(l,r,rt<<1|1);
else {
Tree ll = query(l,tree[rt].mid(),rt<<1), rr = query(tree[rt].mid() + 1,r,rt<<1|1);
Tree ret;
ret.sum = ll.sum + rr.sum;
ret.lmax = max(ll.lmax, ll.sum + rr.lmax);
ret.rmax = max(rr.rmax, rr.sum + ll.rmax);
ret.tmax = max(ll.rmax + rr.lmax , max(ll.tmax , rr.tmax));
return ret;
}
}
int main(){
scanf("%d %d",&n,&m);
for(int i = 1;i<=n;++i) scanf("%d",&num[i]);
build(1,n,1);
int k,a,b;
for(int i = 0;i<m;++i){
scanf("%d %d %d",&k,&a,&b);
if(k == 1){
if(a>b) swap(a,b);
Tree tt = query(a,b,1);
printf("%d\n",tt.tmax);
}else{
change(b,a,1);
}
}
return 0;
}