【数据结构模版】可持久化线段树 && 主席树

　　浙江集训Day4，从早8:00懵B到晚21：00，只搞懂了可持久化线段树以及主席树的板子。今天只能记个大概，以后详细完善讲解。

　　可持久化线段树指的是一种基于线段树的可回溯历史状态的数据结构。我们想要保存某一段序列的历史信息，可以朴素地开m倍的空间存储；而线段树是一种优秀的数据结构，它每次修改操作只把原先的数据修改log个节点，这就意味着我们只增加这些被修改的节点存储新信息就好了。传统的线段树用二叉树存储，但是因为涉及加点，可持久化线段树必须维护一个自增的tot动态成为新点的编号。对于每个节点，我们要维护两个指针指向它的左右子节点。对于每次修改，我们把新节点的不被修改的一半区间指向上个状态该区间的节点编号，等同于利用了上次的区间信息；另一半新开一个点继续递归下去即可。

　　可持久化线段树一般不能维护区间修改，因为多个根共享子节点，我们没有办法对每个状态的树独立维护它的lazy_tag。

代码（动态开点&&区间最大值）：

1. #include <iostream>
2. #include <cstdio>
3. #include <cstring>
4. #include <algorithm>
5. #include <climits>
6. #define maxn 1000100
7. #define maxm 1000100
8. #define LG 21
9. #define BUG putchar('*')
10. using namespace std;
11. template <typename T>
12. void read(T &x) {
13. x = 0;
14. int f = 1;
15. char ch = getchar();
16. while (!isdigit(ch)) {
17. if (ch == '-')
18. f = -1;
19. ch = getchar();
20. }
21. while (isdigit(ch)) {
22. x = x * 10 + (ch ^ 48);
23. ch = getchar();
24. }
25. x *= f;
26. return;
27. }
28. int n, m, a[maxn];
29. namespace Segment_tree {
30. int root[maxm], tot;
31. struct node {
32. int dat, lc, rc;
33. } seg[maxn * 20];
34. inline void update(int nd) {
35. seg[nd].dat = max(seg[seg[nd].lc].dat, seg[seg[nd].rc].dat);
36. }
37. int build(int l, int r) {
38. int nd = ++tot;
39. if (l == r) {
40. seg[nd].dat = a[l];
41. return nd;
42. }
43. int mid = (l + r) >> 1;
44. seg[nd].lc = build(l, mid), seg[nd].rc = build(mid + 1, r);
45. update(nd);
46. return nd;
47. }
48. int modify(int nd, int l, int r, int x, int val) {
49. int now = ++tot;
50. seg[now] = seg[nd];
51. if (l == r) {
52. seg[now].dat = val;
53. return now;
54. }
55. int mid = (l + r) >> 1;
56. if (x <= mid)
57. seg[now].lc = modify(seg[nd].lc, l, mid, x, val);
58. else seg[now].rc = modify(seg[nd].rc, mid + 1, r, x, val);
59. update(now);
60. return now;
61. }
62. int query(int nd, int l, int r, int ql, int qr) {
63. if (ql <= l && qr >= r)
64. return seg[nd].dat;
65. if (ql > r || qr < l)
66. return INT_MIN;
67. int mid = (l + r) >> 1;
68. return max(query(seg[nd].lc, l, mid, ql, qr), query(seg[nd].rc, mid + 1, r, ql, qr));
69. }
70. } using namespace Segment_tree;
71. int main() {
72. read(n), read(m);
73. for (int i = 1; i <= n; ++i)
74. read(a[i]);
75. root[0] = build(1, n);
76. int k, op, val, loc, l, r;
77. for (int i = 1; i <= m; ++i) {
78. read(k), read(op);
79. if (op == 1) {
80. read(loc), read(val);
81. root[i] = modify(root[k], 1, n, loc, val);
82. } else {
83. read(l), read(r);
84. printf("%d\n", query(root[k], 1, n, l, r));
85. root[i] = root[k];
86. }
87. }
88. return 0;
89. }

　　听说主席树是由fotile96主席发明的（总之是位神犇就对了）。主要用于维护数据区间的值域问题。

