JZOJ4605. 排序（线段树合并与分裂）

题目大意：

---

每次把一个区间升序或降序排序，最后问一个点是什么。

题解：

---

如果只是问一个点，这确乎是个经典题，二分一下答案然后线段树维护01排序。

从pty那里get到了可以用线段树的合并与分裂实时地维护整个序列。

考虑一次排序就把这个区间的数搞到一个线段树上，在根处标记是正的还是反的。

如果想搞到一棵树上就需要用到分裂与合并，根据势能分析，复杂度还是\(O(n~log~n)\)。

Code：

#include<bits/stdc++.h>
#define fo(i, x, y) for(int i = x, B = y; i <= B; i ++)
#define ff(i, x, y) for(int i = x, B = y; i <  B; i ++)
#define fd(i, x, y) for(int i = x, B = y; i >= B; i --)
#define ll long long
#define pp printf
#define hh pp("\n")
using namespace std;

const int N = 1e5 + 5;

int n, m, a[N];
int op, l, r;

struct tree {
	int l, r, x;
} t[N * 120];
#define i0 t[i].l
#define i1 t[i].r
int tot, pl, pr, px;
int rt[N * 120], rev[N * 120];

void upd(int i) {
	t[i].x = t[i0].x + t[i1].x;
}

void add(int &i, int x, int y) {
	i = ++ tot;
	t[i].x = 1;
	if(x == y) return;
	int m = x + y >> 1;
	if(pl <= m) add(i0, x, m); else add(i1, m + 1, y);
}

void merge(int &i, int j) {
	if(!i || !j) { i = i + j; return;}
	merge(t[i].l, t[j].l);
	merge(t[i].r, t[j].r);
	upd(i);
}

void ft(int i, int x, int y) {
	if(x == y) { if(t[i].x == pl) px = x; return;}
	int m = x + y >> 1;
	if(t[i0].x >= pl) ft(i0, x, m); else
		pl -= t[i0].x, ft(i1, m + 1, y);
}

void sp(int i, int &j1, int &j2, int x, int y, int k) {
	if(!i) return;
	if(y <= k) { j1 = i, j2 = 0; return;}
	if(x > k)  { j1 = 0, j2 = i; return;}
	int m = x + y >> 1;
	if(!j1) j1 = ++ tot;
	if(!j2) j2 = ++ tot;
	sp(i0, t[j1].l, t[j2].l, x, m, k);
	sp(i1, t[j1].r, t[j2].r, m + 1, y, k);
	upd(j1); upd(j2);
}

int bz[N];
void dg(int i, int x, int y) {
	if(!t[i].x) return;
	if(x == y) { bz[x] ++; return;}
	int m = x + y >> 1;
	dg(i0, x, m); dg(i1, m + 1, y);
}

void split(int &i, int &j, int k) {
	pl = rev[i] ? t[i].x - k : k;
	px = 0; ft(i, 0, n);
	if(rev[i]) sp(i, j, i, 0, n, px); else sp(i, i, j, 0, n, px);
}

set<int> s;
set<int> :: iterator st, en;
int d[N], d0;

int main() {
	scanf("%d %d", &n, &m);
	fo(i, 1, n) scanf("%d", &a[i]);
	fo(i, 0, n + 1) s.insert(i);
	fo(i, 0, n) pl = pr = a[i], add(rt[i], 0, n);
	fo(i, 1, m) {
		scanf("%d %d %d", &op, &l, &r);
		en = s.upper_bound(r);
		st = -- s.lower_bound(l);
		d0 = 0;
		for(; st != en; st ++) d[++ d0] = *st;
		d[++ d0] = *en;
		int rt1 = 0, rt2 = 0;
		split(rt[d[1]], rt1, l - d[1]);
		if(t[rt1].x) {
			rt[l] = rt1;
			rev[rt[l]] = rev[rt[d[1]]];
			d[1] = l;
		}
		if(d[d0] != r + 1) {
			split(rt[d[d0 - 1]], rt2, r + 1 - d[d0 - 1]);
			rt[r + 1] = rt2;
			rev[rt[r + 1]] = rev[rt[d[d0 - 1]]];
		}
		fo(j, 2, d0 - 1) {
			if(l != d[j]) merge(rt[l], rt[d[j]]), rt[d[j]] = 0;
			s.erase(d[j]);
		}
		rev[rt[l]] = op;
		s.insert(l); s.insert(r + 1);
	}
	scanf("%d", &l);
	l ++;
	for(st = s.begin(); st != s.end(); st ++) {
		int i = *st;
		if(t[rt[i]].x >= l) {
			pl = rev[rt[i]] ? t[rt[i]].x - l + 1 : l;
			px = 0;
			ft(rt[i], 0, n);
			pp("%d\n", px);
			return 0;
		} else l -= t[rt[i]].x;
	}
}