#include #include using namespace std; class range {private: struct I{int x;int operator*(){return x;}bool operator!=(I& lhs){return x> &res) { int i = 0, j = n / 2; for(int k : range(n*n)) { res[i][j] = k + 1; int ni = (i - 1 + n) % n; int nj = (j + 1) % n; if(res[ni][nj] > 0) { ni = (i + 1) % n; nj = j; } i = ni, j = nj; } } void solve2(int n, vector> &res) { int blocks = n / 4; for(int k : range(n*n)) { int i = k / n; int j = k % n; if((i + j) % 4 == 3 || (i - j + 4) % 4 == 0) { res[i][j] = k + 1; } } for(int k : range(n*n)) { int i = n-1 - k / n; int j = n-1 - k % n; if((i + j) % 4 != 3 && (i - j + 4) % 4 != 0) { res[i][j] = k + 1; } } } // (ΦωΦ)<関孝和の方法 void solve3(int n, vector> &res) { int n2 = n - 2; vector> inner(n2, vector(n2)); solve2(n2, inner); for(int i : range(n2)) { for(int j : range(n2)) { res[i+1][j+1] = inner[i][j] + 2*n - 2; } } for(int k : range(n2)) { res[n-1][k+2] = k + 1; } for(int k : range(n2)) { res[n-1-k][0] = n2 + k + 1; } res[n-1][1] = 2 * n2 + 1; res[1][0] = 2 * n2 + 2; for(int k : range(n)) { int target = (k==0||k==n-1) ? res[n-1][n-1-k] : res[n-1][k]; res[0][k] = n*n+1 - target; } for(int k : range(1, n)) { int target = (k==n-1) ? res[0][0] : res[k][0]; res[k][n-1] = n*n+1 - target; } swap(res[0][1], res[n-1][1]); for(int j : range(2, n-1)) { if((j-2)%4 < 2) { swap(res[0][j], res[n-1][j]); } } for(int i : range(1, n-1)) { if(i%4 < 2) { swap(res[i][0], res[i][n-1]); } } } int main(void) { int n; scanf("%d", &n); vector> res(n, vector(n)); if(n & 1) { solve1(n, res); } else if(n % 4 == 0) { solve2(n, res); } else { solve3(n, res); } for(int i : range(n)) { for(int j : range(n)) { printf("%d%c", res[i][j], j==n-1?'\n':' '); } } return 0; }