划分树 求区间第k小问题

划分树是用来求区间第k小问题，区别于主席树，划分树的的内存消耗更小，主席树能够接树状数组或者线段树求动态区间第k小。

划分树的原理基于线段树，在给定的大n区间内部求任意区间第k大，我们将n个数不断按照大小顺序并且与原序列对比，再把比较后的值按下表保存在左右子树里，比中间值小的在左子树，比中间值大的在右子树，同时记录每个层每个点i和i点之前进入左子树点的个数（重点，这方便后面的查询）。

首先我们用网上流传最广的一张图（其实我懒得做新图）

原序列15638442，按照大小先排序，这里是另存一个数组，很多人都不明白这里建树的过程，问如果有sort的话，不是会超时吗，因为他只是排了一次序，然后每次都按照这个序列和原序列进行比较。

例如我排完序12344568，找到中间值4，那么我第一个进入树的值1，比4小，那么1就是左子树第一个，他和他之前进入左子树点的个数就是1（他自己），假设我们用一个数组num[][]来存，num[1][1]就是1，然后原数组下表2的位置是5，比4大那么他就进入右子树第一个，那么他之前还是只有一个1在左子树里，num[1][2]还是1，然后第三个位置6比4大，右子树第二个，num[1][3]=1、、以此类推，线段树的标准递归建树。
ps:中间值是需要特判的，这里有2个4，正好一个在左，一个在又，如果我们多加几个4，那么我就们就要判断有多少个4在左孩子里面，这里不能像线段树那样根据下标规则无脑放，这里的下标是不满足规则的。

那么建树看完了，说查找，这里有点小复杂，我们假设要求原数组第二位到第六位第3大的数(5)，那么我们就要从根（原序列）开始找，首先找到第6位比第2位进入左子树的个数多多少，假设多了c个，因为左子树值肯定是小于右子树的，所以满足第k小关系，假如这个k要小于c,那么我这个k小值就在左子树里，然后我们要缩小区间，这里也给一张图

我们找到k是在左区间，那么怎么缩小区间呢？首先我们要确定L--l里面在左子树的个数，这部分是我们要减去的，比如我现在求的l是这l-r区间最小的，但是不代表他从L-l是最小的，我们要相对的减去L-l在当前l之前的左子树个数，所以新的大区间起点应该是
L-num[dep][l-1]，相对的，如果k>c,那么第k小就是在右子树，我们的就要减去L-l中左子树里面的值，因为我们在求的在右子树里，这个k要变成
k-c,因为已经有c个在左子树了，所以只要求右子树第k-c小即可，同样的道理在重置缩小区间，直到查找到最后一层只有一个值的时候返回，这个值就是要求的第k小的值

完结上代码：模板题

---

#include<iostream>
#include<cstdio>
#include<algorithm> 
using namespace std;
const int maxn = 100010;
int tree[20][maxn],sorted[maxn],nleft[20][maxn];
void build(int l,int r,int dep){
	if(l==r) return;
	int mid = (l+r)>>1;
	int same = mid-l+1;  //same是有多少个和smid一样大小的值进入左孩子 
	int smid = sorted[mid];
	for(int i=l;i<=r;i++)
		if(tree[dep][i]<smid) same--;
	int lpos = l,rpos=mid+1;
	for(int i = l;i<=r;i++){
		nleft[dep][i]=(i==l)?0:nleft[dep][i-1];
		if( tree[dep][i] < smid ||tree[dep][i] == smid && same > 0){ //这里有个优先级关系，不能写成 tree[dep][i] <= smid,必须是满足same大于0，才可以进行前面判断 
			tree[dep+1][lpos++] = tree[dep][i];
			nleft[dep][i]++; //当前点进入左孩子的个数是这区间之前进入的总数+这个区间内进入左孩子的个数 
			if(tree[dep][i]==smid) same--;
		}else{
			tree[dep+1][rpos++] = tree[dep][i];
		} 
	}
	build(l,mid,dep+1); 
	build(mid+1,r,dep+1);
}
// L R是大区间 l,r是查询区间,查询第k大值
int search(int L,int R,int l,int r,int dep,int k){
	if(L==R) return tree[dep][R];
	int mid = (L+R)>>1;
	int ly = (l==L)?0:nleft[dep][l-1];//ly是l之前有多少个在左子树 
	int cnt = nleft[dep][r]-ly; //cnt 是这一层这个区间有多少个在左子树 
	if(k <= cnt){
		int newl = L+ly; //新的左起点 
		int newr = L+nleft[dep][r]-1;
		return search(L,mid,newl,newr,dep+1,k);
	}else{
		int newl=mid+1+(l-L-ly); //先去除掉L-l中左子树的个数 
		//	int newr=newl+(r-l+1)-cnt-1; //总个数减去这个区间左子树的个数，再去右边一个 
		int newr=newl+r-l-cnt;
		return search(mid+1,R,newl,newr,dep+1,k-cnt);
	}
}
int main(){
	int n,m,l,r,k;
	while(scanf("%d%d",&n,&m)!=EOF)
   {
     for(int i=1;i<=n;i++)
     {
       scanf("%d",&tree[0][i]);
       sorted[i]=tree[0][i];
     }
     sort(sorted+1,sorted+n+1);
     build(1,n,0);
     while(m--)
     {
       scanf("%d%d%d",&l,&r,&k);
       printf("%d\n",search(1,n,l,r,0,k));
     }
  }
	return 0;
}