【POJ 2104】 K-th Number 主席树模板题

达神主席树讲解传送门：http://blog.csdn.net/dad3zz/article/details/50638026

---

2016-02-23：真的是模板题诶，主席树模板水过。今天新校网不好，没有评测，但我立下flag这个代码一定能A。我的同学在自习课上考语文，然而机房党都跑到机房来避难了\(^o^)/~

2016-02-24：果然A乐~~~，我们机房党又躲过啦数学考试。良心不安啊

#include<cstdio>
#include<cstring>
#include<algorithm>
#define for1(i,a,n) for(int i=(a);i<=(n);i++)
#define for2(i,a,n) for(int i=(a);i<(n);i++)
#define for3(i,a,n) for(int i=(a);i>=(n);i--)
#define for4(i,a,n) for(int i=(a);i>(n);i--)
#define CC(i,a) memset(i,a,sizeof(i))
#define MID (l+r)>>1
#define read(x) x=getint()
using namespace std;
const int N=100010;
inline int max(const int &a,const int &b){return a>b?a:b;}
inline int min(const int &a,const int &b){return a<b?a:b;}
inline void swapp(int &a,int &b){int c=a;a=b;b=c;}
inline int getint(){char c;int ret=0,k=1;for(c=getchar();c<'0'||c>'9';c=getchar())if(c=='-')k=-1;for(;c>='0'&&c<='9';c=getchar())ret=ret*10+c-'0';return ret*k;}
int key,a[N],id[N],b[N],n,m,root[N],cnt=0;
struct node{int l,r,s;}T[N*50];
inline bool cmp(int x,int y){return a[x]<a[y];}
inline void update(int l,int r,int &pos){
	T[++cnt]=T[pos]; pos=cnt; T[pos].s++;
	if (l==r) return;
	int m=MID;
	if (key<=m) update(l,m,T[pos].l); else update(m+1,r,T[pos].r);
}
inline int query(int l,int r,int x,int y,int k){
	if (l==r) return l;
	int m=MID,s=T[T[y].l].s-T[T[x].l].s;
	if (k<=s) return query(l,m,T[x].l,T[y].l,k); else return query(m+1,r,T[x].r,T[y].r,k-s);
}
int main(){
	read(n); read(m);
	for1(i,1,n) read(a[i]),id[i]=i;
	sort(id+1,id+n+1,cmp);
	for1(i,1,n) b[id[i]]=i;
	for1(i,1,n){
		root[i]=root[i-1]; key=b[i];
		update(1,n,root[i]);
	}int x,y;
	while (m--){
		read(x); read(y); read(key);
		printf("%d\n",a[id[query(1,n,root[x-1],root[y],key)]]);
	}return 0;
}

这样就可以了，主席树水过了，然而这是DaD3zZ几年前就虐的东西~~(╯﹏╰)b，本蒟蒻还得努力啊。