【poj2104-求区间第k大数（不修改）】主席树/可持续化线段树

第一道主席树~然而是道比较水的。。。因为它不用修改。。。

转载一个让我看懂的主席树的讲解吧：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\n",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\n",num[query(root[l-1],root[r],k)]);
122         }
123     }
124     return 0;
125 }