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

时间:2024-01-21 11:28:28

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\n", ans[i] );
return 0;
}