#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define MN 500
#define MX 5000
int u[MX+5],p[MX+5],pn;
ll a[MN+5][MN+5];
map<ll,int> mp;
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
int main()
{
	int n,i,j,k,x,y;
	for(i=2;i<=MX;++i)
	{
		if(!u[i])p[++pn]=i;
		for(j=1;i*p[j]<=MX;++j)if(u[i*p[j]]=1,i%p[j]==0)break;
	}
	for(i=2;i<=MN+MN;i+=2)for(j=-MN;j<=MN;j+=2)
	{
		x=(i+j)/2;y=(i-j)/2;
		if(x<1||x>MN||y<1||y>MN)continue;
		for(k=1;mp[2*p[i/2]*p[(j+MN)/2]*k]++;++k);
		a[x][y]=2*p[i/2]*p[(j+MN)/2]*k;
	}
	for(i=1;i<=MN;++i)for(j=1;j<=MN;++j)if(!a[i][j])
	{
		ll x=1;
		if(i>1)x=x/gcd(x,a[i-1][j])*a[i-1][j];
		if(i<MN)x=x/gcd(x,a[i+1][j])*a[i+1][j];
		if(j>1)x=x/gcd(x,a[i][j-1])*a[i][j-1];
		if(j<MN)x=x/gcd(x,a[i][j+1])*a[i][j+1];
		for(k=1;mp[x*k+1]++;++k);
		a[i][j]=x*k+1;
	}
	scanf("%d",&n);
	for(i=1;i<=n;++i,puts(""))for(j=1;j<=n;++j)printf("%lld ",a[i][j]);
}