[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

luogu SP1043 GSS1 - Can you answer these queries I

「SP1043」GSS1 - Can you answer these queries I