#include using namespace std; #define rep(i,n) for(ll i=0;i=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define pb push_back #define ins insert #define pqueue(x) priority_queue,greater> #define all(x) (x).begin(),(x).end() #define CST(x) cout<> #define rev(x) reverse(x); using ll=long long; using vl=vector; using vvl=vector>; using pl=pair; using vpl=vector; using vvpl=vector; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; const ll dy[9]={0,1,-1,0,1,1,-1,-1,0}; const ll dx[9]={1,0,0,-1,1,-1,1,-1,0}; template inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } struct SCC{//多重辺、自己辺ok int n,k=0; vector> g,rg; vector> ret; vector vs,cmp;//cmpは各頂点が縮約後のグラフでどの頂点に属するか、0~k-1の番号がトポロジカルソートして振られる。 vector used; SCC(int _n){ n=_n; g=vector>(n);rg=g; cmp=vector(n);used=vector(n); } void add_edge(int f, int t){ g[f].push_back(t);rg[t].push_back(f); } int init(){ rep(i,n)if(!used[i])dfs(i); used=vector(n,false); per(i,n){ if(!used[vs[i]]){ rdfs(vs[i],k); k++; } } return k; } vector> build(){//自己辺はないが多重辺あり ret.resize(k); rep(i,n){ for(auto to:g[i]){ if(cmp[i]==cmp[to])continue; ret[cmp[i]].push_back(cmp[to]); } } return ret; } private: void dfs(int v){ used[v]=true; rep(i,g[v].size()){ if(!used[g[v][i]])dfs(g[v][i]); } vs.push_back(v); } void rdfs(int v, int k){ used[v]=true;cmp[v]=k; rep(i,rg[v].size()){ if(!used[rg[v][i]]){ rdfs(rg[v][i],k); } } } }; struct twoSAT{ vector ans; twoSAT(int _n):n(_n),graph(2*_n),ans(_n){} void add_edge(int i,bool x,int j,bool y){//( i V j ),falseは否定 if(!x)i+=n;if(!y)j+=n; graph.add_edge((i+n)%(2*n),j); graph.add_edge((j+n)%(2*n),i); } bool exe(){ graph.init(); rep(i,n){ if(graph.cmp[i]==graph.cmp[i+n])return false; else if(graph.cmp[i]>graph.cmp[i+n])ans[i]=true; else ans[i]=false; } return true; } private: int n;SCC graph; }; class osak{ vector res; public: vector prime; osak(long long n){//sieve(n),O(NloglogN); res=vector(n); iota(res.begin(), res.end(), 0); for (int i = 2; i*i < n; ++i) { if (res[i] < i) continue; for (int j = i*i; j < n; j += i) if (res[j] == j) res[j] = i; } for(int i = 2; i < n; ++i){ if(res[i]==i)prime.emplace_back(i); } } map factor(int n) { map factor; while (n > 1) { factor[res[n]]++; n /= res[n]; } return factor; } vector divisor(int N){//O(div.size()) map facs=factor(N); vector ret={1}; for(auto p:facs){ ll range=ret.size(); ll now=1; for(int i=0;i> n; vector> a(n,vector(2)); rep(i,n)cin >> a[i][0] >> a[i][1]; osak os(1e6+10); twoSAT ts(n); /*rep(i,n){ ts.add_edge(i,false,i+n,false); }*/ rep(i,n){ rep(j,n){ rep(k,2){ rep(t,2){ ll q=stoll(a[i][k]+a[j][t^1]); if(binary_search(all(os.prime),q)){ ts.add_edge(i,k^1,j,t^1); } } } } } if(ts.exe())cout << "Yes" << endl; else cout << "No" << endl; }