HDU2665 POJ2104 K-th Number(主席树)

给出一串长度为n的序列，并给出m个询问，每个询问针对序列中一段区间[s,t]，求这区间中第k小的数字是多少。

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 10;

struct node {
    int ls, rs, sum;
} ns[MAXN * 20];

int ct;
int rt[MAXN * 20];

void cpy(int& now, int old) {
    now = ++ct;
    ns[now] = ns[old];
}

void pushUp(int& now) {
    ns[now].sum = ns[ns[now].ls].sum + ns[ns[now].rs].sum;
}

//建树 
void build(int& now, int l, int r) {
    now = ++ct;
    ns[now].sum = 0;
    if (l == r) return;
    int m = (l + r) >> 1;
    build(ns[now].ls, l, m);
    build(ns[now].rs, m + 1, r);
}

//开始调用： update(rt[i], rt[i - 1], 1, sz, a[i])
//sz是去重后的数据个数 
//a[i]是原第i个数在去重后的序号 
void update(int& now, int old, int l, int r, int x) {
    cpy(now, old);
    if (l == r) {
        ns[now].sum++;
        return;
    }
    int m = (l + r) >> 1;
    if (x <= m) update(ns[now].ls, ns[old].ls, l, m, x);
    else update(ns[now].rs, ns[old].rs, m + 1, r, x);
    pushUp(now);
}

//查询第s棵到第t棵 第k大的树 
int query(int s, int t, int l, int r, int k) {
    if (l == r) return l;
    int m = (l + r) >> 1;
    int cnt = ns[ns[t].ls].sum - ns[ns[s].ls].sum;
    //cout << s << " " << t << " " << cnt << endl;
    if (k <= cnt) return query(ns[s].ls, ns[t].ls, l, m, k);
    return query(ns[s].rs, ns[t].rs, m + 1, r, k - cnt);
}

void init(int n) {
    ct = 0;
    build(rt[0], 1, n);
}

int a[MAXN], b[MAXN];

int main() {
    freopen("a.txt", "r", stdin);
    
    int T;
    scanf("%d", &T);
    while (T--) {
        int n, m;
        scanf("%d%d", &n, &m);
        for (int i = 1; i <= n; i++) {
            scanf("%d", &a[i]);
            b[i] = a[i];
        }
        sort(b + 1, b + n + 1);		//原数组排序 
        int sz = unique(b + 1, b + 1 + n) - b - 1;	//去重，得到不重复数据个数 
        init(sz);		//用序号建线段树 
        
        for (int i = 1; i <= n; i++) {
        	
            //求出原数在去重后的数据中的序号 
			a[i] = lower_bound(b + 1, b + 1 + sz, a[i]) - b;
            
            update(rt[i], rt[i - 1], 1, sz, a[i]);
        }
        
        printf("node count=%d\n",ct);
        for (int i = 0; i < ct; i++) {
            printf("%d, ls = %d, rs = %d, sum = %d\n", i,  ns[i].ls, ns[i].rs, ns[i].sum);
        }
        
        while (m--) {
            int s, t, k;
            scanf("%d%d%d", &s, &t, &k);
            printf("%d\n", b[query(rt[s - 1], rt[t], 1, sz, k)]);
        }
    }
    return 0;
}