cogs775 山海经 线段树

链接：http://cogs.pro/cogs/problem/problem.php?pid=775

题意：维护区间最大值及其最小字典序来源。

细节巨多……多的狗死人了……

首先我们要建出一棵线段树，这棵线段树要存放以下几个东西：最长区间，起点，终点，最长前缀，前缀终点，最长后缀，后缀起点。（所以维护打了足足六十行……查询打了三十行……）

如何查询呢？最长序列无外乎三种情况：

1、全在左子树中。我们设$len[i]$为$i$节点最长序列，则此时$len[i]=len[i<<1]$；

2、全在右子树中，同理可以得出$len[i]=len[i<<1|1]$。

3、横跨两个子树。设前缀长为$pre[i]$，后缀长为$suf[i]$，那么$len[i]=suf[i<<1]+pre[i<<1|1]$。

维护时顺带按照长度第一字典序第二要求维护区间起终点即可。

查询时按照相同方法进行查询，需要注意的是，最后一次合并时左右可能并未完全更新，上传完成后还要再找一次最优解。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=100005,inf=0x3f3f3f3f;
struct node
{
    int l,r,lmax,rmax,st,en,totmax,pre,suf,sum;
    node()
    {
        lmax=rmax=totmax=-inf;sum=0;
    }
}segtree[maxn<<2];
#define mid ((l+r)>>1)
#define lc root<<1
#define rc root<<1|1
#define lson lc,l,mid
#define rson rc,mid+1,r
void pushup(int root)
{
    segtree[root].sum=segtree[lc].sum+segtree[rc].sum;
    if(segtree[lc].lmax<segtree[lc].sum+segtree[rc].lmax)
    {
        segtree[root].lmax=segtree[lc].sum+segtree[rc].lmax;
        segtree[root].pre=segtree[rc].pre;
    }
    else 
    {
        segtree[root].lmax=segtree[lc].lmax;
        segtree[root].pre=segtree[lc].pre;
    }
    if(segtree[rc].rmax<segtree[rc].sum+segtree[lc].rmax)
    {
        segtree[root].rmax=segtree[rc].sum+segtree[lc].rmax;
        segtree[root].suf=segtree[lc].suf;
    }
    else 
    {
        segtree[root].rmax=segtree[rc].rmax;
        segtree[root].suf=segtree[rc].suf;
    }
    if(segtree[lc].totmax<segtree[rc].totmax)
    {
        if(segtree[rc].totmax<=segtree[lc].rmax+segtree[rc].lmax)
        {
                segtree[root].totmax=segtree[lc].rmax+segtree[rc].lmax;
                segtree[root].st=segtree[lc].suf;
                segtree[root].en=segtree[rc].pre;
        }
        else 
        {
            segtree[root].totmax=segtree[rc].totmax;
            segtree[root].st=segtree[rc].st;
            segtree[root].en=segtree[rc].en;
        }
    }
    else if(segtree[lc].totmax<segtree[lc].rmax+segtree[rc].lmax)
    {
        segtree[root].totmax=segtree[lc].rmax+segtree[rc].lmax;
        segtree[root].st=segtree[lc].suf;
        segtree[root].en=segtree[rc].pre;
    }
    else if(segtree[lc].totmax==segtree[lc].rmax+segtree[rc].lmax)
    {
        segtree[root].totmax=segtree[lc].totmax;
        if(segtree[lc].st>segtree[lc].suf)
        {
            segtree[root].st=segtree[lc].suf;
            segtree[root].en=segtree[rc].pre;
        }
        else 
        {
            segtree[root].st=segtree[lc].st;
            segtree[root].en=segtree[root].en;
        }
    }
    else 
    {
        segtree[root].totmax=segtree[lc].totmax;
        segtree[root].st=segtree[lc].st;
        segtree[root].en=segtree[lc].en;
    }
}
void cmp(node x,node y,node &z)
{
    z.l=x.l,z.r=y.r;z.sum=x.sum+y.sum;
    if(x.lmax<x.sum+y.lmax)
    {
        z.lmax=x.sum+y.lmax;
        z.pre=y.pre;
    }
    else z.lmax=x.lmax,z.pre=x.pre;
    if(y.rmax<y.sum+x.rmax)
    {
        z.rmax=y.sum+x.rmax;
        z.suf=x.suf;
    }
    else 
    {
        z.rmax=y.rmax;
        z.suf=y.suf;
    }
    if(x.totmax<y.totmax)
    {
        if(y.totmax<=x.rmax+y.lmax)
        {
            z.totmax=x.rmax+y.lmax;
            z.st=x.suf;
            z.en=y.pre;
        }
        else 
        {
            z.totmax=y.totmax;
            z.st=y.st;
            z.en=y.en;
        }
    }
    else if(x.totmax<x.rmax+y.lmax)
    {
        z.totmax=x.rmax+y.lmax;
        z.st=x.suf;
        z.en=y.pre;
    }
    else if(x.totmax==x.rmax+y.lmax)
    {
        z.totmax=x.totmax;
        if(x.suf<x.st)z.st=x.suf,z.en=y.pre;
        else z.st=x.st,z.en=x.en;
    }
    else z.totmax=x.totmax,z.st=x.st,z.en=x.en;
}
void build(int root,int l,int r)
{
    segtree[root].l=l,segtree[root].r=r;
    if(l==r)
    {
        int tmp;scanf("%d",&tmp);
        segtree[root].sum=segtree[root].lmax=segtree[root].rmax=segtree[root].totmax=tmp;
        segtree[root].st=segtree[root].en=segtree[root].pre=segtree[root].suf=l;
        return;
    }
    build(lson),build(rson);
    pushup(root);
}
#undef mid
void print(node t)
{
    if(t.lmax>=t.rmax)
        if(t.lmax>=t.totmax)printf("%d %d %d
",t.l,t.pre,t.lmax);
        else printf("%d %d %d
",t.st,t.en,t.totmax);
    else if(t.rmax>t.totmax)printf("%d %d %d
",t.suf,t.r,t.rmax);
    else printf("%d %d %d
",t.st,t.en,t.totmax);
}
node query(int root,int l,int r)
{
    if(segtree[root].l>=l&&segtree[root].r<=r)
        return segtree[root];
    int mid=segtree[root].l+segtree[root].r>>1;node t1,t2,t3;
    if(l<=mid)t1=query(lc,l,r);if(r>mid)t2=query(rc,l,r);
    t1.l=l,t1.r=mid,t2.l=mid+1,t2.r=r;
    cmp(t1,t2,t3);return t3;
}
int n,q;
int haha()
{
    freopen("hill.in","r",stdin);
    freopen("hill.out","w",stdout);
    scanf("%d%d",&n,&q);
    build(1,1,n);
    for(int i=1;i<=q;i++)
    {
        int x,y;scanf("%d%d",&x,&y);
        node t=query(1,x,y);
        print(t);
    }
}
int sb=haha();
int main(){;}

.