N, X, Y, Z = map(int, raw_input().split())

diag = set([(0, 0), (0, 3), (1, 1), (1, 2), (2, 1), (2, 2), (0, 3), (3, 3)])

for d in xrange(N * N):
        A = [[0] * N for i in xrange(N)]
        for r in xrange(N):
                for w in xrange(N):
                        if (w % 4, r % 4) in diag:
                                A[r][w] = (w + N * r + d) % (N * N) + 1
        for r in xrange(N):
                for w in xrange(N):
                        if A[N - r - 1][N - w - 1] == 0:
                                A[N - r - 1][N - w - 1] = (w + N * r + d) % (N * N) + 1
        if A[Y - 1][X - 1] == Z:
                break
for line in A:
        print " ".join(map(str, line))