#include<bits/stdc++.h>
using namespace std;
#define ll long long
map<int,int>mat;
int n,p[10005],vis[10005];
ll a[505][505];
ll gcd(ll x,ll y){
	if (!y)return x;
	return gcd(y,x%y);
}
ll lcm(ll x,ll y){
	if ((!x)||(!y))return x+y;
	return x/gcd(x,y)*y;
}
int nex(int x){
	if (!mat[x])mat[x]=p[p[0]--];
	return mat[x];
}
int main(){
	scanf("%d",&n);
	if (n==2){
		printf("4 7\n23 10\n");
		return 0;
	}
	for(int i=2;i<=10000;i++){
		if (!vis[i])p[++p[0]]=i;
		for(int j=1;(j<=p[0])&&(i*p[j]<=10000);j++){
			vis[i*p[j]]=1;
			if (i%p[j]==0)break;
		}
	}
	for(int i=1;i<=p[0]/2;i++)swap(p[i],p[p[0]-i+1]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if ((i+j)%2==0)a[i][j]=nex(i+j);
	mat.clear();
	for(int i=n;i;i--)
		for(int j=n;j;j--)
			if ((i+j)%2==0)a[i][j]*=nex(i-j);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if ((i+j)%2)a[i][j]=lcm(lcm(a[i-1][j],a[i][j-1]),lcm(a[i+1][j],a[i][j+1]))+1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)printf("%lld ",a[i][j]);
		printf("\n");
	}
}