【例题】【线段树】连续区间

1、
NKOJ 2753 区间连续值
时间限制 : 10000 MS 空间限制 : 65536 KB

问题描述
有一数列只有0和1构成，数列中数字个数为为n。
现在有m个形式为x y的提问，询问区间[x,y]中，最多有多少个连续的1。
对于每个询问，请你快速做出回答

输入格式
第一行，两个整数n和m
第二行，n个空格间隔的数字，表示数列
接下来m行，每行两个空格间隔的整数x和y,表示一个询问(x<=y)

输出格式
m行，每行一个整数，对应询问的答案

样例输入
9 3
1 0 0 0 1 1 1 0 1
1 9
3 6
2 4

样例输出
3
2
0

提示
1<=n,m<=100000

#include<iostream>
#include<cstdio>
using namespace std;
const int need=100004;

struct fy{int a,b,fl,fr,v;}t[need*2];
//a,b表边界；fl从左往右可得到的最大；fr从右往左；；v该段内的最大
int le[need*2],ri[need*2],c[need];
int tot=0,x,y;

int max3(int a,int b,int c)
{
    return max(a,max(b,c));
}

void build(int x,int y)
{
    int s=++tot;
    t[s].a=x,t[s].b=y;
    if(x==y)
    {
        t[s].fl=t[s].fr=t[s].mid=t[s].v=c[x];
        return ;
    }
    le[s]=tot+1;build(x,(x+y)/2);
    ri[s]=tot+1;build((x+y)/2+1,y);
    int llen=t[le[s]].b-t[le[s]].a+1,rlen=t[ri[s]].b-t[ri[s]].a+1;
    t[s].fl=t[le[s]].fl==llen?t[le[s]].fl+t[ri[s]].fl:t[le[s]].fl;
    t[s].fr=t[ri[s]].fr==rlen?t[ri[s]].fr+t[le[s]].fr:t[ri[s]].fr;
    //如果左/右儿子全是1，则fl/fr应为左/右儿子长度加上右/左儿子的fl/fr
    t[s].v=max3(t[s].v,t[ri[s]].v,t[le[s]].fr+t[ri[s]].fl);
    //该段内的最大应在左右儿子的最大和中间的最大中选，中间的最大为左儿子右侧加上右儿子左侧
}

int cnt(int s)
{
    if(x>t[s].b||y<t[s].a) return 0;
    if(x<=t[s].a&&t[s].b<=y) return t[s].v;//如果覆盖完了直接返回整个区间的最大值
    int la,ra,ma=0;
    la=cnt(le[s]);
    ra=cnt(ri[s]);
    //如果没覆盖玩应该找左右儿子在区间内的最大值
    if(x<=t[le[s]].b&&t[ri[s]].a<=y)//如果所求区间于左右儿子都有交集，就要找中间的最大值
    {
        int temp=t[le[s]].b-x+1;
        ma+=min(t[le[s]].fr,temp);
        temp=y-t[ri[s]].a+1;
        ma+=min(t[ri[s]].fl,temp);
    }
    return max3(la,ra,ma);//该段内的最大应在左右儿子的最大和中间的最大中选
}

int main()
{
    int n,m;scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&c[i]);
    build(1,n);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        printf("%d\n",cnt(1));
    }
}

2、
NKOJ 1316 小白逛公园
时间限制 : 20000 MS 空间限制 : 65536 KB

问题描述
小新经常陪小白去公园玩，也就是所谓的遛狗啦…在小新家附近有一条“公园路”，路的一边从南到北依次排着n个公园，小白早就看花了眼，自己也不清楚该去哪些公园玩了。
一开始，小白就根据公园的风景给每个公园打了分-.-。小新为了省事，每次遛狗的时候都会事先规定一个范围，小白只可以选择第a个和第b个公园之间（包括a、b两个公园）选择连续的一些公园玩。小白当然希望选出的公园的分数总和尽量高咯。同时，由于一些公园的景观会有所改变，所以，小白的打分也可能会有一些变化。
那么，就请你来帮小白选择公园吧。

