第一道主席树~然而是道比较水的。。。因为它不用修改。。。
转载一个让我看懂的主席树的讲解吧:http://blog.csdn.net/regina8023/article/details/41910615 (未授权,侵权删)
---------------------------------------------------------------------------------------------------------------------
那么如果要询问i-j之间数字出现的次数怎么办呢?
因为每一棵线段树的区间都是相同的,所以要求l-r之间的数字的出现次数只要用前r位出现的次数减去前l-1位出现的次数,就是ans
但是如果有修改操作怎么办?
如果沿用上面的做法,那么修改操作是O(nlogn)的,查询是O(1)的,修改要花好长时间。。。
前缀和联想到了树状数组,那么将前缀和用树状数组维护的话修改是O(logn*logn),查询时O(logn),查询的时间虽然变长,但是修改的时间缩短许多!!
注意:
函数式线段树的数组要开大一点!!
---------------------------------------------------------------------------------------------------------------------
这题就是模版题啦,先离散,求区间第k大的时候lx=root[l-1],rx=root[r],两边同时走不断作差,看看左孩子的数量,如果k更大就减掉左孩子的到有右孩子中找。
代码:
1 #include<cstdio>
2 #include<cstdlib>
3 #include<cstring>
4 #include<cmath>
5 #include<iostream>
6 #include<algorithm>
7 using namespace std;
8
9 const int N=100010,INF=(int)1e9+100;
10 struct trnode{
11 int lc,rc,cnt;
12 }t[30*N];
13 struct node{
14 int d,id;
15 }p[N];
16 int n,m,tl,mx;
17 int a[N],num[N],root[N];
18
19 bool cmp(node x,node y){return x.d<y.d;}
20
21 int bt(int l,int r)
22 {
23 int x=++tl;
24 // a[x].l=l;a[x].r=r;
25 t[x].lc=t[x].rc=0;
26 t[x].cnt=0;
27 if(l<r)
28 {
29 int mid=(l+r)/2;
30 t[x].lc=bt(l,mid);
31 t[x].rc=bt(mid+1,r);
32 }
33 return x;
34 }
35
36 int add(int rt,int x)
37 {
38 int now=++tl,tmp=now;
39 t[now].cnt=t[rt].cnt+1;
40 int l=1,r=mx,mid;
41 while(l<r)
42 {
43 mid=(l+r)/2;
44 if(x<=mid)
45 {
46 t[now].lc=++tl;
47 t[now].rc=t[rt].rc;
48 rt=t[rt].lc;
49 now=tl;
50 r=mid;
51 }
52 else
53 {
54 t[now].lc=t[rt].lc;
55 t[now].rc=++tl;
56 rt=t[rt].rc;
57 now=tl;
58 l=mid+1;
59 }
60 t[now].cnt=t[rt].cnt+1;
61 }
62 return tmp;
63 }
64
65 int query(int lx,int rx,int k)
66 {
67 int l=1,r=mx,mid;
68 while(l<r)
69 {
70 mid=(l+r)/2;
71 if(t[t[rx].lc].cnt-t[t[lx].lc].cnt>=k)
72 {
73 r=mid;
74 lx=t[lx].lc;
75 rx=t[rx].lc;
76 }
77 else
78 {
79 l=mid+1;
80 k-=t[t[rx].lc].cnt-t[t[lx].lc].cnt;
81 lx=t[lx].rc;
82 rx=t[rx].rc;
83 }
84 }
85 return l;
86 }
87
88 void output(int x)
89 {
90 printf("x = %d lc = %d rc = %d cnt = %d
",x,t[x].lc,t[x].rc,t[x].cnt);
91 if(t[x].lc) output(t[x].lc);
92 if(t[x].rc) output(t[x].rc);
93 }
94
95 int main()
96 {
97 freopen("a.in","r",stdin);
98 while(scanf("%d%d",&n,&m)!=EOF)
99 {
100 tl=0;mx=0;
101 for(int i=1;i<=n;i++)
102 {
103 scanf("%d",&a[i]);
104 p[i].d=a[i];p[i].id=i;
105 }
106 sort(p+1,p+1+n,cmp);
107 p[0].d=INF;
108 for(int i=1;i<=n;i++)
109 {
110 if(p[i].d!=p[i-1].d) mx++,num[mx]=p[i].d;
111 a[p[i].id]=mx;
112 }
113 root[0]=bt(1,mx);
114 for(int i=1;i<=n;i++)
115 root[i]=add(root[i-1],a[i]);
116 // output(root[5]);
117 for(int i=1;i<=m;i++)
118 {
119 int l,r,k;
120 scanf("%d%d%d",&l,&r,&k);
121 printf("%d
",num[query(root[l-1],root[r],k)]);
122 }
123 }
124 return 0;
125 }