题目链接:http://www.spoj.com/problems/GSS5/

题意:给出一个数列。每次查询最大子段和Sum[i,j],其中i和j满足x1<=i<=y1,x2<=j<=y2,x1<=x2,y1<=y2。

思路:线段树的节点[L,R]保存LMax,RMax,Max,sum,表示左起最大值、右起最大值、区间最大值、区间数字和。更新比较简单。下面说查询。另外设置三个函数,可以查询任意区间[L,R]的最大值,以L开始向右最多到R的最大值、以R开始向左最多到L的最大值,设其分别为Q(L,R),QL(L,R),QR(L,R)。查询分两种情况讨论(设输入数组为d,S[i]表示d的前i项和):

(1)x2>y1: ans=QR(x1,y1)+S[x2-1]-S[y1]+QL(x2,y2);

(2)x2<=y1,那么有x1<=x2<=y1<=y2,则答案为下面三种情况的最大值:

a、x1<=i<=x2<=j<=y2,ans=QR(x1,x2)+QL(x2,y2)-d[x2];

b、x2<=i<=j<=y1,ans=Q(x2,y1);

c、x1<=i<=y1<=j<=y2,ans=QR(x1,y1)+QL(y1,y2)-d[y1]。

struct node
{
    int LMax,RMax,Max,sum;
    int L,R;
};

node a[N<<2];
int d[N],S[N];
int n;

void pushUp(int t)
{
    if(a[t].L==a[t].R) return;

a[t].LMax=max(a[t*2].LMax,a[t*2].sum+a[t*2+1].LMax);
    a[t].RMax=max(a[t*2+1].RMax,a[t*2+1].sum+a[t*2].RMax);
    a[t].sum=a[t*2].sum+a[t*2+1].sum;
    a[t].Max=max3(a[t*2].Max,a[t*2+1].Max,a[t*2].RMax+a[t*2+1].LMax);
}

void build(int t,int L,int R)
{
    a[t].L=L;
    a[t].R=R;
    if(L==R)
    {
        a[t].sum=d[L];
        a[t].LMax=a[t].RMax=a[t].Max=d[L];
        return;
    }
    int mid=(L+R)>>1;
    build(t*2,L,mid);
    build(t*2+1,mid+1,R);
    pushUp(t);
}

int queryL(int t,int L,int R)
{
    if(a[t].L==L&&a[t].R==R) return a[t].LMax;

int mid=(a[t].L+a[t].R)>>1;
    int ans;
    if(R<=mid) ans=queryL(t*2,L,R);
    else if(L>mid) ans=queryL(t*2+1,L,R);
    else
    {
        ans=max(queryL(t*2,L,mid),S[mid]-S[L-1]+queryL(t*2+1,mid+1,R));
    }
    return ans;
}

int queryR(int t,int L,int R)
{
    if(a[t].L==L&&a[t].R==R) return a[t].RMax;

int mid=(a[t].L+a[t].R)>>1;
    int ans;
    if(R<=mid) ans=queryR(t*2,L,R);
    else if(L>mid) ans=queryR(t*2+1,L,R);
    else
    {
        ans=max(queryR(t*2,L,mid)+S[R]-S[mid],queryR(t*2+1,mid+1,R));
    }
    return ans;
}

int query(int t,int L,int R)
{
    if(a[t].L==L&&a[t].R==R) return a[t].Max;

int mid=(a[t].L+a[t].R)>>1;
    int ans;
    if(R<=mid) ans=query(t*2,L,R);
    else if(L>mid) ans=query(t*2+1,L,R);
    else
    {
        int x=query(t*2,L,mid);
        int y=query(t*2+1,mid+1,R);
        ans=max(x,y);
        ans=max(ans,queryR(t*2,L,mid)+queryL(t*2+1,mid+1,R));
    }
    return ans;
}

int cal(int x1,int y1,int x2,int y2)
{
    int a,b,c;
    if(x2>y1)
    {
        a=queryR(1,x1,y1);
        b=queryL(1,x2,y2);
        return a+b+S[x2-1]-S[y1];
    }
    int ans;
    a=queryR(1,x1,x2);
    b=queryL(1,x2,y2);
    ans=a+b-d[x2];

a=query(1,x2,y1);
    ans=max(ans,a);

a=queryR(1,x1,y1);
    b=queryL(1,y1,y2);
    ans=max(ans,a+b-d[y1]);

return ans;
}

