【题解】【LibreOJ Beta Round #5】游戏 LOJ 531 基环树 博弈论

Prelude

题目链接：萌萌哒传送门♪(^∇^*)

---

Subtask 1 & 2

这是什么鬼题面。。。
首先要看出，这就是一个基环树博弈。
具体题意：给出一个基环内向树，一个棋子初始在(1)号节点，双方轮流操作，设棋子所在节点为(u)，每次可以从所有指向(u)的节点中选择一个，把棋子移动过去，不能操作者输，问先手是否有必胜策略，或者是否平局。
但是，这个图是可以变的，每次询问，假如我选择环上的两个点(u)和(v)，把(u)的出边指向(v)，那么游戏的结果如何？
考虑如果只有一次询问，不修改的话怎么做。
这么裸的博弈随便做吧？先假设把环上的所有边断开，然后dp求出每个节点的状态，最后再考虑环上的情况，当环上所有点都是必败态的时候就是平局，否则两个人沿着环走，第一个走到必胜节点的人获胜。
求解一次是(O(n))的，一共(q)次询问，总时间复杂度(O(nq))。

---

Subtask 3 & 4

首先我们可以预处理求出离环上每个点最近的必胜点是哪个，我这里用(left_{i})表示。
每次对原图做修改，相当于把环缩小到原来的环的一个片段。
假如我们选了环上的(u)和(v)两个节点，并且把(u)的出边连接到(v)，那么(f_{u})到(v)这一段的dp值都可能会改变。
我们可以用倍增预处理，用(jump[x][s][k])表示假如点(x)的状态为(s)，那么(x)之后第(2^{k})个点的状态是什么，这样，查询一个点的状态就是(O(log n))的了。
对于点(1)不在环上的情况，我们可以直接查询，(O(log n))解决。
如果点(1)在环上，我们因为有之前预处理的结果，所以只需要求出点(v)的新的dp值（(O(log n))），然后大力分类讨论（(O(1))）解决。
总体复杂度就是(O(n log n))，能过87分辣~

---

Subtask 5

LCA告诉我可以线性做，然而我并不会。
瓶颈是倍增，怎么去掉倍增？怎么降低查询的复杂度？
我们离线所有的询问，这样就可以用并查集，在(O(n alpha(n)))的时间内处理出所有的“对点(v)的dp值的询问”。
这样总体复杂度就变成了(O(n alpha(n)))辣~
因为我比较懒，所以离线询问的时候对询问排了序，复杂度还是(O(n log n))的，但是这个(O(n log n))比倍增的(O(n log n))要小得多了，也是能AC的。

---

Conclusion

为了处理方便，我们通常把环拆开，拆成一条链来处理。
于是就少不了大量的分类讨论，我写了200行。
这件事告诉我们，LCA题不可做(╯‵□′)╯︵┻━┻。

---

Code

#include <cstring>
#include <algorithm>
#include <cstdio>
#include <queue>
#include <cctype>

using namespace std;
const int MAXN = 1000010;
const int LOGN = 21;
int _w;

int read() {
	int x = 0, ch;
	while( isspace(ch = getchar()) );
	do x = x * 10 + ch - '0';
	while( isdigit(ch = getchar()) );
	return x;
}

int n, q, f[MAXN];

int lp[MAXN], lpsz, pos[MAXN]; // lp为环上的每个节点，lpsz为环的大小，pos为原图中的节点在环中的位置
void find_loop() {
	int u = 1;
	while( !pos[u] ) {
		pos[u] = ++lpsz;
		lp[lpsz] = u;
		u = f[u];
	}
	int sz = 0;
	for( int i = pos[u]; i <= lpsz; ++i )
		lp[++sz] = lp[i];
	lpsz = sz;
	for( int i = 1; i <= n; ++i ) pos[i] = 0; // 不在环中的节点pos == 0
	for( int i = 1; i <= lpsz; ++i )
		pos[lp[i]] = i;
}

int pa[MAXN<<1]; // 并查集相关
int find( int u ) {
	return pa[u] == u ? u : pa[u] = find(pa[u]);
}