　　给定一段区间，每次询问l,r区间内的第k小值。主席树是一个可持久化权值线段树：我们将右端点r抽象成时刻，建立（数据长度）个根，以root[t]表示区间[1, t]所维护的权值线段树的根。（模板题的权值需要离散化）如果查询[l, r]，我们找出第l-1和第r棵线段树（将root[0]建成一棵空树），这两棵树对应节点信息做差，实质上就得到了一棵仅维护区间[l, r]的虚拟线段树的信息。我们动态跑出这棵树的同时二分它的左右子树：如果这棵虚拟树的左端点权值和sum大于等于k，说明第k大数在[l, mid]中，我们转而查询这段的第k大：如果sum小于k，说明所求数是区间[mid + 1, r]的第k-sum大数。实际上，在区间Kth数问题中，主席树以类似前缀和的思想维护了虚拟的维护[l, r]信息的权值线段树来查询对应区间信息的。

　代码：

1. #include <iostream>
2. #include <cstdio>
3. #include <cstring>
4. #include <algorithm>
5. #include <climits>
6. #define maxn 200100
7. #define maxm 200100
8. #define LG 21
9. #define BUG putchar('*')
10. using namespace std;
11. int n, m;
12. template <typename T>
13. void read(T &x) {
14. x = 0;
15. int f = 1;
16. char ch = getchar();
17. while (!isdigit(ch)) {
18. if (ch == '-')
19. f = -1;
20. ch = getchar();
21. }
22. while (isdigit(ch)) {
23. x = x * 10 + (ch ^ 48);
24. ch = getchar();
25. }
26. x *= f;
27. return;
28. }
29. namespace Contra {
30. int st[maxn];
31. int work(int *org, int *st) {
32. for (int i = 1; i <= n; ++i)
33. st[i] = org[i];
34. sort(st + 1, st + 1 + n);
35. int len = unique(st + 1, st + n + 1) - (st + 1);
36. for (int i = 1; i <= n; ++i)
37. org[i] = lower_bound(st + 1, st + len + 1, org[i]) - st;
38. return len;
39. }
40. }
41. int a[maxn], N;
42. namespace President_tree {
43. #define lc(i) seg[(i)].lc
44. #define rc(i) seg[(i)].rc
45. using namespace Contra;
46. int tot, root[maxm];
47. struct node {
48. int cnt, lc, rc;
49. } seg[maxm * 20];
50. void update(int nd) {
51. seg[nd].cnt = seg[lc(nd)].cnt + seg[rc(nd)].cnt;
52. }
53. int build(int l, int r) {//root[0]:empty_tree
54. int nd = ++tot;
55. if (l == r) return nd;
56. int mid = (l + r) >> 1;
57. seg[nd].lc = build(l, mid);
58. seg[nd].rc = build(mid + 1, r);
59. return nd;
60. }
61. int modify(int pre, int l, int r, int x) {
62. int nd = ++tot;
63. seg[nd] = seg[pre];
64. if (l == r) {
65. ++seg[nd].cnt;
66. return nd;
67. }
68. int mid = (l + r) >> 1;
69. if (x <= mid)
70. lc(nd) = modify(lc(pre), l, mid, x);
71. else
72. rc(nd) = modify(rc(pre), mid + 1, r, x);
73. update(nd);
74. return nd;
75. }
76. int query(int ql, int qr, int l, int r, int k) {
77. if (l == r) {
78. return st[l];
79. }
80. int mid = (l + r) >> 1, sum = seg[lc(qr)].cnt - seg[lc(ql)].cnt;
81. if (k <= sum)
82. return query(lc(ql), lc(qr), l, mid, k);
83. return query(rc(ql), rc(qr), mid + 1, r, k - sum);
84. }
85. void init() {
86. N = work(a, st);
87. root[0] = build(1, N);
88. for (int i = 1; i <= n; ++i) {
89. root[i] = modify(root[i - 1], 1, N, rk[i]);
90. }
91. }
92. } using namespace President_tree;
93. int main() {
94. read(n), read(m);
95. for (int i = 1; i <= n; ++i)
96. read(a[i]);
97. init();
98. int l, r, k;
99. for (int i = 1; i <= m; ++i) {
100. read(l), read(r), read(k);
101. printf("%d\n", query(root[l - 1], root[r], 1, N, k));
102. }
103. return 0;
104. }