可持久化线段树

可持久化线段树（又称函数式线段树）是一种可持久化数据结构（）。这种数据结构在普通线段树的基础之上支持查询某个历史版本，同时时间复杂度与线段树是同级，空间复杂度相较而言更高。
原理.
与大部分可持久化数据结构相似，可持久化线段树在更新时尽可能与之前某一个旧版本共用一部分结点，从而节省空间。
例如，有一颗维护着区间和的线段树，现在以这颗线段树作为基础，将下标为 formula_1 的元素数值减去一，同时另存为一个新的版本。按照一般线段树中使用的思路，当位于根结点时我们应该递归操作它的右子树，如果是暴力实现的可持久化线段树，那么则需要再次递归将一整颗左子树复制一遍然后保存指针，但是复制的这一整颗子树是完全一模一样的，因此可以先只复制一个根结点，然后将它的左子树直接指向原先版本根结点的左子树，代表上一个版本和这一个新版本这颗子树保存的信息是完全一样的，然后再按照相似的方法，递归地处理下标 formula_1 存在的右子树。
性能分析.
静态可持久化线段树在每次更新版本时，总是最大程度上减少结点的复制，这样不仅减少了时间的开销，更避免了不必要的空间浪费。与线段树相同，可持久化线段树由于在每一次更新版本时没有访问到不必要的结点，所以每一次查询和修改（即建立一个新的版本）时间复杂度为 formula_3, 在这个过程中，会同时建立 formula_3 个新的结点。如果总操作数量为 formula_5, 那么可持久化线段树可以以 formula_6 的时空复杂度解决问题。
实践.
静态区间第k大数值.
这类问题需要求解在一个长度为 formula_7 的数列中，第 formula_8 个数为 formula_9. 现在给定一些形如 formula_10 的三元组作为询问，设计程序计算数列第 formula_11 这些元素中出现次数排在第 formula_12 位的是多少。
利用可持久化线段树，维护区间 formula_13 代表区间 formula_14 中的元素出现了多少次，以此作为原始版本 formula_15，此后每次建立一个新版本 formula_16，代表去掉原数列中 formula_17 的元素之后建立的线段树，维护目标与上述相同。具体过程可以每次将 formula_9 的出现次数减一，并保存此时生成的新版本。
参考程序.
C++.
constexpr auto MAXN = (int)2e5 + 500;
std::map&lt;int, int&gt; val, mp;
struct Node {
 int fr, to, sum;
 Node *lft, *rgt;
 Node&amp; Copy(Node *targ) {
 fr = targ-&gt;fr; to = targ-&gt;to; sum = targ-&gt;sum;
 lft = targ-&gt;lft; rgt = targ-&gt;rgt;
 return *this;
};
Node* NewNode() {
 Node* ret = (Node*)malloc(sizeof(Node));
 ret-&gt;lft = ret-&gt;rgt = nullptr;
 return ret;
std::vector&lt;Node*&gt; version;
int num[MAXN], cnt[MAXN], orig[MAXN], fr[MAXN], to[MAXN], rank[MAXN];
std::queue&lt;Node*&gt; que, add;
signed main(void)
 int totNums, totQuery, cnt;
 //Read
 scanf("%d%d", &amp;totNums, &amp;totQuery);
 for (int i = 0; i &lt; totNums; i++)scanf("%d", num + i);
 for (int i = 0; i &lt; totQuery; i++)scanf("%d%d%d", fr + i, to + i, rank + i);
 for (int i = 0; i &lt; totNums; i++) orig[i] = num[i];
 std::sort(num, num + totNums);
 cnt = 0; int tmp;
 mp[val[0] = *num] = cnt++;
 for (int i = 1; i &lt; totNums; i++)
 if (num[i] != num[i - 1])
 tmp = cnt,mp[val[tmp] = num[i]] = cnt++; 
 for (int i = 0; i &lt; totNums; i++)::cnt[num[i] = mp[orig[i]]]++;
 //Build
 Node *a, *b, *t;
 for (int i = 0; i &lt; cnt; i++) {
 t = NewNode(); t-&gt;fr = t-&gt;to = i;
 t-&gt;sum = ::cnt[i];
 que.push(t);
 for (; que.size() &gt;= 2; std::swap(que, add)) {
 while (que.size() &gt;= 2) {
 a = que.front(); que.pop(); b = que.front(); que.pop();
 t = NewNode(); t-&gt;fr = a-&gt;fr; t-&gt;to = b-&gt;to; t-&gt;sum = a-&gt;sum + b-&gt;sum;
 t-&gt;lft = a; t-&gt;rgt = b;
 add.push(t);
 }version.push_back(que.front());
 //New Versions
 for (int del = 0; del &lt; totNums; del++) {
 const int &amp;target = num[del];
 t = a = NewNode(); b = version.back();
 while (true) {
 a-&gt;Copy(b); a-&gt;sum--;
 if (a-&gt;fr == target &amp;&amp; a-&gt;to == target)break;
 if (a-&gt;lft-&gt;fr &lt;= target &amp;&amp; target &lt;= a-&gt;lft-&gt;to) {
 a-&gt;lft = NewNode(); a = a-&gt;lft; b = b-&gt;lft;
 } else {
 a-&gt;rgt = NewNode(); a = a-&gt;rgt; b = b-&gt;rgt;
 }version.push_back(t);
 //Query
 int rnk;
 for (int i = 0; i &lt; totQuery; i++)
 a = version[--fr[i]]; b = version[to[i]]; rnk = rank[i];
 while (true) {
 if (a-&gt;lft-&gt;sum - b-&gt;lft-&gt;sum &gt;= rnk) {
 a = a-&gt;lft; b = b-&gt;lft;
 else {
 rnk -= a-&gt;lft-&gt;sum - b-&gt;lft-&gt;sum;
 a = a-&gt;rgt; b = b-&gt;rgt;
 
 //system("pause");
 return 0;