#include <bits/stdc++.h>
using namespace std;
#define int long long
#ifndef GammaRays
#define debug //
#define endl "\n"
#endif
const int N = 3e5 + 10;
const int mod = 998244353;
void solve() {
	int n;
	cin >> n;
	
	vector ans(n + 1, vector<int>(n + 1));
	for (int i = 1; i <= n - 2; i ++) {
		for (int j = 1; j <= n; j ++) ans[i][j] = i;
	}
	for (int i = n - 1; i <= n; i ++) {
		for (int j = 1; j <= n - 2; j ++) ans[i][j] = n - 2 + j;
	}
	
	ans[n - 1][n - 1] = 2 * n - 3;	
	ans[n - 1][n] = 2 * n - 2;	
	ans[n][n - 1] = 2 * n - 1;	
	ans[n][n] = 2 * n;
	
	cout << "Yes" << '\n';
	for (int i = 1; i <= n; i ++) {
		for (int j = 1; j <= n; j ++) {
			cout << ans[i][j] << " \n"[j == n];
		}
	}
}
signed main() {
    cin.tie(nullptr)->sync_with_stdio(false);
    int t = 1;
    // cin >> t;
    while (t --) {
        solve();
    }
    return 0;
}