ARC079F - Namori Grundy(构造,基环树)
ARC079F - Namori Grundy
Solution
首先這是一個NNN個點NNN條邊的有向圖,所以它的基圖是一棵基環樹,其次這個圖的所有點入度為111,因此這是一棵基環外向樹。
然后對于aia_iai?,假設我們求出S={aj∣(i,j)∈E}S=\{a_j|(i,j)\in E\}S={aj?∣(i,j)∈E},即iii的所有出邊的aaa的集合,那么顯然ai=mexSa_i=mex\;Sai?=mexS,aia_iai?的值是可以通過其出邊唯一確定的。
我們先考慮一棵樹的情況,我們發現葉子結點必然為000,因此每一個結點iii的aia_iai?都可以從下到上通過其兒子結點遞推得到。
現在考慮基環樹,對于環上一點,我們可以通過它的子樹唯一確定去掉環時它的aaa,然后可以發現當且僅當環長為奇數且環上結點的aaa都相等時無解。
然后直接找到環,遞推求出每一個aaa,判一判就行了。
時間復雜度O(n)O(n)O(n)。
Proof
現在我們可以把子樹點都扔掉,只考慮一個kkk元環。
我們設環上點依次為0,1,2...k?10,1,2...k-10,1,2...k?1,(i+1,i)∈E(i+1,i)\in E(i+1,i)∈E,(之后的下標都在模kkk意義下)。
首先我們不會讓所有點都加111,因為這和所有點不變一樣,所以若存在方案,一定有一種有至少一個點ppp不動的方案。
如果我們知道ppp,那么我們可以從ppp開始判斷apa_pap?是否等于ap+1a_{p+1}ap+1?,若相等則ap+1+1a_{p+1}+1ap+1?+1。
我們令一段連續的+1+1+1的下標區間為[l+1,r][l+1,r][l+1,r],顯然有al=al+1,al+1+1=al+2,al+2+1=al+3...a_{l}=a_{l+1},a_{l+1}+1=a_{l+2},a_{l+2}+1=a_{l+3}...al?=al+1?,al+1?+1=al+2?,al+2?+1=al+3?...,因此這段aaa一定是al,al,al+1,al+2...al+(r?l?1)a_l,a_l,a_l+1,a_l+2...a_l+(r-l-1)al?,al?,al?+1,al?+2...al?+(r?l?1),點lll不動且為其中的最小值。
所以我們有結論:若能找到一個iii,使得ai=?ai?1a_i\not = a_{i-1}ai??=ai?1?且?jai≤aj\forall_j a_i\leq a_j?j?ai?≤aj?,iii可以作為ppp,方案存在。
那么最后還剩下找不到iii的環,其上的點必然是所有aaa相等的,這個就很容易考慮了,因為所有點都等價,所以可以任取一個作為ppp,最后的序列一定長成ap,ap+1,ap,ap+1,ap,ap+1...a_p,a_p+1,a_p,a_p+1,a_p,a_{p}+1...ap?,ap?+1,ap?,ap?+1,ap?,ap?+1...的形式。
這就相當于一個二分圖染色,只有當kkk為奇數時,存在相鄰兩個aaa相等,不合法;當aaa為偶數時存在合法方案。
綜上,當且當且僅當環長為奇數且環上結點的aaa都相等時無解。
Code
#include <vector> #include <list> #include <map> #include <set> #include <deque> #include <queue> #include <stack> #include <bitset> #include <algorithm> #include <functional> #include <numeric> #include <utility> #include <sstream> #include <iostream> #include <iomanip> #include <cstdio> #include <cmath> #include <cstdlib> #include <cctype> #include <string> #include <cstring> #include <ctime> #include <cassert> #include <string.h> //#include <unordered_set> //#include <unordered_map> //#include <bits/stdc++.h>#define MP(A,B) make_pair(A,B) #define PB(A) push_back(A) #define SIZE(A) ((int)A.size()) #define LEN(A) ((int)A.length()) #define FOR(i,a,b) for(int i=(a);i<(b);++i) #define fi first #define se secondusing namespace std;template<typename T>inline bool upmin(T &x,T y) { return y<x?x=y,1:0; } template<typename T>inline bool upmax(T &x,T y) { return x<y?x=y,1:0; }typedef long long ll; typedef unsigned long long ull; typedef long double lod; typedef pair<int,int> PR; typedef vector<int> VI;const lod eps=1e-11; const lod pi=acos(-1); const int oo=1<<30; const ll loo=1ll<<62; const int mods=998244353; const int MAXN=600005; const int INF=0x3f3f3f3f;//1061109567 /*--------------------------------------------------------------------*/ inline int read() {int f=1,x=0; char c=getchar();while (c<'0'||c>'9') { if (c=='-') f=-1; c=getchar(); }while (c>='0'&&c<='9') { x=(x<<3)+(x<<1)+(c^48); c=getchar(); }return x*f; } vector<int> e[MAXN]; int fa[MAXN],vis[MAXN],instk[MAXN],flag[MAXN],tag[MAXN],f[MAXN],n; void dfs(int x) {vis[x]=1,instk[x]=1;for (auto v:e[x]){if (instk[v]) {for (int p=x;p!=v;p=fa[p]) flag[p]=1;flag[v]=1;}else fa[v]=x,dfs(v);}instk[x]=0; }void tree_dp(int x,int father) {for (auto v:e[x]) if (v!=father&&!flag[v]) tree_dp(v,x);for (auto v:e[x]) if (v!=father&&!flag[v]) tag[f[v]]=x;for (int i=0;i<n;i++) if (tag[i]!=x) { f[x]=i; return; } }signed main() {n=read();for (int i=1,x;i<=n;i++) x=read(),e[x].PB(i);for (int i=1;i<=n;i++) if (!vis[i]) dfs(i);int p=0,num=0;for (int i=1;i<=n;i++) if (flag[i]) tree_dp(i,0),p=f[i],num++;if (!(num&1)) { puts("POSSIBLE"); return 0; } for (int i=1;i<=n;i++)if (flag[i]&&p!=f[i]) { puts("POSSIBLE"); return 0; }puts("IMPOSSIBLE");return 0; }總結
以上是生活随笔為你收集整理的ARC079F - Namori Grundy(构造,基环树)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: ARC078F - Mole and A
- 下一篇: ARC080F - Prime Flip