[SP1043] GSS1 - Can you answer these queries I
Posted qixingzhi
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[SP1043] GSS1 - Can you answer these queries I相关的知识,希望对你有一定的参考价值。
传送门:>Here<
题意:求区间最大子段和 $N leq 50000$ 包括多组询问(不需要支持修改)
解题思路
线段树的一道好题
我们可以考虑,如果一组数据全部都是正数,那么问题等同于是查询区间和。然而如果有负数的存在,问题就不一样了
考虑对于每一个节点,维护四个信息:ls(代表当前区间一定顶着左端点的最大子段和),rs(同理,一定顶着右端点的),sum(区间和),val(最大子段和,也就是答案)
考虑进行转移——一个节点的信息由它的两个子节点转移而来
$ls[rt] = Max(ls[rt*2], sum[rt*2] + ls[rt*2+1])$。子段和之所以不包括整段区间是由于右端有负数。因此再往右扩展不会更优
rs同理转移。sum就不说了
$val[rt] = Max{ val[rt*1], val[rt*1+1], ls[rt], rs[rt], rs[rt*2]+ls[rt*2+1] }$. 最难理解的是最后一部分。
想象一下,当前区间的最大子段和要么有一头顶住端点,要么两头都不碰到端点。
对于有一头一定碰到的情况,直接用$ls[rt]和rs[rt]$转移即可。(注意,这里所说的是一定碰到,当然最大子段也有可能碰到,但是不一定)
对于都不碰到的情况,如果其不跨过中间,那么分别用两个子节点的val转移。如果恰好跨过中间,那我们需要把它拼接起来——为了使答案最优,我们考虑拼接$rs[rt*2]和ls[rt*2+1]$ (仔细思考)
查询的时候也一样,还是通过递归来完成转移。这里需要对线段树的query有一个较为深刻的理解——不同于build,query(l,r)表示的是区间$[l, r]$中包含在查询区间的那一部分,而不是真的$[l, r]$。因为在递归的时候我们会判断超界。另外,这里的转移需要刚才的四个参数,因此query的返回值应当是一个结构体,而不是单单一个数值。我暂时还没有想出非结构体的做法……
Code
/*By DennyQi*/ #include <cstdio> #include <queue> #include <cstring> #include <algorithm> #define r read() #define Max(a,b) (((a)>(b)) ? (a) : (b)) #define Min(a,b) (((a)<(b)) ? (a) : (b)) using namespace std; typedef long long ll; const int MAXN = 100010; const int MAXM = 27010; const int INF = 1061109567; inline int read(){ int x = 0; int w = 1; register int c = getchar(); while(c ^ ‘-‘ && (c < ‘0‘ || c > ‘9‘)) c = getchar(); if(c == ‘-‘) w = -1, c = getchar(); while(c >= ‘0‘ && c <= ‘9‘) x = (x << 3) + (x << 1) + c - ‘0‘, c = getchar(); return x * w; } int N,M,x,y,opt,a[MAXN]; struct Data{ int ls,rs,sum,val; }; struct SegmentTree{ int ls[MAXN<<2], rs[MAXN<<2], sum[MAXN<<2], val[MAXN<<2]; inline void Pushup(int rt){ sum[rt] = sum[rt<<1] + sum[rt<<1|1]; ls[rt] = Max(ls[rt<<1], sum[rt<<1]+ls[rt<<1|1]); rs[rt] = Max(rs[rt<<1|1], sum[rt<<1|1]+rs[rt<<1]); val[rt] = Max(Max(val[rt<<1], val[rt<<1|1]), Max(Max(ls[rt], rs[rt]), rs[rt<<1] + ls[rt<<1|1])); } void build(int L, int R, int rt){ if(L >= R){ val[rt] = sum[rt] = ls[rt] = rs[rt] = a[L]; return; } int Mid = (L + R) >> 1; build(L, Mid, rt<<1); build(Mid+1, R, rt<<1|1); Pushup(rt); } Data query(int L, int R, int rt, int x, int y){ if(x<=L && R<=y) return (Data){ls[rt],rs[rt],sum[rt],val[rt]}; int Mid = (L + R) >> 1; if(y <= Mid) return query(L, Mid, rt<<1, x, y); if(x >= Mid+1) return query(Mid+1, R, rt<<1|1, x, y); Data res, t_1 = query(L, Mid, rt<<1, x, y), t_2 = query(Mid+1, R, rt<<1|1, x, y); res.sum = t_1.sum + t_2.sum; res.ls = Max(t_1.ls, t_1.sum + t_2.ls); res.rs = Max(t_2.rs, t_1.rs + t_2.sum); res.val = Max(Max(t_1.val, t_2.val), Max(Max(res.ls, res.rs), t_1.rs + t_2.ls)); return res; } }qxz; int main(){ N=r; for(int i = 1; i <= N; ++i) a[i] = r; qxz.build(1, N, 1); M=r; for(int i = 1; i <= M; ++i){ x = r, y = r; printf("%d ", qxz.query(1, N, 1, x, y).val); } return 0; }
以上是关于[SP1043] GSS1 - Can you answer these queries I的主要内容,如果未能解决你的问题,请参考以下文章
$SP1043$ $GSS1$ $-$ $Can$ $you$ $answer$ $these$ $queries$ $I$
SP1043 GSS1 - Can you answer these queries I 线段树
SP1043 GSS1 - Can you answer these queries I(线段树,区间最大子段和(静态))
luogu SP1043 GSS1 - Can you answer these queries I