“玲珑杯”ACM比赛 Round #20 E -- 造物主的戒律（主席树）

给你一个序列，每次查询区间中小于等于x的所有数字里面第k1小的值以及大于x的所有数字里面第k2小的值，如果不存在，输出-1
每次输出两个数，对于每个数如果不存在，则单独输出-1
INPUT
第一行两个数n,m
第二行n个数表示序列a
后面m行每行五个数l,r,x,k1,k2
OUTPUT
对于每个查询输出答案
SAMPLE INPUT
5 5
1 2 3 4 5
1 5 3 1 2
2 5 2 1 3
2 3 3 3 3
3 3 3 3 3
3 3 3 3 3
SAMPLE OUTPUT
1 5
2 5
-1 -1
-1 -1
-1 -1
题解：很明显需要用主席树来维护。先看“区间中小于等于x的所有数字里面第k1小的值”这个问，我们先正常求出区间第k1小（设为cnt），则cnt和x只有两种关系，1：cnt<=x,则cnt就是我们要求的。2：cnt>x,则这个区间没有我们所求的数，就是-1。再看“大于x的所有数字里面第k2小的值“这个问，我们先把区间内<=x的数的个数先求出来（设为tot），则这个问实际上就是求区间第k2+tot小。
代码：

#include<iostream>#include<stdio.h>#include<stdlib.h>#include<string.h>#include<vector>#include<queue>#include<set>#include<algorithm>#include<map>#include<math.h>using namespace std;typedef pair<int,int>pa;const int N=5e5+100;int read(){    int x=0;    char ch = getchar();    while('0'>ch||ch>'9')ch=getchar();    while('0'<=ch&&ch<='9')    {        x=(x<<3)+(x<<1)+ch-'0';        ch=getchar();    }    return x;}/***********************************************************/vector<int>p;int cnt,n,m,tmp,k,root[N],a[N];int ll[N],rr[N],x[N],k1[N],k2[N];struct node{    int l,r,sum;} T[N*40];int getid(int x) {return lower_bound(p.begin(),p.end(),x)-p.begin()+1;}int ans1,ans2,ans;void update(int l,int r,int &x,int y,int pos){    T[++cnt]=T[y],T[cnt].sum++,x=cnt;    if(l==r) return ;    int mid=(r+l)>>1;    if(mid>=pos)        update(l,mid,T[x].l,T[y].l,pos);    else        update(mid+1,r,T[x].r,T[y].r,pos);}int query(int l,int r,int x,int y,int k){    if(l==r) return l;    int mid=(r+l)>>1;    int sum=T[T[y].l].sum-T[T[x].l].sum;    if(sum>=k)        return query(l,mid,T[x].l,T[y].l,k);    else        return query(mid+1,r,T[x].r,T[y].r,k-sum);}int query_2(int l,int r,int x,int y,int L,int R){    if(L<=l&&R>=r) return T[y].sum-T[x].sum;    int mid=(r+l)>>1,ans=0;    if(L<=mid) ans+=query_2(l,mid,T[x].l,T[y].l,L,R);    if(R>mid) ans+=query_2(mid+1,r,T[x].r,T[y].r,L,R);    return ans;}int main(){    scanf("%d%d",&n,&m);    for(int i=1; i<=n; i++) scanf("%d",&a[i]),p.push_back(a[i]);    for(int i=1; i<=m; i++) scanf("%d%d%d%d%d",&ll[i],&rr[i],&x[i],&k1[i],&k2[i]),p.push_back(x[i]);    sort(p.begin(),p.end()),p.erase(unique(p.begin(),p.end()),p.end());    int sz=p.size();    for(int i=1; i<=n; i++) update(1,sz,root[i],root[i-1],getid(a[i]));    for(int i=1; i<=m; i++)    {        ans1=ans2=-1;        if(rr[i]-ll[i]+1>=k1[i]&&k1[i]>0)        {            ans=query(1,sz,root[ll[i]-1],root[rr[i]],k1[i]);            if(ans<=getid(x[i]))                ans1=p[ans-1];        }        ans=query_2(1,sz,root[ll[i]-1],root[rr[i]],1,getid(x[i]));        if(rr[i]-ll[i]+1>=ans+k2[i]&&k2[i]>0)        {            tmp = query(1,sz,root[ll[i]-1],root[rr[i]],ans+k2[i]);            ans2=p[tmp-1];        }        printf("%d %d\n",ans1,ans2);    }    return 0;}