N,X,Y,Z=map(int,raw_input().split())
sq=[[-1 for i in range(N)] for j in range(N)]
cur=0
for x in range(N):
    for y in range(N):
        a,b=x%4,y%4
        if a==b or a+b==3:
            sq[x][y]=cur
        cur+=1
cur=0
for x in range(N-1,-1,-1):
    for y in range(N-1,-1,-1):
        if sq[x][y]==-1:
            sq[x][y]=cur
        cur+=1

dif= (Z-1)^sq[Y-1][X-1]
for x in range(N):
    for y in range(N):
        sq[x][y]^=dif
        print sq[x][y]+1,
    print