int dp[MAXN], deg[MAXN], left[MAXN];
queue<int> bfsq;
int toend[MAXN][2], to1[MAXN][2]; // to1[x][s]表示当x点的状态为s的时候，点1的状态是什么，toend类似，表示序列上最后一个点的信息
void prelude() {
	find_loop();
	for( int i = 1; i <= n; ++i )
		++deg[f[i]];
	for( int i = 1; i <= n; ++i )
		if( !deg[i] )
			bfsq.push(i);
	while( !bfsq.empty() ) {
		int u = bfsq.front(); bfsq.pop();
		if( !dp[u] ) dp[f[u]] = 1;
		if( --deg[f[u]] == 0 )
			bfsq.push(f[u]);
	}
	for( int i = 1; i <= lpsz; ++i )
		if( dp[lp[i]] ) left[i] = i;
		else left[i] = left[i-1];
	for( int i = 1; i < lpsz; ++i ) {
		pa[i] = i+1+lpsz;
		pa[i+lpsz] = i+1 + dp[lp[i+1]] * lpsz;
	}
	pa[lpsz] = lpsz;
	pa[lpsz+lpsz] = lpsz+lpsz;
	for( int i = 1; i <= lpsz; ++i ) {
		toend[i][0] = find(i) > lpsz;
		toend[i][1] = find(i+lpsz) > lpsz;
	}
	if( pos[1] ) {
		for( int i = 1; i != pos[1]; ++i ) {
			pa[i] = i+1+lpsz;
			pa[i+lpsz] = i+1 + dp[lp[i+1]] * lpsz;
		}
		pa[pos[1]] = pos[1];
		pa[pos[1]+lpsz] = pos[1]+lpsz;
		for( int i = 1; i != pos[1]; ++i ) {
			to1[i][0] = find(i) > lpsz;
			to1[i][1] = find(i+lpsz) > lpsz;
		}
		to1[pos[1]][0] = 0;
		to1[pos[1]][1] = 1;
	}
	for( int i = 1; i <= lpsz; ++i ) {
		pa[i] = i;
		pa[i+lpsz] = i+lpsz;
	}
}

int calc( int t, int v, int p ) {
	static int right = 1;
	if( t > p ) { // 把环拆成序列就要大力分类讨论。。。
		v = toend[t][v];
		t = 1;
		v = !v || dp[lp[1]];
	}
	if( p == pos[1] )
		return to1[t][v];
	while( right < p ) { // 对所有询问按照查询位置排序
		pa[right] = right+1+lpsz;
		pa[right+lpsz] = right+1 + dp[lp[right+1]] * lpsz;
		++right;
	}
	return find(t + v*lpsz) > lpsz;
}

int solve( int l, int r ) {
	if( !pos[1] ) return dp[1]; // 点1不在环上
	l = pos[l], r = pos[r];
	int p = pos[1];
        if( l >= r ) { // 下面都是分类讨论。。。
		swap(l, r);
		if( p < l || p > r ) { // 点1不在环上
			int t = r == lpsz ? 1 : r+1;
			return calc(t, dp[lp[t]], p);
		} else {
			int dpl = dp[lp[l]];
			if( l > 1 || r < lpsz ) {
				int t = r == lpsz ? 1 : r+1;
				dpl = calc(t, dp[lp[t]], l);
			}
			if( left[p] >= l+1 ) {
				int q = left[p];
				return (p-q+1)&1;
			} else if( dpl ) {
				return (p-l+1)&1;
			} else if( left[r] >= l+1 ) {
				int q = left[r];
				return (p-l+1+r-q+1)&1;
			} else {
				return 2;
			}
		}
	} else {
		if( p > l && p < r ) {
			int t = l+1;
			return calc(t, dp[lp[t]], p);
		} else if( p <= l ) {
			int t = l+1;
			int dpr = calc(t, dp[lp[t]], r);
			if( left[p] ) {
				int q = left[p];
				return (p-q+1)&1;
			} else if( left[lpsz] >= r+1 ) {
				int q = left[lpsz];
				return (lpsz-q+1+p)&1;
			} else if( dpr ) {
				return (p+lpsz-r+1)&1;
			} else if( left[l] ) {
				int q = left[l];
				return (lpsz-r+1+p+l-q+1)&1;
			} else {
				return 2;
			}
		} else {
			int t = l+1;
			int dpr = calc(t, dp[lp[t]], r);
			if( left[p] >= r+1 ) {
				int q = left[p];
				return (p-q+1)&1;
			} else if( dpr ) {
				return (p-r+1)&1;
			} else if( left[l] ) {
				int q = left[l];
				return (p-r+1+l-q+1)&1;
			} else if( left[lpsz] >= r+1 ) {
				int q = left[lpsz];
				return (p-r+1+l+lpsz-q+1)&1;
			} else {
				return 2;
			}
		}
	}
}

int qu[MAXN], qv[MAXN]; // 存放询问的地方
int rk[MAXN], ans[MAXN];

bool cmp_qv( int i, int j ) { // 对询问按照点在序列上的位置排序
	i = pos[qv[i]], j = pos[qv[j]];
	return i < j;
}

int main() {
	n = read(), q = read();
	for( int i = 1; i <= n; ++i )
		f[i] = read();
	prelude();
	for( int i = 0; i < q; ++i )
		qu[i] = read(), qv[i] = read(), rk[i] = i;
	sort(rk, rk+q, cmp_qv);
	for( int i = 0; i < q; ++i )
		ans[rk[i]] = solve(qu[rk[i]], qv[rk[i]]);
	for( int i = 0; i < q; ++i )
		printf( "%d
", ans[i] );
	return 0;
}