主席树（可持久化线段树）学习

前言：

可持久化线段树包括了主席树
主席树 其实就是可持久化权值线段树
个人建议 最好搞懂线段树和权值线段树再来学习主席树
个人看到网上的主席树只是解决区间第k大的问题，不知道还能有其他的应用吗，希望知道的能不吝指出

Reference：
权值线段树学习
https://blog.csdn.net/creatorx/article/details/75446472
https://blog.csdn.net/qq_39565901/article/details/81782739

---

首先先从权值线段树说起

可以看看这篇，我个人觉得写得很好易懂
权值线段树学习
简单来说就是，线段树维护的是一个数出现的个数，我们称之为权值，每出现一个数字，就在其所在的区间+1
权值线段树能解决的问题：
1.求查询一个数出现的次数
2.查询一段区间的数出现的次数
3.查询所有数的第k大值

---

那么什么是主席树

可以看看这个
主席树介绍
简单来说就是多棵权值线段树，第一棵树只维护了第一个节点，第二棵维护前2个节点…
第n棵维护前n个节点，很明显，我们的普通权值线段树就是主席树的第n棵节点
那么我们需要建n棵线段树吗？
很明显，我们第一棵树与第二棵树就改变了一个数，有很多重复的节点，所以我们只需建一棵线段树，其他的线段树由第一棵数拓展而来。
我们来看图：
（图片来自上面那篇博客）
加入有2，4，2 ，3这4个数，那么就有4颗线段树，那么如何变成主席树呢


上面是4棵线段树，下面就是一棵主席树了

主席树解决的问题：

很明显就是比权值线段树强一点，可以求区间的第k大
貌似就只能解决这种问题了，我还不知道能解决其他什么类型的问题，希望知道的可以不吝指出

---

先模拟一遍：

如果我们求[1,n]的第k大：
从根节点开始，判断右节点的权值是否大于等于k，如果是就进入右子树，说明第k大在右边。
否则进入左子树，但是k需要变，如果m=右子树的权值，k=k-m，进入左子树后就变为了在左子树中找第k-m大，这个很好理解。
求[l,r]区间的第k大：
那么我们就需要一个维护[l,r]区间的线段树，如何得到呢？就是用主席树的 第r棵线段树 - 第l棵线段树 得到 [l,r]的线段树
为什么呢？
因为假如[1,L]区间大于5的有3个，[1,R]区间大于5 的有10个，那么[L,R]区间大于5的就有10-3=7个数

---

看图模拟：

自己画的，推荐大家也收到画一遍，更加好理解
蓝色代表的是第3棵线段树维护的是[1,3]区间
红色代表的是第8棵线段树维护的是[1,8]区间
绿色代表的是维护[3,8]区间的线段树，树上的权值都用颜色表面了
（我这里为了方便，就把三棵树画在一起了）
不同的颜色代表不同数的权值，绿色的权值等于红色的权值-蓝色的权值
当得到绿色的线段树后，我们就直接使用普通的权值线段树求第k大，即可。
如图求[3,8]第5大，得到的是3

代码实现

以poj204为例，里面多了个离散化

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
using namespace std;
vector<int> v;
const int maxn=1e5+5;
int a[maxn];
int root[maxn];//多颗线段树的root节点
int m,n;
int cnt;//节点编号和个数
struct{int l,r,sum;}tree[maxn*20];
int getid(int x)
{
    return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}
void update(int l,int r,int &x,int y,int pos)
{
    tree[++cnt]=tree[y];
    tree[cnt].sum++;
    x=cnt;
    if(l==r) return;
    int mid=(l+r)/2;
    if(mid >= pos) update(l,mid,tree[x].l,tree[y].l,pos);
      else    update(mid+1,r,tree[x].r,tree[y].r,pos);

}
int query(int l,int r,int x,int y,int k)
{
    if(l==r) return l;
    int mid=(l+r)/2;
    int sum=tree[tree[y].l].sum  - tree[tree[x].l].sum;
    if(sum>=k) return query(l,mid,tree[x].l,tree[y].l,k);
    else return query(mid+1,r,tree[x].r,tree[y].r,k-sum);
}
int main()
{
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++)
    {
        scanf("%d", &a[i]);
        v.push_back(a[i]);
    }
    sort(v.begin(),v.end());
    v.erase( unique(v.begin(),v.end()) , v.end());//去重
    for(int i=1;i<=n;i++)
    {
        update(1,n,root[i],root[i-1],getid(a[i]));
    }
    int x,y,k;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&x,&y,&k);
        printf("%d\n",v[query(1,n,root[x-1],root[y],k) -1]);
    }
    return 0;
}