#include #include #include using namespace std; using ll = long long; #define rep(i, s, t) for (ll i = s; i < (ll)(t); i++) #define all(x) begin(x), end(x) template bool chmin(T& x, T y) { return x > y ? (x = y, true) : false; } template bool chmax(T& x, T y) { return x < y ? (x = y, true) : false; } struct IOST { IOST() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(20); } } IOST; int timer() { static auto st = chrono::system_clock::now(); auto en = chrono::system_clock::now(); return chrono::duration_cast(en - st).count(); } void solve() { timer(); int n; cin >> n; vector> a(n, vector(n)); rep(i, 0, n) rep(j, 0, n) cin >> a[i][j]; int m = min(4, n); auto print_p = [&](vector p) { if ((int)p.size() < n * n) { p.resize(n * n); rep(i, m * m, n * n) { p[i] = i; } } rep(i, 0, n) { rep(j, 0, n) cout << p[i * n + j] << " "; cout << "\n"; } }; auto check_p = [&](vector p) { int ad = 0; rep(i, 0, m * m) { if (p[i] == 0) { ad = (i / n + i % n) & 1; } } vector rp(m * m); rep(i, 0, m * m) rp[p[i]] = i; rep(i, 0, m * m) { if (p[i] != i) { swap(p[rp[i]], p[i]); rp[p[rp[i]]] = rp[i]; rp[i] = i; ad = (ad + 1) & 1; } } return ad; }; vector p(m * m); iota(p.begin(), p.end(), 0); map> mp; do { if (check_p(p)) continue; ll tmp = 0; rep(i, 0, m) { rep(j, 0, m) tmp += a[i][j] * p[i * n + j]; } if (mp.count(tmp)) { cout << "Yes\n"; print_p(p); print_p(mp[tmp]); return; } mp[tmp] = p; if (timer() > 3500) break; } while (next_permutation(p.begin(), p.end())); cout << "No\n"; } int main() { int t = 1; // cin >> t; rep(i, 0, t) solve(); }