#include<bits/stdc++.h>
using namespace std;
const int N=2005;
int f[N][N],h[N];
int main(){
    int n;cin>>n;
    cout<<"Yes\n";
    int l=!n?1:__lg(n-1)+1;
    for(int i=0;i<n;i++){
        for(int j=0;j<l;j++) h[i+1]|=((i>>j&1)^(i>>j+1&1))<<j;
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)if(i+j+1&1){
        f[i][j]=h[i+j>>1]<<l|h[abs(i-j)+2>>1];
    }
    int ful=(1<<2*l)-1;
    for(int i=1;i<=n;i++,cout<<'\n')
    for(int j=1;j<=n;j++){
        if(i+j&1) cout<<((f[i][j-1]|f[i-1][j]|f[i+1][j]|f[i][j+1])^ful)<<' ';
        else cout<<f[i][j]<<' ';
    }
}