（Sicily OJ）可持续化线段树

限制条件

时间限制: 1000 ms, 空间限制: 256 MB

题目描述

线段树有个非常经典的应用是处理RMQ问题，即区间最大/最小值询问问题。现在我们把这个问题可持久化一下:

Q k l r 查询数列在第k个版本时，区间[l, r]上的最大值

M k p v 把数列在第k个版本时的第p个数修改为v，并产生一个新的数列版本

最开始会给你一个数列，作为第1个版本。每次M操作会导致产生一个新的版本。

输入格式

第一行两个整数N, Q。N是数列的长度，Q表示询问数

第二行N个整数，是这个数列

之后Q行，每行以0或者1开头，0表示查询操作Q，1表示修改操作M，格式为

0 k l r 查询数列在第k个版本时，区间[l, r]上的最大值 或者

1 k p v 把数列在第k个版本时的第p个数修改为v，并产生一个新的数列版本

输出格式

对于每个M询问，输出正确答案

样例

input  

4 5
1 2 3 4
0 1 1 4
1 1 3 5
0 2 1 3
0 2 4 4
0 1 2 4

output  

4
5
4
4

 大致题意

两个操作

1在第k个版本的基础上修改第p个数形成新的版本

2查询第k个版本的（l,r）区间的最大值 

核心思路

构建线段树时，因为每次修改只需要一条边。所以只需要新增一条边就够了。而这一条边就是，你修改某一个点是遍历的路径。

于是，我们就可以在遍历时对该路径进行新增一条边，这一条边左（右）子节点连接新建的节点，另一个不变。

所以，我们需要一个数组ver记录第k代的线段树的根节点。

因为，每一次修改时只有遍历的路径变了，其他的没变，以此形成一个新的线段树。所以，从不同的根节点出发，一定会得到不同的线段树，只是这些线段树有相同的地方作为公共点（即节省空间和时间），但是在以这个根节点出发的关键路径，一定是区别于其他线段树的。

总结：可持续化线段树即通过每一代线段树之间相同的点，直接进行连边。并且新增一条关键路径（多个点），构成一颗新的线段树。

AC代码

#include <bits/stdc++.h>

using namespace std;

const int N = 10000 + 10;
const int Q = 100000 + 10;
const int LOG = 15;

struct Node {
	int val, lc, rc;
};

int a[N];
int ver[Q];
Node mem[2*N+Q*LOG];
int tot;

int build(int begin, int end) {
	if (begin == end) {
		mem[tot].val = a[begin];
		return tot ++;
	}
	int mid = (begin + end) / 2, u = tot ++;
	mem[u].lc = build(begin, mid);
	mem[u].rc = build(mid + 1, end);
	mem[u].val = max(mem[mem[u].lc].val, mem[mem[u].rc].val);
	return u;
}

int update(int root, int begin, int end, int loc, int val) {
	if (begin == end) {
		mem[tot].val = val;
		return tot ++;
	}
	int mid = (begin + end) / 2, u = tot ++;
	if (loc <= mid) {
		mem[u].lc = update(mem[root].lc, begin, mid, loc, val);
		mem[u].rc = mem[root].rc;
	}
	else {
		mem[u].lc = mem[root].lc; 
		mem[u].rc = update(mem[root].rc, mid + 1, end, loc, val);
	}
	mem[u].val = max(mem[mem[u].lc].val, mem[mem[u].rc].val);
	return u;
}

int query(int root, int begin, int end, int left, int right) {
	if (left <= begin && right >= end) return mem[root].val;
	int mid = (begin + end) / 2;
	if (right <= mid) return query(mem[root].lc, begin, mid, left, right);
	else if (left > mid) return query(mem[root].rc, mid + 1, end, left, right);
	else return max(query(mem[root].lc, begin, mid, left, mid), query(mem[root].rc, mid + 1, end, mid + 1, right));
}

int main() {
	int n, q, s = 0, op, k, p, v, l, r;
	scanf("%d %d", &n, &q);
	for (int i=1; i<=n; i++) scanf("%d", &a[i]);
	ver[++s] = build(1, n);
	while (q --) {
		scanf("%d %d", &op, &k);
		if (op == 0) {
			scanf("%d %d", &l, &r);
			printf("%d\n", query(ver[k], 1, n, l, r));
		}
		else {
			scanf("%d %d", &p, &v);
			ver[++s] = update(ver[k], 1, n, p, v);
		}
	}
	return 0;
}