输入格式
第一行，两个整数N和M，分别表示表示公园的数量和操作（遛狗或者改变打分）总数。
接下来N行，每行一个整数，依次给出小白 开始时对公园的打分。
接下来M行，每行三个整数。第一个整数K，1或2。K=1表示，小新要带小白出去玩，接下来的两个整数a和b给出了选择公园的范围（1≤a,b≤N）；K=2表示，小白改变了对某个公园的打分，接下来的两个整数p和s，表示小白对第p个公园的打分变成了s（1≤p≤N）。
其中，1≤N≤500 000，1≤M≤100 000，所有打分都是绝对值不超过1000的整数。

输出格式
小白每出去玩一次，都对应输出一行，只包含一个整数，表示小白可以选出的公园得分和的最大值。

样例输入
5 3
1 2 -3 4 5
1 2 3
2 2 -1
1 2 3

样例输出
2
-1

来源 vijos

思路：
该题即为求求任意区间的最大连续和

#include<cstdio>
#include<iostream>
using namespace std;
const int need=500004;
const int inf=1e9;

struct fy{int a,b,fl,fr,fmax,sum;}t[need*2];

int le[need*2],ri[need*2],c[need];
int tot=0,x,y,k,d;

int max3(int a,int b,int c){return max(a,max(b,c));}

void update(int s)
{
    t[s].sum=t[le[s]].sum+t[ri[s]].sum;
    t[s].fl=max(t[le[s]].fl,t[le[s]].sum+t[ri[s]].fl);
    t[s].fr=max(t[ri[s]].fr,t[ri[s]].sum+t[le[s]].fr);
    t[s].fmax=max3(t[le[s]].fmax,t[ri[s]].fmax,t[le[s]].fr+t[ri[s]].fl); 
}

void build(int x,int y)
{
    int s=++tot;
    t[s].a=x,t[s].b=y;
    if(x==y)
    {
        t[s].fl=t[s].fr=t[s].fmax=t[s].sum=c[x];
        return ;
    }
    le[s]=tot+1;build(x,(x+y)>>1);
    ri[s]=tot+1;build((x+y)/2+1,y);
    update(s);
}

void change(int s)
{
    if(t[s].a==t[s].b)
    {
        t[s].fl=t[s].fr=t[s].fmax=t[s].sum=d;
        return ;
    }
    if(t[le[s]].a<=k&&k<=t[le[s]].b) change(le[s]);
    else change(ri[s]);
    update(s);
}
int getr(int s)//查找s向右在x，y区间内可得到的max，注意x<=s.a<=y恒成立 
{ 
    if(x<=t[s].a&&t[s].b<=y) return t[s].fl;//全覆盖 
    int mid=(t[s].a+t[s].b)>>1;
    if(y>mid) return max(t[le[s]].fl,t[le[s]].sum+getr(ri[s])); //覆盖全部s的左儿子和部分右儿子 
    else return getr(le[s]);//覆盖s的部分左儿子 
}
int getl(int s)//查找s向左在x，y区间内可得到的max，注意x<=s.b<=y恒成立
{
    if(x<=t[s].a&&t[s].b<=y) return t[s].fr;
    int mid=(t[s].a+t[s].b)>>1;
    if(mid>=x) return max(t[ri[s]].fr,t[ri[s]].sum+getl(le[s])); //全部右，部分左 
    else return getl(ri[s]);//部分右 
}
int solve(int s)
{
    if(x>t[s].b||y<t[s].a) return -inf;
    if(x<=t[s].a&&t[s].b<=y) return t[s].fmax;
    int la=solve(le[s]),ra=solve(ri[s]),ma=-inf;
    if(x<=t[le[s]].b&&t[ri[s]].a<=y)
    {
        ma=getl(le[s])+getr(ri[s]);
    }
    return max3(la,ra,ma);
}

int main()
{
    int n,m;scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) 
    {
        scanf("%d",&c[i]);
        t[i].fl=t[i].fr=t[i].fmax=-inf;
        t[i+n].fl=t[i+n].fr=t[i+n].fmax=-inf;
    }
    build(1,n);
    for(int i=1,e;i<=m;i++)
    {
        scanf("%d",&e);
        if(e==1)
        {
            scanf("%d%d",&x,&y);
            if(x>y) swap(x,y);//可能出现x>y的情况
            printf("%d\n",solve(1));
        }
        else 
        {
            scanf("%d%d",&k,&d);
            change(1);
        }
    }
} 

3、http://blog.csdn.net/Y__XV/article/details/51913412