middle（bzoj 2653）

Description

　　一个长度为n的序列a，设其排过序之后为b，其中位数定义为b[n/2]，其中a,b从0开始标号,除法取下整。
　　给你一个长度为n的序列s。
　　回答Q个这样的询问：s的左端点在[a,b]之间,右端点在[c,d]之间的子序列中，最大的中位数。
　　其中a<b<c<d。
　　位置也从0开始标号。
　　我会使用一些方式强制你在线。

Input

　　第一行序列长度n。
　　接下来n行按顺序给出a中的数。
　　接下来一行Q。
　　然后Q行每行a,b,c,d，我们令上个询问的答案是x(如果这是第一个询问则x=0)。
　　令数组q={(a+x)%n,(b+x)%n,(c+x)%n,(d+x)%n}。
　　将q从小到大排序之后，令真正的要询问的a=q[0],b=q[1],c=q[2],d=q[3]。
　　输入保证满足条件。

Output

　　Q行依次给出询问的答案。

Sample Input

5
170337785
271451044
22430280
969056313
206452321
3
3 1 0 2
2 3 1 4
3 1 4 0

271451044
271451044
969056313

Sample Output

 

HINT

　　0:n,Q<=100

　　1,...,5:n<=2000

　　0,...,19:n<=20000,Q<=25000

/*
  二分+主席树 （感觉主席树越来越神奇）
  首先要明白中位数怎么求，可以二分，当我们二分出一个答案时，把区间中大于这个数标为1，小于则标为-1，求和。。。
  对于这道题来说，要是中位数越大越好，如果按照上面的方法，1越多，代表中位数越大，所以二分出答案后，用GSS的方法求RGSS(a,b)+sum(b,c)+LGSS(c,d)，
  判断是否大于0。
  至于怎么快速求出1/-1数组，用主席树维护（？）
*/
#include<cstdio>
#include<iostream>
#include<algorithm>
#define N 20010
using namespace std;
int n,m,cnt,ans;
int a[N],id[N],rt[N],sum[N*],lx[N*],rx[N*],son[N*][];
bool cmp(const int &x1,const int &x2){
    return a[x1]<a[x2];
}
void pushup(int x){
    sum[x]=sum[son[x][]]+sum[son[x][]];
    lx[x]=max(lx[son[x][]],sum[son[x][]]+lx[son[x][]]);
    rx[x]=max(rx[son[x][]],sum[son[x][]]+rx[son[x][]]);
}
void build(int l,int r,int &rt){
    if(l==r){
        rt=++cnt;sum[rt]=lx[rt]=rx[rt]=;
        return;
    }
    rt=++cnt;
    int mid=l+r>>;
    build(l,mid,son[rt][]);
    build(mid+,r,son[rt][]);
    pushup(rt);
}
void insert(int x,int &y,int l,int r,int pos,int val){
    y=++cnt;
    sum[y]=sum[x];lx[y]=lx[x];rx[y]=rx[x];
    son[y][]=son[x][];son[y][]=son[x][];
    if(l==r){
        lx[y]=rx[y]=sum[y]=val;
        return;
    }
    int mid=l+r>>;
    if(pos<=mid) insert(son[x][],son[y][],l,mid,pos,val);
    else insert(son[x][],son[y][],mid+,r,pos,val);
    pushup(y);
}
int get_all(int rt,int l,int r,int x,int y){
    if(l==x&&r==y) return sum[rt];
    int mid=l+r>>;
    if(y<=mid) return get_all(son[rt][],l,mid,x,y);
    else if(x>mid) return get_all(son[rt][],mid+,r,x,y);
    else return get_all(son[rt][],l,mid,x,mid)+get_all(son[rt][],mid+,r,mid+,y);
}
int get_lx(int rt,int l,int r,int x,int y){
    if(l==x&&r==y) return lx[rt];
    int mid=l+r>>;
    if(y<=mid) return get_lx(son[rt][],l,mid,x,y);
    else if(x>mid) return get_lx(son[rt][],mid+,r,x,y);
    else return max(get_lx(son[rt][],l,mid,x,mid),get_all(son[rt][],l,mid,x,mid)+get_lx(son[rt][],mid+,r,mid+,y));
}
int get_rx(int rt,int l,int r,int x,int y){
    if(l==x&&r==y) return rx[rt];
    int mid=l+r>>;
    if(y<=mid) return get_rx(son[rt][],l,mid,x,y);
    else if(x>mid) return get_rx(son[rt][],mid+,r,x,y);
    else return max(get_rx(son[rt][],mid+,r,mid+,y),get_all(son[rt][],mid+,r,mid+,y)+get_rx(son[rt][],l,mid,x,mid));
}
bool check(int k,int a,int b,int c,int d){
    int sum=;
    if(c>b+) sum+=get_all(rt[k],,n-,b+,c-);
    sum+=get_rx(rt[k],,n-,a,b);
    sum+=get_lx(rt[k],,n-,c,d);
    return sum>=;
}
int main(){
    scanf("%d",&n);
    for(int i=;i<n;i++){
        scanf("%d",&a[i]);id[i]=i;
    }
    sort(id,id+n,cmp);
    build(,n-,rt[]);
    for(int i=;i<n;i++)
        insert(rt[i-],rt[i],,n-,id[i-],-);
    int ord[];
    scanf("%d",&m);
    while(m--){
        scanf("%d%d%d%d",&ord[],&ord[],&ord[],&ord[]);
        for(int i=;i<;i++)
            ord[i]=(ord[i]+ans)%n;
        sort(ord,ord+);
        int l=,r=n,mid;
        while(l+<r){
            mid=l+r>>;
            if(check(mid,ord[],ord[],ord[],ord[]))l=mid;
            else r=mid;
        }
        ans=a[id[l]];
        printf("%d\n",ans);
    }
    return ;
}