#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define fi first #define se second #define rep(i,n) for(int i = 0; i < n; ++i) #define rrep(i,n) for(int i = 1; i <= n; ++i) #define drep(i,n) for(int i = n-1; i >= 0; --i) #define gep(i,g,j) for(int i = g.head[j]; i != -1; i = g.e[i].next) #define each(it,c) for(__typeof((c).begin()) it=(c).begin();it!=(c).end();it++) #define rng(a) a.begin(),a.end() #define maxs(x,y) x = max(x,y) #define mins(x,y) x = min(x,y) #define pb push_back #define sz(x) (int)(x).size() #define pcnt __builtin_popcount #define snuke srand((unsigned)clock()+(unsigned)time(NULL)); using namespace std; typedef long long int ll; typedef pair P; typedef vector vi; typedef vector vvi; inline int in() { int x; scanf("%d",&x); return x;} inline void priv(vi a) { rep(i,sz(a)) printf("%d%c",a[i],i==sz(a)-1?'\n':' ');} const int MX = 100005, INF = 1000010000; const ll LINF = 1000000000000000000ll; const double eps = 1e-10; int n, s; vvi a; int calc() { int res = 0; rep(i,n) { int x = 0; rep(j,n) x += a[i][j]; res += abs(x-s); x = 0; rep(j,n) x += a[j][i]; res += abs(x-s); } { int x = 0; rep(i,n) x += a[i][i]; res += abs(x-s); } { int x = 0; rep(i,n) x += a[n-1-i][i]; res += abs(x-s); } return res; } int main(){ scanf("%d",&n); a = vvi(n,vi(n)); rep(i,n)rep(j,n) a[i][j] = i*n+j+1; int sum = 0; rep(i,n)rep(j,n) sum += a[i][j]; s = sum/n; int score = calc(); rep(ti,MX) { vvi p = a; int ai = rand()%n; int aj = rand()%n; int bi = rand()%n; int bj = rand()%n; swap(a[ai][aj], a[bi][bj]); int sc = calc(); if (sc > score) { a = p; } else { score = sc; // cerr << score << endl; // rep(i,n) priv(a[i]); } } rep(i,n) priv(a[i]); return 0; }