uoj218火车管理 可持久化线段树+线段树

uoj218火车管理

题目连接

分析

如果没有弹栈操作，就是一个普通的线段数+区间覆盖。
有了弹栈操作之后，一个非常巧妙的思路是，用主席树维护当前每个节点的插入时间。
当弹出一个节点的时候，用主席树找到这个位置的插入时间。调出这个插入时间的前一个时刻版本的线段树，这个版本记录了当前位置的栈的下一个个数的插入时间。记录一下插入操作的吨位，这样就可以在对应的线段树和主席树上修改了。
有一个神奇的小操作来实现主席树的区间覆盖标记永久化。那就是遇到了一个被覆盖的节点，把这个节点的左右儿子指向这个节点本身。这个操作可以用于任意的线段树（当然要记录左右子树，所以在动态开点的线段树比较吃香），可以省时间省空间。

代码

#include<bits/stdc++.h>
const int N = 5e5 + 10, T = 4e7 + 10;
int ri() {
	char c = getchar(); int x = 0, f = 1; for(;c < '0' || c > '9'; c = getchar()) if(c == '-') f = -1;
	for(;c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) - '0' + c; return x * f;
}
int val[T], ls[T], rs[T], rt[N], u[N], sz, n, m, ty;
long long Lans;
namespace seg {
	const int Ts = N << 2;
	#define lc p << 1
	#define rc p << 1 | 1
	int tg[Ts]; long long s[Ts];
	void Tag(int p, int L, int R, int v) {tg[p] = v; s[p] = 1LL * (R - L + 1) * v;}
	void Push(int p, int L, int m, int R) {
		if(tg[p]) {
			Tag(lc, L, m, tg[p]);
			Tag(rc, m + 1, R, tg[p]);
			tg[p] = 0;
		}
	}
	void Mdf(int p, int L, int R, int st, int ed, int v) {
		if(L == st && ed == R)
			return Tag(p, L, R, v);
		int m = L + R >> 1; Push(p, L, m, R);
		if(st <= m) Mdf(lc, L, m, st, std::min(ed, m), v);
		if(ed > m) Mdf(rc, m + 1, R, std::max(st, m + 1), ed, v);
		s[p] = s[lc] + s[rc];
	}
	long long Qry(int p, int L, int R, int st, int ed) {
		if(L == st && ed == R)
			return s[p];
		int m = L + R >> 1; long long ans = 0; Push(p, L, m, R);
		if(st <= m) ans += Qry(lc, L, m, st, std::min(ed, m));
		if(ed > m) ans += Qry(rc, m + 1, R, std::max(st, m + 1), ed);
		return ans;
	}
}
int Qry(int rt, int L, int R, int v) {
	if(L == R) return val[rt]; int m = L + R >> 1;
	return v <= m ? Qry(ls[rt], L, m, v) : Qry(rs[rt], m + 1, R, v);
}
void Mdf(int &cp, int lp, int L, int R, int st, int ed, int x) {
	cp = ++sz;
	if(L == st && ed == R)
		return val[cp] = x, ls[cp] = rs[cp] = cp, void();
	int m = L + R >> 1; ls[cp] = ls[lp]; rs[cp] = rs[lp];
	if(st <= m) Mdf(ls[cp], ls[lp], L, m , st, std::min(ed, m), x);
	if(ed > m) Mdf(rs[cp], rs[lp], m + 1, R, std::max(st, m + 1), ed, x);
}
void Modify(int &rt, int l, int r, int v, int tm) {
	seg::Mdf(1, 1, n, l, r, v);
	Mdf(rt, rt, 1, n, l, r, tm);
}
int code(int x) {return (x + Lans * ty) % n + 1;}
int main() {
	n = ri(); m = ri(); ty = ri();
	for(int i = 1;i <= m; ++i) {
		int op = ri(); rt[i] = rt[i - 1];
		if(op == 2) {
			int x = code(ri());
			int tm = Qry(rt[i], 1, n, x);
			int nw = Qry(rt[tm - 1], 1, n, x);
			Modify(rt[i], x, x, u[nw], nw);
		}
		else {
			int l = code(ri()), r = code(ri());
			if(l > r) std::swap(l, r);
			if(op == 1)
				printf("%lld\n", Lans = seg::Qry(1, 1, n, l, r));
			else {
				u[i] = ri();
				Modify(rt[i], l, r, u[i], i);
			}
		}
	}
	return 0;
}