【bzoj 2653】middle（二分+可持久化线段树）

传送门biu~
先二分答案，对于二分出的答案x，在序列中把大于等于x的值改为1，小于x的值改为-1，那么求区间[a,d]包含[b,c]的最大连续子段和，如果大于等于0说明答案大于等于x，如果小于0说明答案小于x。可以证明经过二分最后得到的答案一定在区间中。
我们可以对每个x序列建一棵线段树维护最大子段和。假设x排好序，每个x的序列跟前一个相比只有一个点的1改为-1，所以可以对线段树进行可持久化。
时间复杂度$O(n\log ^2n)$。

#include<bits/stdc++.h>
#define N 20005
using namespace std;
int n,b[4],T,ans;
struct data{int v,id;}a[N];
inline bool cmp(data a,data b){return a.v<b.v;}
struct Node{
    Node *ch[2];
    int sum,lsum,rsum;
    inline void maintain(){
        sum=ch[0]->sum+ch[1]->sum;
        lsum=max(ch[0]->lsum,ch[0]->sum+ch[1]->lsum);
        rsum=max(ch[1]->rsum,ch[1]->sum+ch[0]->rsum);
    }
}*root[N];
void buildtree(Node *&o,int l,int r){
    o=new Node;
    if(l==r){
        o->sum=o->lsum=o->rsum=1;
        return;
    }
    int mid=l+r>>1;
    buildtree(o->ch[0],l,mid);
    buildtree(o->ch[1],mid+1,r);
    o->maintain();
}
void Modify(Node *&o,Node *u,int l,int r,int x){
    o=new Node;
    if(l==r){
        o->sum=-1;
        o->lsum=o->rsum=0;
        return;
    }
    int mid=l+r>>1;
    if(x<=mid){
        o->ch[1]=u->ch[1];
        Modify(o->ch[0],u->ch[0],l,mid,x);
    }
    else{
        o->ch[0]=u->ch[0];
        Modify(o->ch[1],u->ch[1],mid+1,r,x);
    }
    o->maintain();
}
int Query(Node *o,int l,int r,int L,int R){
    if(l==L && r==R)    return o->sum;
    int mid=l+r>>1;
    if(R<=mid)  return Query(o->ch[0],l,mid,L,R);
    else if(L>mid)  return Query(o->ch[1],mid+1,r,L,R);
    else    return Query(o->ch[0],l,mid,L,mid)+Query(o->ch[1],mid+1,r,mid+1,R);
}
int QueryL(Node *o,int l,int r,int L,int R){
    if(l==L && r==R)    return o->lsum;
    int mid=l+r>>1;
    if(R<=mid)  return QueryL(o->ch[0],l,mid,L,R);
    else if(L>mid)  return QueryL(o->ch[1],mid+1,r,L,R);
    else return max(QueryL(o->ch[0],l,mid,L,mid),Query(o->ch[0],l,mid,L,mid)+QueryL(o->ch[1],mid+1,r,mid+1,R));
}
int QueryR(Node *o,int l,int r,int L,int R){
    if(l==L && r==R)    return o->rsum;
    int mid=l+r>>1;
    if(R<=mid)  return QueryR(o->ch[0],l,mid,L,R);
    else if(L>mid)  return QueryR(o->ch[1],mid+1,r,L,R);
    else return max(QueryR(o->ch[1],mid+1,r,mid+1,R),Query(o->ch[1],mid+1,r,mid+1,R)+QueryR(o->ch[0],l,mid,L,mid));
}
inline bool check(int k){
    return QueryR(root[k],0,n-1,b[0],b[1]-1)+Query(root[k],0,n-1,b[1],b[2])+QueryL(root[k],0,n-1,b[2]+1,b[3])>=0;
}
int main(){
    scanf("%d",&n);
    for(int i=0;i<n;++i)    scanf("%d",&a[i].v),a[i].id=i;
    sort(a,a+n,cmp);
    buildtree(root[0],0,n-1);
    for(int i=1;i<n;++i)    Modify(root[i],root[i-1],0,n-1,a[i-1].id);
    scanf("%d",&T);
    while(T--){
        for(int i=0;i<4;++i)    scanf("%d",&b[i]),b[i]=(b[i]+ans)%n;
        sort(b,b+4);
        int l=0,r=n-1,res;
        while(l<=r){
            int mid=l+r>>1;
            if(check(mid))  res=mid,l=mid+1;
            else    r=mid-1;
        }
        printf("%d\n",ans=a[res].v);
    }
    return 0;
}