int main()
{
    rush()
    {
        RD(n);
        int i;
        FOR1(i,n) RD(d[i]),S[i]=S[i-1]+d[i];
        build(1,1,n);
        int x1,y1,x2,y2;
        int m;
        RD(m);
        while(m--)
        {
            RD(x1,y1); RD(x2,y2);
            PR(cal(x1,y1,x2,y2));
        }
    }
}

SPOJ 2916 Can you answer these queries V(线段树-分类讨论)的更多相关文章

  1. GSS5 spoj 2916. Can you answer these queries V 线段树

    gss5 Can you answer these queries V 给出数列a1...an,询问时给出: Query(x1,y1,x2,y2) = Max { A[i]+A[i+1]+...+A[ ...

  2. SPOJ GSS5 Can you answer these queries V ——线段树

    [题目分析] GSS1上增加区间左右端点的限制. 直接分类讨论就好了. [代码] #include <cstdio> #include <cstring> #include & ...

  3. SPOJ GSS3 Can you answer these queries III[线段树]

    SPOJ - GSS3 Can you answer these queries III Description You are given a sequence A of N (N <= 50 ...

  4. SPOJ 1557. Can you answer these queries II 线段树

    Can you answer these queries II Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 https://www.spoj.com/pr ...

  5. bzoj 2482: [Spoj GSS2] Can you answer these queries II 线段树

    2482: [Spoj1557] Can you answer these queries II Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 145 ...

  6. SPOJ GSS1 - Can you answer these queries I(线段树维护GSS)

    Can you answer these queries I SPOJ - GSS1 You are given a sequence A[1], A[2], -, A[N] . ( |A[i]| ≤ ...

  7. SPOJ GSS1 Can you answer these queries I[线段树]

    Description You are given a sequence A[1], A[2], ..., A[N] . ( |A[i]| ≤ 15007 , 1 ≤ N ≤ 50000 ). A q ...

  8. 【SP2916】Can you answer these queries V - 线段树

    题面 You are given a sequence \(a_1,a_2,...,a_n\). (\(|A[i]| \leq 10000 , 1 \leq N \leq 10000\)). A qu ...

  9. SPOJ GSS3 Can you answer these queries III ——线段树

    [题目分析] GSS1的基础上增加修改操作. 同理线段树即可,多写一个函数就好了. [代码] #include <cstdio> #include <cstring> #inc ...

随机推荐

  1. Python大数据依赖包安装

    一.安装 先安装python2.7.6,win下的numpy这些包需要直接匹配版本,然后安装“numpy-1.8.1-win32-superpack-python2.7”和“scipy-0.16.0- ...

  2. noi2006day2_最大获利 网络流

    这道题是上一题的数据加强版,dinic表示毫无压力: #include<iostream> #include<cstdio> #include<cstring> # ...

  3. [工作积累] Android dynamic library & JNI_OnLoad

    Bionic libc doesn't load dependencies for current .so file (diff from Windows or Linux) so a explici ...

  4. HDU 5014Number Sequence

    思路: 对于一个二进制100011: 尽量将填满:填成111111: 然后有一个很好算的方法 gets(n)表示二进制下N有多少位,N^X=(111111)2 X=111111^N; 其实答案可以直接 ...

  5. 【六】PHP正则表达式方法

    PHP中正则表达式的声明格式有两种方式,一种是POSIX老版模式,已经不常用.还有一种是其他语言中常见的PCRE方法. 1.正则表达式的匹配方法并返回匹配的项:array preg_grep(stri ...

  6. HDU 1397 Goldbach's Conjecture(二分,查找素数)

    题目意思很简单,就是找n=p1+p2的种类数,具体看题目吧. 次重点是查找一定范围内的素数: 重点是用二分查找,不然会超时!!! #include<stdio.h> #include< ...

  7. iOS16进制设置颜色

    UIColor+Hex.h // // UIColor+Hex.h // 16进制颜色类别 // // Created by apple on 15-4-3. // Copyright (c) 201 ...

  8. POJ 1978

    #include <iostream> #define MAXN 55 using namespace std; int _m[MAXN]; int tem[MAXN]; void cop ...

  9. [C++]iostream的几种输入形式

    做ACM题的时候,发现cin并不能满足所有输入要求.比如说: 每行给出一款运动鞋信息,若该款写不打折,则先后给出每种运动鞋单价P,所购买的数量Q:若打折,则先后给出每种运动鞋单价P,所购买的数量Q,折 ...

  10. 01 - 开发成功的Oracle应用

    笔记 1. 开发数据库应用,不能把数据库当黑盒.需要了解数据库的一下内容 数据库的体系结构 并发控制 开发的时候就要调优你的代码 数据库有哪些特性,不要在你的代码里重复实现 深入的学习SQL 2. 我 ...