#include <bits/stdc++.h>
using namespace std;
typedef int (*uniopi)(int);
typedef long long (*uniopl)(long long);
typedef int (*binopi)(int, int);
typedef long long (*binopl)(long long, long long);
typedef int (*tupopi)(int, int, int);
typedef long long (*tupopl)(long long, long long, long long);
const long double PI = acos(0) * 2;
const int INF = 0x3f2f1f0f;
const long long LINF = 1ll * INF * INF;
int X4[] = {-1, 0, 1, 0}, Y4[] = {0, -1, 0, 1};
int X8[] = {-1, -1, -1, 0, 1, 1, 1, 0}, Y8[] = {-1, 0, 1, 1, 1, 0, -1, -1};
template <typename T>
inline bool uax(T& a, T b) {
  return a < b ? (a = b, 1) : 0;
}
template <typename T>
inline bool uin(T& a, T b) {
  return a > b ? (a = b, 1) : 0;
}
template <typename T1, typename T2>
istream& operator>>(istream& is, pair<T1, T2>& p) {
  return is >> p.first >> p.second;
}
template <typename T1, typename T2>
ostream& operator<<(ostream& os, const pair<T1, T2>& p) {
  return os << '(' << p.first << ", " << p.second << ")";
}
template <typename T>
ostream& operator<<(ostream& os, const vector<T>& v) {
  for (int i = 0; i < (int)v.size(); i++) {
    os << v[i];
    if (i + 1 < (int)v.size()) os << ' ';
  }
  return os;
}
void filll(int xs, int xf, int ys, int yf, int u, vector<vector<int> >& a) {
  int c = 0;
  for (int i = xs; i < xf; i++)
    for (int j = ys; j < yf; j++) a[i][j] = 4 * (c++) + u;
}
int main() {
  ios::sync_with_stdio(0);
  cin.tie(nullptr);
  cout << std::setprecision(10);
  cout << fixed;
  int n;
  cin >> n;
  int m = n / 2;
  vector<vector<int> > a(n, vector<int>(n, 0));
  filll(0, m, 0, m, 0, a);
  filll(0, m, m, n, 1, a);
  filll(m, n, 0, m, 2, a);
  filll(m, n, m, n, 3, a);
  for (auto& v : a) {
    cout << v << endl;
  }
  return 0;
}
