#include typedef long long ll; typedef unsigned long long ull; #define FOR(i,a,b) for(int (i)=(a);i<(b);i++) #define REP(i,n) FOR(i,0,n) #define RANGE(vec) (vec).begin(),(vec).end() using namespace std; class SeriesAmidakuji { public: bool dfs(vector &s, const vector &used, const vector &A) { int N = s.size(); int n = count_if(RANGE(s), [&](int v){return (v>=0);}); if (n==N) { REP(i,N) // 整合性チェック if (A[i] != s[s[i]]) return false; return true; } int beg; REP(i,N) { if (s[i] < 0) { beg = i; break; } } // O(N^2) REP(i,N) { if (used[i]) continue; vector w(used); bool ok = true; // ii -> i に決め打ちしてためしてみる if (beg == i) // s[beg] = beg のとき { s[beg] = i; if (s[s[beg]] != A[beg]) continue; w[i] = true; } else { int cur = beg; s[cur] = i; w[i] = true; do { if (s[cur] == A[cur]) { ok = false; break; } s[s[cur]] = A[cur]; w[A[cur]] = true; cur = s[cur]; } while (cur != beg); } if (ok && dfs(s, w, A)) return true; } return false; } void solve(void) { int N; cin>>N; vector A(N); REP(i,N) { cin>>A[i]; --A[i]; } // 1,...,N の置換 s のうち // s^2 = A なるものを見つける // // 置換は部分置換の直積として表現できる。 // s = s1*s2*...*sr // 部分置換の個数は最大 N なので全探索でできる。 // vector s(N,-1); vector used(N,false); if (dfs(s, used, A)) cout<<"Yes"<solve(); delete obj; return 0; } #endif