def solve():
    N=int(input())

    if N==1:
        print(-1)
        return

    X=[[7,14,0,8],[4,12,2,11],[15,9,6,1],[13,10,5,3]]
    Y=[[0,1,2,3],[4,5,6,7],[8,9,10,11],[12,13,14,15]]

    A=[[-1]*(1<<N) for _ in range(1<<N)]
    m=0
    for s in range(pow(2,N)//4):
        for t in range(pow(2,N)//4):
            if s==t:
                for i in range(4):
                    for j in range(4):
                        A[4*s+i][4*t+j]=X[i][j]^m
            else:
                for i in range(4):
                    for j in range(4):
                        A[4*s+i][4*t+j]=Y[i][j]^m
            m+=16

    for Ai in A:
        print(*Ai)
    return

#==================================================
import sys
input=sys.stdin.readline
write=sys.stdout.write

solve()