主席树解决区间第k大

【填坑】可持久化线段树解决无修改的区间k大问题 - FAreStorm的博客 - 博客频道 - CSDN.NET
http://blog.csdn.net/farestorm/article/details/50337527
看上面这个学会的

我的理解是，一共有n个数，理想的方法是依次创建n颗权值线段树，第i颗线段树代表1到i的各个区间内存在的数的个数，那么，求解区间[s,t],便可通过第t颗减去第s-1颗得到，类似于前缀和，第k大的数二分区间得到，（别忘了线段树的结点定义），然而这样建树太浪费空间，应为有些区间在加入一个数后，并没有改变，所以，动态开点，保存下各个节点左右子树所对应的编号，最后为了进一步压缩空间，离散化原来的值。

例题

K-th Number POJ - 2104 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int maxn=5000000;

int ls[maxn]; //左儿子的编号
int rs[maxn]; //右儿子的编号
int num[maxn]; //有两次变换，第一次为原序列的值，第二次在lower_bound之后，为离散化后的值，也就是下标
int san[maxn]; //排序后原序列的值
int sum[maxn]; //线段树里保存的值
int T[maxn]; //第i颗线段树的编号
int n,m,tot = 0;

void build(int l,int r,int &x){
x=++tot;
sum[x]=0;
if(l==r) return;
int mid=(l+r)/2;
build(l,mid,ls[x]);
build(mid+1,r,rs[x]);
}

void update(int last,int p,int l,int r,int &x){
x=++tot;
ls[x]=ls[last];
rs[x]=rs[last];
sum[x]=sum[last]+1;
if(l==r) return;
int mid=(l+r)/2;
if(p<=mid) update(ls[last],p,l,mid,ls[x]);
else update(rs[last],p,mid+1,r,rs[x]);
}

int query(int s,int t,int l,int r,int k){
if(l==r) return l;
int mid=(l+r)/2;
int cnt=sum[ls[t]]-sum[ls[s]];
if(k<=cnt) return query(ls[s],ls[t],l,mid,k);
else return query(rs[s],rs[t],mid+1,r,k-cnt);
}

int main()
{
/* ios_base::sync_with_stdio(0);
cin.tie(0);*/
int n,m;
while(scanf("%d %d",&n,&m)!=EOF){
tot=0;
for(int i=0;i<n;i++){
scanf("%d",&san[i]);
num[i]=san[i];
}
sort(san,san+n);
int cnt=unique(san,san+n)-san;
build(1,cnt,T[0]);
for(int i=0;i<n;i++){
num[i]=lower_bound(san,san+cnt,num[i])-san+1;
update(T[i],num[i],1,cnt,T[i+1]);
}
while(m--){
int s,t,k;
scanf("%d %d %d",&s,&t,&k);
printf("%d\n",san[query(T[s-1],T[t],1,cnt,k)-1]);
}
}
return 0;
}