#include <cstdio>
#include <algorithm>
#define MAXN (100000)
#define DEEP (20)
using namespace std;
int sorted[DEEP][MAXN], a[MAXN];
void build(int deep, int lft, int rht){//传入深度,区间左右界,目标就是更新好SORTED数组
if(lft == rht){//区间只有一个数就是结束条件,SORTED[当前深度][位置编号]就是A数组值
sorted[deep][lft] = a[lft];
return ;
}
int mid = (lft + rht) >> 1;//得到区间中间位置
build(deep+1, lft, mid);//左半区间建树
build(deep+1, mid+1, rht);//右半区间建树
int p = lft, q = mid+1, k = lft;//读出左半区间左界,右半区间左界,整个区间的左界
while(p <= mid && q <= rht){//归并排序
if(sorted[deep+1][p] <= sorted[deep+1][q])//下一深度的P位数值比Q位小
sorted[deep][k++] = sorted[deep+1][p++];//本深度第K位的数就记为下一深度P位的数
else sorted[deep][k++] = sorted[deep+1][q++];//否则本深度第K位的数就记为下一深度Q位的数
}
while(p <= mid) sorted[deep][k++] = sorted[deep+1][p++];//如果Q先排完,就再把P以后的数打入SORTED
while(q <= rht) sorted[deep][k++] = sorted[deep+1][q++];//如果P先排完,就再把Q以后的数打入SORTED
}
int query(int deep, int lft, int rht, int qlft, int qrht, int key){//本层深度,搜索区间左右界,询问区间左右界,本层中间值
if(qrht < lft || qlft > rht) return 0;//只要区间不是完全包含,就表进行统计,结束递归
if(qlft <= lft && rht <= qrht)//搜索区间完全处于询问区间内时返回询问值在当前搜索区间中的排名
return lower_bound(&sorted[deep][lft],&sorted[deep][rht]+1,key) - &sorted[deep][lft];
int mid = (lft + rht) >> 1;//取出区间中点位置
return query(deep+1, lft, mid, qlft, qrht, key) + query(deep+1, mid+1, rht, qlft, qrht, key);//递归左右半区间共有多少个数小于KEY值
}
int solve(int n, int qlft, int qrht, int k){//总点数(即区间数),待询问左右界,第K大数
int low = 1, high = n+1;//二分初始化左右界
while(low+1 < high){//未结束
int mid= (low + high) >> 1;//二分答案
int cnt = query(0, 1, n, qlft, qrht, sorted[0][mid]);//sorted[0][]是对整个区间排好序的!,所以就是问整个区间的第MID大的值在询问区间中第几大
if(cnt <= k) low = mid;//得到的数的排名小于K,说明太小了,所以下界要变大
else high = mid;//否则上界变大
}
return sorted[0][low];
}
int main(){
int n, m;
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++)
scanf("%d", &a[i]);
build(0, 1, n);//建树
while(m--){
int qlft, qrht, k;
scanf("%d%d%d", &qlft, &qrht, &k);
printf("%d\n", solve(n, qlft, qrht, k-1));
}
return 0;
}
归并树
POJ2104第K小数
Sample Input
7 3
1 5 2 6 3 7 4
2 5 3
4 4 1
1 7 3
Sample Output
5
6
3