主席树模板（区间第k大！）hdu2665

http://www.cnblogs.com/zyf0163/p/4749042.html

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<algorithm>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<string>
using namespace std;
typedef long long ll;
const double PI=acos(-1.0);
const double EXP=1e-8;
int Rank[100010],root[100010];
int cnt;
struct node1
{
int l,r,sum;
}node[100010*18];

struct point
{
int id,x;
}a[100010];

bool cmp(point p1,point p2)
{
return p1.x<p2.x;
}

void build(int &num,int &x,int l,int r)
{
node[cnt++]=node[x];
x=cnt-1;
node[x].sum++;
if(l==r)
return;
int mid=(l+r)/2;
if(num<=mid)
build(num,node[x].l,l,mid);
else
build(num,node[x].r,mid+1,r);
}

int query(int p,int q,int c,int l,int r)
{
if(l==r)
return l;
int mid=(l+r)/2;
int t=node[node[q].l].sum-node[node[p].l].sum;
if(c<=t)
return query(node[p].l,node[q].l,c,l,mid);
else
return query(node[p].r,node[q].r,c-t,mid+1,r);
}

int n,m;
int main()
{
node[0].l=node[0].r=node[0].sum=0;
root[0]=0;
int t;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i].x);
a[i].id=i;
}
sort(a+1,a+n+1,cmp);
for(int i=1;i<=n;i++)
Rank[a[i].id]=i;
cnt=1;
for(int i=1;i<=n;i++)
{
root[i]=root[i-1];
build(Rank[i],root[i],1,n);
}
while(m--)
{
int l,r,k;
scanf("%d%d%d",&l,&r,&k);
printf("%d\n",a[query(root[l-1],root[r],k,1,n)].x);
}
}
return 0;
}