【主席树】区间查询第k大

洛谷P3834：

#include<bits/stdc++.h> 
using namespace std;
const int maxn = 1e5+5;
int a[maxn],d[maxn];
int T[maxn],tot;
struct SegTree{
    int l, r, sum;
}t[maxn*40];

int build(int l,int r){
    int p=++tot,mid=(l+r)>>1;
    if(l<r){
        t[p].l=build(l,mid);
        t[p].r=build(mid+1,r);
    }
    t[p].sum=0;return p;
}

int update(int pre,int l,int r,int pos){
    int p=++tot,mid=(l+r)>>1;
    t[p].l=t[pre].l;t[p].r=t[pre].r;t[p].sum=t[pre].sum+1;
    if(l<r){
        if(pos<=mid) t[p].l=update(t[pre].l,l,mid,pos);
        else t[p].r=update(t[pre].r,mid+1,r,pos); 
    }
    return p;
}

int query(int x,int y,int l,int r,int k){
    if(l==r) return l;
    int remain=t[t[y].l].sum - t[t[x].l].sum;int mid=(l+r)>>1;
    if(k<=remain) return query(t[x].l,t[y].l,l,mid,k);
    else return query(t[x].r,t[y].r,mid+1,r,k-remain);
}

int main(){
	int n,m;scanf("%d%d",&n,&m);tot=0;
    for (int i = 1; i <= n; i++)
        scanf("%d",&a[i]),d[i]=a[i];
    sort(d+1,d+1+n);
    int len=unique(d+1,d+1+n)-(d+1);
    for(int i=1;i<=n;i++) 
        a[i]=lower_bound(d+1,d+1+len,a[i])-d;

    T[0]=build(1, len);
    for(int i=1;i<=n;i++)
        T[i]=update(T[i-1],1,len,a[i]);
    
    while(m--){
        int l,r,k; scanf("%d%d%d",&l,&r,&k);
        int ans = query(T[l-1],T[r],1,len,k);
        printf("%d\n",d[ans]);
    }
    return 0;
}

#include<bits/stdc++.h>
#pragma GCC optimize(2)
#define ll long long
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define per(i,n,a) for(int i=n;i>=a;i--)
#define endl '\n'
#define eps 0.000000001
#define pb push_back
#define mem(a,b) memset(a,b,sizeof(a))
#define IO ios::sync_with_stdio(false);cin.tie(0);
using namespace std;
const int INF=0x3f3f3f3f;
const ll inf=0x3f3f3f3f3f3f3f3f;
const int mod=1e9+7;
const int maxn=2e5+5;
int a[maxn];
vector<int> v;
int getid(int x){
	return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}
struct node{
	int l,r,sum;
}hjt[maxn*40];
int cnt,root[maxn];
void _insert(int l,int r,int pre,int &now,int p){
	hjt[++cnt]=hjt[pre];
	now=cnt;
	hjt[now].sum++;
	if(l==r) return ;
	int mid=(l+r)>>1;
	if(p<=mid) _insert(l,mid,hjt[pre].l,hjt[now].l,p);
	else _insert(mid+1,r,hjt[pre].r,hjt[now].r,p);
}
int _query(int l,int r,int L,int R,int k){
	if(l==r) return l;
	int mid=(l+r)>>1;
	int tmp=hjt[hjt[R].l].sum-hjt[hjt[L].l].sum;
	if(k<=tmp) return _query(l,mid,hjt[L].l,hjt[R].l,k);
	else return _query(mid+1,r,hjt[L].r,hjt[R].r,k-tmp);
}
int main(){
	int n,m;scanf("%d%d",&n,&m);
	rep(i,1,n){
		scanf("%d",&a[i]);
		v.push_back(a[i]);
	}
	sort(v.begin(),v.end());
	v.erase(std::unique(v.begin(),v.end()),v.end());
	for(int i=1;i<=n;i++){
		_insert(1,n,root[i-1],root[i],getid(a[i]));
	}
	while(m--){
		int l,r,k;scanf("%d%d%d",&l,&r,&k);
		cout<<v[_query(1,n,root[l-1],root[r],k)-1]<<endl;
	}
}