CH4301 Can you answer on these queries III

最新推荐文章于 2022-01-26 18:47:02 发布

strategist_614

最新推荐文章于 2022-01-26 18:47:02 发布

阅读量213

点赞数 1

本文链接：https://blog.csdn.net/strategist_614/article/details/86662787

版权

####

Can you answer on these queries III

题目
CH4301
题意
- 1 操作查询区间 $[x, y]$ 最大连续子段和
- 0 操作单点修改
分析
$0$ 操作单点修改板子
$1$ 操作区间操作，首先分析区间可加性，最大连续子段和满足区间可加性。与普通求区间和不同，要对区间和进行判断，找到最大的和。直接维护区间最大连续子段和是不行的，数据中出现负数，在两个区间合并时，不只是简单的 $t r e e [r t] . d a t = m a x (t r e e [r t < < 1] . d a t, t r e e [r t < < 1 ∣ 1] . d a t)$ ,而是两个区间交界的部分形成的子段可能的和会最大,那么引入两个还需要维护的值，紧靠左端的最大连续和，紧靠右段的最大连续和。现在的答案就是 $t r e e [r t] . d a t = m a x (m a x (t r e e [r t < < 1] . d a t, t r e e [r t < < 1 ∣ 1] . d a t), t r e e [r t < < 1] . r m a x + t r e e [r t < < 1 ∣ 1] . l m a x)$ .想要维护这两个值，同样还是由这个点的左右两个子节点的 $l m a x$ 贡献答案，但是右子节点的 $l m a x$ 要加上左子节点的 $s u m$ (即区间和)，因为区间是连续的。所以还要在维护一个 $s u m$ 区间和.
在执行 $a s k$ 操作时要注意：需要先更新左右子节点区间记为 $a, b$ ，答案记为 $a n s$ ,那么先更新左子节点区间和 $a n s . s u m = a s k (r t < < 1, L, R)$ ,右子节点同理。再累计左右子节点的 $l m a x, r m a x$ ,以 $l a m x$ 为例，当 $L > m i d$ 时，只递归右子节点，答案更新为 $a n s . l m a x = b . l m a x$ ,当 ${\leq} mid$ 左右都会递归，答案更新为 $a n s . l m a x = m a x (a . l m a x, a . s u m + b . l m a x)$ .最后 $a n s . d a t = m a x (m a x (a . d a t, b . d a t), a . r m a x + b . l m a x)$ ,就是最终答案。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 5e5+5;
const int INF = 0x3fffffff;
int a[N];
struct segment_tree
{
  	int l;
  	int r;
    int dat;//区间最大连续子段和
    int lmax;//紧靠左端的最大连续和
    int rmax;//紧靠右段的最大连续和
    int sum;//区间和
}tree[N<<2];
void pushup(int rt)
{
   tree[rt].sum = tree[rt<<1].sum + tree[rt<<1|1].sum;
   tree[rt].lmax = max(tree[rt<<1].lmax,tree[rt<<1].sum+tree[rt<<1|1].lmax);
   tree[rt].rmax = max(tree[rt<<1|1].rmax,tree[rt<<1|1].sum+tree[rt<<1].rmax);
   tree[rt].dat = max(max(tree[rt<<1].dat,tree[rt<<1|1].dat),tree[rt<<1].rmax+tree[rt<<1|1].lmax);
}
void build(int rt,int l,int r)
{
	tree[rt].l = l;
	tree[rt].r = r;
	if(l == r)
	{
		tree[rt].dat = tree[rt].rmax = tree[rt].lmax = tree[rt].sum = a[l];
		return ;
	}
	int mid = (l+r)>>1;
	build(rt<<1,l,mid);
	build(rt<<1|1,mid+1,r);
	pushup(rt);
}
void update(int rt,int x,int val)
{
	int l = tree[rt].l;
	int r = tree[rt].r;
	if(l == r) {tree[rt].dat = tree[rt].sum = tree[rt].lmax = tree[rt].rmax = val;return;}
	int mid = (l+r)>>1;
	if(x <= mid) update(rt<<1,x,val);
	else update(rt<<1|1,x,val);
	pushup(rt);
}
segment_tree ask(int rt,int L,int R)
{
	int l = tree[rt].l;
	int r = tree[rt].r;
	if(L <= l && r <= R) return tree[rt];
	segment_tree a,b,ans;
	a.sum = a.dat = a.lmax = a.rmax = b.sum = b.dat = b.lmax = b.rmax = -INF;
    ans.sum = 0;
	int mid = (l+r)>>1;
	if(L <= mid)
	{
		a = ask(rt<<1,L,R);
		ans.sum += a.sum;
	}
	if(R > mid)
	{
		b = ask(rt<<1|1,L,R);
		ans.sum += b.sum;
	}
	if(L > mid)
		ans.lmax = b.lmax;
	else ans.lmax = max(a.lmax,a.sum+b.lmax);
	if(R <= mid)
		ans.rmax = a.rmax;
	else ans.rmax = max(b.rmax,b.sum+a.rmax);
	ans.dat = max(max(a.dat,b.dat),a.rmax+b.lmax);
	return ans;
}
int main ()
{
    int n,m;
    cin>>n>>m;
    for(int i = 1;i <= n;i++) cin>>a[i];
    build(1,1,n);
    while(m--)
    {
    	int op,x,y;
    	cin>>op>>x>>y;
    	if(op == 1)
    	{
    	   if(x > y) swap(x,y);
           cout<<ask(1,x,y).dat<<endl;
    	}
    	else
    	{
          update(1,x,y);
    	}
    }
	return 0 ;
}