可持久化线段树|主席树 POJ 2104 区间第k大的数



http://www.cnblogs.com/oyking/p/3230296.html

ps：整理了一个静态主席树的模板

struct Node {
    int L, R, sum;                              //左节点编号，右结点编号，和所有元素在区间内出现的总次数 
};
Node T[MAXN * 20];
int T_cnt;

struct A {
    int x, idx;
    bool operator < (const A &rhs) const {
        return x < rhs.x;
    }
};

A a[MAXN];
int rank[MAXN], root[MAXN];                    //rank[i]记录离散化后a[i]的值 

struct ZXT {
	void init() {
		T[0].L = T[0].R = T[0].sum = 0;
	   	root[0] = 0;
		for(int i = 1; i <= n; ++i) {
	        scanf("%d", &a[i].x);
    	    a[i].idx = i;
	    }
	    sort(a + 1, a + n + 1);
    	for(int i = 1; i <= n; ++i) rank[a[i].idx] = i;
	    T_cnt = 1;
    	for(int i = 1; i <= n; ++i) {
        	root[i] = root[i - 1];                   //root记录根的编号 
        	insert(rank[i], root[i], 1, n);
    	}
	}
	void insert(int &num, int &x, int L, int R) {   //num是插入的位置，即离散化后的值，x是节点的编号 
    	T[T_cnt++] = T[x]; x = T_cnt - 1;
    	++T[x].sum;
    	if(L == R) return ;
    	int mid = (L + R) >> 1;
    	if(num <= mid) insert(num, T[x].L, L, mid);
    	else insert(num, T[x].R, mid + 1, R);
	}

	int query(int i, int j, int k, int L, int R) {   //返回离散化后的值，即位置 
    	if(L == R) return L;
    	int t = T[T[j].L].sum - T[T[i].L].sum;
    	int mid = (R + L) >> 1;
    	if(k <= t) return query(T[i].L, T[j].L, k, L, mid);
    	else return query(T[i].R, T[j].R, k - t, mid + 1, R);
	}
	
	void print_ans() {
	//	printf("%d\n", a[query(root[i - 1], root[j], k, 1, n)].x);
	}
}