#include using namespace std; using ll=long long; const ll ILL=2167167167167167167; const int INF=2100000000; #define rep(i,a,b) for (int i=(int)(a);i<(int)(b);i++) #define all(p) p.begin(),p.end() template using _pq = priority_queue, greater>; template int LB(vector &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();} template int UB(vector &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();} template bool chmin(T &a,T b){if(b bool chmax(T &a,T b){if(a void So(vector &v) {sort(v.begin(),v.end());} template void Sore(vector &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});} bool yneos(bool a,bool upp=false){if(a){cout<<(upp?"YES\n":"Yes\n");}else{cout<<(upp?"NO\n":"No\n");}return a;} template void vec_out(vector &p,int ty=0){ if(ty==2){cout<<'{';for(int i=0;i<(int)p.size();i++){if(i){cout<<",";}cout<<'"'< T vec_min(vector &a){assert(!a.empty());T ans=a[0];for(auto &x:a) chmin(ans,x);return ans;} template T vec_max(vector &a){assert(!a.empty());T ans=a[0];for(auto &x:a) chmax(ans,x);return ans;} template T vec_sum(vector &a){T ans=T(0);for(auto &x:a) ans+=x;return ans;} int pop_count(long long a){int res=0;while(a){res+=(a&1),a>>=1;}return res;} template T square(T a){return a * a;} namespace po167{ template struct seg_sum{ int n; int seg_size; std::vector seg; T e; //要素数、単位元 typeは+演算が定義されているもの seg_sum(int k,T e_in){ e=e_in; n=k; seg_size=1; while(seg_size &q){ int n=q.size(); po167::seg_sum S(n,0); long long ans=0; for(int i=0;i //pのなかみは、重複ありで内容が同じならおけ //比較ができる必要がある。 long long inversion_vector(std::vector &p){ int n=p.size(); std::vector q(n); for(int i=0;i> t; rep(i, 0, t) solve(); } void solve(){ int N; cin >> N; vector A(N, vector(N)); rep(i, 0, N) rep(j, 0, N) cin >> A[i][j]; if (N < 4){ vector p(N * N); auto f = [&]() -> int { vector q; rep(i, 0, N * N) if (p[i]) q.push_back(p[i] - 1); return (po167::inversion_permutation(q)) & 1; }; rep(i, 0, N * N) p[i] = i; reverse(p.begin() + N, p.begin() + N * 2); auto base = f(); map> m; So(p); do{ if (f() == base){ ll val = 0; reverse(p.begin() + N, p.begin() + N + N); rep(i, 0, N) rep(j, 0, N) val += p[i * N + j] * A[i][j]; reverse(p.begin() + N, p.begin() + N + N); if (m.count(val)){ cout << "Yes\n"; auto out = [&](vector v) -> void { reverse(v.begin() + N, v.begin() + N + N); rep(i, 0, N){ rep(j, 0, N){ cout << v[i * N + j]; if (j + 1 == N) cout << "\n"; else cout << " "; } } }; out(m[val]); out(p); return; } m[val] = p; } }while (next_permutation(all(p))); cout << "No\n"; return; } ll Z = 12; vector num; vector X; vector ans(N, vector(N)); vector> pos; rep(i, 0, N){ if (i & 1){ for (int j = N - 1; j >= 0; j--){ ans[i][j] = i * N + j; if (ans[i][j] < Z) pos.emplace_back(i, j), num.push_back(ans[i][j]), X.push_back(A[i][j]); } } else{ rep(j, 0, N){ ans[i][j] = i * N + j; if (ans[i][j] < Z) pos.emplace_back(i, j), num.push_back(ans[i][j]), X.push_back(A[i][j]); } } } vector p; auto f = [&]() -> int { vector q; rep(i, 0, Z) if (p[i]) q.push_back(p[i] - 1); return (po167::inversion_permutation(q)) & 1; }; p = num; auto base = f(); unordered_map> m; So(p); do{ if (f() == base){ ll val = 0; rep(i, 0, Z) val += X[i] * p[i]; if (m.count(val)){ cout << "Yes\n"; auto out = [&](vector v) -> void { rep(i, 0, Z) ans[pos[i].first][pos[i].second] = v[i]; rep(i, 0, N) vec_out(ans[i]); }; out(m[val]); out(p); return; } m[val] = p; } }while (next_permutation(all(p))); cout << "No\n"; return; }