bzoj 2653 二分答案+可持久化线段树

　　首先离散化，然后我们知道如果对于一个询问的区间[l1,r1],[l2,r2]，我们二分到一个答案x，将[l1,r2]区间中的元素大于等于x的设为1，其余的设为-1，那么如果[l1,r1]的最大右区间和加上[r1,l2]的区间和加上[l2,r2]的最大左区间和大于等于0，那么最大的中位数一定大于等于x。因为这个区间中大于等于x的数量超过了一半，那么我们可以二分答案，然后判断最大的合法（见上文）区间和是否大于等于0。

　　那么对于每个我们二分的值的区间-1,1情况我们不能建立n颗线段树，我们可以建立可持久化线段树来维护这个，最开始的初始值都为1，设rot[x]为二分的值为x的时候区间的1，-1情况的线段树，可以由rot[x-1]这颗线段树继承过来。

　　反思：之前写的可持久化线段树都是建立的权值线段树，这次是用线段树维护区间值的，而且之前对于rot[x]只会有一次插入，这道题的rot[x]可能会有多次插入，在这里纠结了半天。还有我就是最后输出的是adr[ans]，但是前面强制在线的时候我加的是ans,忘了改了= =。

/**************************************************************
    Problem: 2653
    User: BLADEVIL
    Language: C++
    Result: Accepted
    Time:2012 ms
    Memory:18196 kb
****************************************************************/

//By BLADEVIL
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 50010

using namespace std;

struct rec {
    int key,ans,num;
    rec() {
        key=ans=num=;
    }
}a[maxn];

struct segment {
    int left,right,maxr,maxl,sum;
    int son[];
    segment() {
        left=right=maxr=maxl=sum=;
        memset(son,,sizeof son);
    }
}t[*maxn];

int n,m,tot;
int rot[maxn],adr[maxn];

bool cmp1(rec x,rec y) {
    return x.ans<y.ans;
}

bool cmp2(rec x,rec y) {
    return x.num<y.num;
}

void update(int x) {
    t[x].sum=t[t[x].son[]].sum+t[t[x].son[]].sum;
    t[x].maxl=max(t[t[x].son[]].sum+t[t[x].son[]].maxl,t[t[x].son[]].maxl);
    t[x].maxr=max(t[t[x].son[]].sum+t[t[x].son[]].maxr,t[t[x].son[]].maxr);
}

void build(int &x,int l,int r) {
    if (!x) x=++tot;
    t[x].left=l; t[x].right=r;
    if (l==r) {
        t[x].sum=t[x].maxl=t[x].maxr=;
        return ;
    }
    int mid=t[x].left+t[x].right>>;
    build(t[x].son[],l,mid); build(t[x].son[],mid+,r);
    update(x);
}

void insert(int &x,int rot,int y) {
    if (!x) x=++tot;
    t[x].left=t[rot].left; t[x].right=t[rot].right;
    if (t[x].left==t[x].right) {
        t[x].sum=t[x].maxl=t[x].maxr=-;
        return ;
    }
    int mid=t[x].left+t[x].right>>;
    if (y>mid) {
        if (!t[x].son[]) t[x].son[]=t[rot].son[];
        if (t[x].son[]==t[rot].son[]) t[x].son[]=;
        insert(t[x].son[],t[rot].son[],y);
    } else {
        if (!t[x].son[]) t[x].son[]=t[rot].son[];
        if (t[x].son[]==t[rot].son[]) t[x].son[]=;
        insert(t[x].son[],t[rot].son[],y);
    }
    update(x);
}

segment combine(segment x,segment y) {
    segment ans;
    ans.sum=x.sum+y.sum;
    ans.maxl=max(x.sum+y.maxl,x.maxl);
    ans.maxr=max(y.sum+x.maxr,y.maxr);
    return ans;
}

segment query(int x,int l,int r) {
    if ((t[x].left==l)&&(t[x].right==r)) return t[x];
    int mid=t[x].left+t[x].right>>;
    if (l>mid) return query(t[x].son[],l,r); else
    if (r<=mid) return query(t[x].son[],l,r); else
        return combine(query(t[x].son[],l,mid),query(t[x].son[],mid+,r));
}

int main(){
    scanf("%d",&n);
    for (int i=;i<=n;i++) scanf("%d",&a[a[i].num=i].ans);
    sort(a+,a++n,cmp1);
    int sum=,cur=a[].ans; adr[]=a[].ans;
    for (int i=;i<=n;i++) if (a[i].ans==cur) a[i].key=sum; else a[i].key=++sum,adr[sum]=cur=a[i].ans;
    //sort(a+1,a+1+n,cmp2);
    //for (int i=1;i<=n;i++) printf("%d ",a[i].key); printf("\n");
    build(rot[],,n);
    for (int i=;i<=n;i++) insert(rot[a[i].key],rot[a[i].key-],a[i].num);
    //for (int i=1;i<=tot;i++) printf("%d %d %d %d %d %d %d\n",i,t[i].left,t[i].right,t[i].son[0],t[i].son[1],t[i].maxl,t[i].maxr);
    //for (int i=1;i<=n;i++) printf("%d ",adr[i]); printf("\n");
    scanf("%d",&m);
    int ans=;
    while (m--) {
        int ask[]; for (int i=;i<=;i++) scanf("%d",&ask[i]);
        for (int i=;i<=;i++) ask[i]=(ask[i]+adr[ans])%n+;
        ans=;
        sort(ask+,ask+);
        int l=,r=n;
        while (l<=r) {
            int mid=l+r>>,TOT=;
            //printf("%d %d %d\n",l,r,a[mid].key);
            //printf("%d %d\n",l,r);
            segment a1=query(rot[a[mid].key-],ask[],ask[]),a2=query(rot[a[mid].key-],ask[],ask[]);
            if (ask[]+<=ask[]-) TOT=query(rot[a[mid].key-],ask[]+,ask[]-).sum;
            //printf("%d\n",TOT);
            TOT+=a1.maxr+a2.maxl; //printf("%d\n",TOT);
            //printf("%d\n",a1.maxr);
            if (TOT>=) ans=a[mid].key,l=mid+; else r=mid-;
        }
        printf("%d\n",adr[ans]);
    }
    return ;
}