主席树模板poj2104

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<utility>
#include<iostream>
#define debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
const int N=1000009;
struct aa
{
	int ls,rs,l,r,cnt;
}t[N*20];
int tot,a[N],b[N],nn,n,rt[N],q;
void build(int &now,int l,int r)
{
	now=++tot;
	t[now].l=l;t[now].r=r;
	if (l==r) return;
	int mid=(l+r)>>1;
	build(t[now].ls,l,mid);
	build(t[now].rs,mid+1,r);
}
int binary(int k)
{
	int l=1,r=nn,mid;
	while (l<=r)
	{
		mid=(l+r)>>1;
		if (b[mid]==k) return mid;
		if (b[mid]>k) r=mid-1;else l=mid+1;
	}
}
void insert(int last,int &now,int k)
{
	now=++tot;
	t[now]=t[last];
	t[now].cnt++;
	if (t[last].l==t[last].r) return;
	int mid=(t[now].l+t[now].r)>>1;
	if (mid>=k) insert(t[last].ls,t[now].ls,k);
	else insert(t[last].rs,t[now].rs,k);
}
int query(int la,int no,int k)
{
	if (t[la].l==t[la].r) return t[la].l;
	int h=t[t[no].ls].cnt-t[t[la].ls].cnt;
	if (h>=k) return query(t[la].ls,t[no].ls,k);
	return query(t[la].rs,t[no].rs,k-h);
}
int main()
{
	scanf("%d%d",&n,&q);
	for (int i=1;i<=n;i++) scanf("%d",&a[i]),b[i]=a[i];
	sort(b+1,b+n+1);
	//nn=n;
	nn=unique(b+1,b+n+1)-b-1;
	build(rt[0],1,nn);
	int x,y,z;
	for (int i=1;i<=n;i++) 
	{
		int k=binary(a[i]);
		insert(rt[i-1],rt[i],k);
	}
	while (q--)
	{
		scanf("%d%d%d",&x,&y,&z);
		printf("%d\n",b[query(rt[x-1],rt[y],z)]);
	}
	return 0;
}