#define _USE_MATH_DEFINES
#include <algorithm>
#include <cstdio>
#include <functional>
#include <iostream>
#include <cfloat>
#include <climits>
#include <cstring>
#include <cmath>
#include <map>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <time.h>
#include <vector>
using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> i_i;
typedef pair<ll, int> ll_i;
typedef pair<double, int> d_i;
typedef pair<ll, ll> ll_ll;
typedef pair<double, double> d_d;
struct edge { int u, v; ll w; };

ll MOD = 1000000007;
ll _MOD = 1000000009;
double EPS = 1e-10;

int a[4][4] = {
	{1, 0, 0, 1},
	{0, 1, 1, 0},
	{0, 1, 1, 0},
	{1, 0, 0, 1}
};

int main() {
	int N, X, Y, Z;
	cin >> N >> X >> Y >> Z;
	X--; Y--; Z--;
	vector<vector<int> > b(N, vector<int>(N));
	for (int y = 0; y < N; y++)
		for (int x = 0; x < N; x++)
			if (a[y % 4][x % 4])
				b[y][x] = y * N + x;
			else
				b[y][x] = N * N - 1 - y * N - x;
	int z = b[Y][X];
	for (int y = 0; y < N; y++) {
		for (int x = 0; x < N; x++) {
			b[y][x] ^= Z ^ z;
			cout << b[y][x] + 1 << ' ';
		}
		cout << endl;
	}
}