摘抄了一段主席樹的解釋：所謂主席樹呢，就是對原來的數列[1..n]的每一個字首[1..i]（1≤i≤n）建立一棵線段樹，線段樹的每一個節點存某個字首[1..i]中屬於區間[L..R]的數一共有多少個（比如根節點是[1..n]，一共i個數，sum[root] = i；根節點的左兒子是[1..(L+R)/2]，若不大於(L+R)/2的數有x個，那麼sum[root.left] = x）。若要查詢[i..j]中第k大數時，設某結點x，那麼x.sum[j] - x.sum[i - 1]就是[i..j]中在結點x內的數字總數。而對每一個字首都建一棵樹，會MLE，觀察到每個[1..i]和[1..i-1]只有一條路是不一樣的，那麼其他的結點只要用回前一棵樹的結點即可，時空複雜度為O(nlogn)。

自己的程式碼醜陋，見諒

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn = 1e5 + 5;
struct node {
    int l, r, cnt;
    /*
    注意！！！
    l,r不是指區間，而是節點的序號
    cnt是有多少個數字出現
    */
} T[maxn * 20];
int len, n, m, a[maxn], t[maxn], root[maxn], tot;
// root 儲存的是所有歷史記錄
 

int build(int l, int r) {
    int rt = ++tot;
    T[rt].l = 0;
    T[rt].r = 0;
    T[rt].cnt = 0;//建立一棵空樹
    if(l == r) return rt;
    int mid = (l + r) >> 1;
    T[rt].l = build(l, mid);//左右遞迴建立
    T[rt].r = build(mid + 1, r);
    return rt;
}
int updata(int l, int r, int pre, int x) {
    int rt = ++tot;
    T[rt] = T[pre];//相當於連線新節點和以前節點的左右子樹
 

    T[rt].cnt++;
    if(l == r) return rt;
    int mid = (l + r) >> 1;//每次二分新增一個新節點
    if(mid >= x)T[rt].l = updata(l, mid, T[pre].l, x);
    else T[rt].r = updata(mid + 1, r, T[pre].r, x);
    return rt;
}
int query(int l, int r, int pre, int rt, int k) {
    if(l == r) return l;
    int mid = (l + r) >> 1;
    int sum = T[T[rt].l].cnt - T[T[pre].l].cnt;//每次比較兩個節點的左子樹
    if(sum >= k) return query(l, mid, T[pre].l, T[rt].l, k);
    else return query(mid + 1, r, T[pre].r, T[rt].r, k - sum);//如果在右子樹上別忘了減去sum
}
int main() {
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++) {
        scanf("%d", &a[i]);
        t[i] = a[i];
    }
    sort(t + 1, t + 1 + n);//必須離散化
    len = unique(t + 1, t + 1 + n) - t - 1;
    for(int i = 1; i <= n; i++) {
        a[i] = lower_bound(t + 1, t + 1 + len, a[i]) - t;
    }
    root[0] = build(1, n);
    for(int i = 1; i <= n; i++) {
        root[i] = updata(1, n, root[i - 1], a[i]);
    }
    for(int i = 1; i <= m; i++) {
        int u, v, w;
        scanf("%d%d%d", &u, &v, &w);
        printf("%d\n", t[query(1, n, root[u - 1], root[v], w)]);
    }
    return 0;
}