第K小数-可持续优化结构线段树-主席树

给定长度为N的整数序列A，下标为 1∼N1∼N。
现在要执行M次操作，其中第i次操作为给出三个整数li,ri,kili,ri,ki，求A[lili],A[lili+1],…,A[ri]ri中第kiki小的数是多少。
输入格式
第一行包含两个整数N和M。
第二行包含N个整数，表示整数序列A。
接下来M行，每行包含三个整数li,ri,kili,ri,ki，用以描述第i次操作。
输出格式
对于每次操作输出一个结果，表示在该次操作中，第k小的数的数值。
每个结果占一行。
数据范围
N≤105,M≤104,|A[i]|≤109N≤105,M≤104,|A[i]|≤109
输入样例：
7 3
1 5 2 6 3 7 4
2 5 3
4 4 1
1 7 3

输出样例：
5
6
3

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int N=100010,M=10010;
int n,m;
int a[N],idx,root[N];
vector<int> nums;
struct Node{
 int l,r;
 int cnt;
}tr[N*4+N*17];
int find(int x){
 return lower_bound(nums.begin(),nums.end(),x)-nums.begin();
}
int build(int l,int r){
 int p=++idx;
 if(l==r)    return p;
 int mid=l+r>>1;
 tr[p].l=build(l,mid),tr[p].r=build(mid+1,r);
 return p;
}
int insert(int p,int l,int r,int x){
 int q=++idx;
 tr[q]=tr[p];
 if(l==r){
  tr[q].cnt++;
  return q;
 }
 int mid=l+r>>1;
 if(x<=mid)   tr[q].l=insert(tr[p].l,l,mid,x);
 else        tr[q].r=insert(tr[p].r,mid+1,r,x);
 tr[q].cnt=tr[tr[q].l].cnt+tr[tr[q].r].cnt;
 return q;
}
int query(int q,int p,int l,int r,int k){
 if(l==r)   return r;
 int cnt=tr[tr[q].l].cnt-tr[tr[p].l].cnt;
 int mid=l+r>>1;
 if(k<=cnt)    return query(tr[q].l,tr[p].l,l,mid,k);
 else          return query(tr[q].r,tr[p].r,mid+1,r,k-cnt);
}
int main(){
 scanf("%d%d",&n,&m);
 for(int i=1;i<=n;i++){
  scanf("%d",&a[i]);
  nums.push_back(a[i]);
 }
 sort(nums.begin(),nums.end());
 nums.erase(unique(nums.begin(),nums.end()),nums.end());
 root[0]=build(0,nums.size()-1);
 for(int i=1;i<=n;i++)
    root[i]=insert(root[i-1],0,nums.size()-1,find(a[i]));
    while(m--){
     int l,r,k;
     scanf("%d%d%d",&l,&r,&k);
     printf("%d\n",nums[query(root[r],root[l-1],0,nums.size()-1,k)]);
    }
    return 0;
}