[河南省队2012] 找第k小的数

★★☆   输入文件:kth.in   输出文件:kth.out   简单对比
时间限制:1 s   内存限制:128 MB

题目描述

看到很短的题目会让人心情愉悦,所以给出一个长度为N的序列A1,A2,A3,...,AN,

现在有M个询问,每个询问都是Ai...Aj中第k小的数等于多少。

输入格式
第一行两个正整数N,M。
第二行N个数,表示序列A1,A2,...,AN。
紧着的M行,每行三个正整数i,j,k(k≤j-i+1),表示

询问Ai...Aj中第k小的数等于多少。

输出格式

共输出M行,第i行输出第i个询问的答案。

样例输入1:
4 3
4 1 2 3
1 3 1
2 4 3

1 4 4

样例输出1:
1
3
4
样例输入2:
5 5
4 2 9 9 10
1 3 1
2 4 3
1 4 4
3 5 2
2 5 2
样例输出2:
2
9
9
9
9
注释:
询问区间的第k小值并非严格第k小,例如样例2中第4个询问,询问3到5中第2小的数,
答案输出9,并不是严格第2小的10。
数据范围:
在50%的数据中,1<=N<=10000,1<=M<=10000,A[i]<=100000;
在100%的数据中,1<=N<=100000,1<=M<=100000,A[i]<=1000000;
思路:主席树+权值线段树
代码实现:
 1 #include<cstdio>
 2 #include<algorithm>
 3 using namespace std;
 4 const int maxn=1e5+10;
 5 int n,k,sz,ts;
 6 int a,b,c;
 7 int s[maxn],hs[maxn],tt[maxn];
 8 struct tree{int s,l,r,mid,lp,rp;}t[maxn<<5];
 9 void build(int l,int r,int k){
10     t[k].l=l,t[k].r=r;
11     t[k].mid=l+r>>1;
12     if(l==r) return;
13     t[k].lp=++ts,t[k].rp=++ts;
14     build(l,t[k].mid,t[k].lp);
15     build(t[k].mid+1,r,t[k].rp);
16 }
17 void put(int l,int r,int k,int nk,int p){
18     t[nk]=(tree){t[k].s+1,t[k].l,t[k].r,t[k].mid};
19     if(l==r) return;
20     if(p<=t[nk].mid){
21         t[nk].lp=++ts,t[nk].rp=t[k].rp;
22         put(l,t[nk].mid,t[k].lp,t[nk].lp,p);
23     }
24     else{
25         t[nk].rp=++ts,t[nk].lp=t[k].lp;
26         put(t[nk].mid+1,r,t[k].rp,t[nk].rp,p);
27     }
28 }
29 int search(int k,int nk,int v){
30     if(t[k].l==t[k].r) return t[k].l;
31     int w=t[t[nk].lp].s-t[t[k].lp].s;
32     if(v<=w) return search(t[k].lp,t[nk].lp,v);
33     else return search(t[k].rp,t[nk].rp,v-w);
34 }
35 int main(){
36     freopen("kth.in","r",stdin);
37     freopen("kth.out","w",stdout);
38     scanf("%d%d",&n,&k);
39     for(int i=1;i<=n;i++){
40         scanf("%d",&s[i]);
41         hs[i]=s[i];
42     }
43     sort(hs+1,hs+n+1);
44     sz=unique(hs+1,hs+n+1)-hs-1;
45     build(1,sz,tt[0]);
46     for(int i=1;i<=n;i++){
47         int pos=lower_bound(hs+1,hs+sz+1,s[i])-hs;
48         tt[i]=++ts;
49         put(1,sz,tt[i-1],tt[i],pos);
50     }
51     for(int i=1;i<=k;i++){
52         scanf("%d%d%d",&a,&b,&c);
53         printf("%d
",hs[search(tt[a-1],tt[b],c)]);
54     }
55     return 0;
56 }

钦定的板子。

题目来源:COGS

原文地址:https://www.cnblogs.com/J-william/p/6639845.html