import sys,time

from itertools import permutations
from heapq import heappop,heappush
from collections import deque
import random
import bisect
from math import log,gcd

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

def find_euler_tour(N,M,edge):
    done_counter = [0] * N
    used = [False] * M

    visit = [False] * N

    

    res = []
    for start in range(N):
        if visit[start]:
            continue

        visit[start] = True
        st = [start]
        tmp_res = []
        while st:
            v = st[-1]
            visit[v] = True

            if done_counter[v] == len(edge[v]):
                tmp_res.append(v)
                st.pop()
                continue

            nv,eid = edge[v][done_counter[v]]
            done_counter[v] += 1
            if used[eid]:
                continue
            used[eid] = True
            st.append(nv)
        res.append(tmp_res)
    
    
    return res

def solve(N):
    if N == 1:
        return ("No",[])
    
    if N & 1 == 0:
        edge = [[] for v in range(1<<N)]
        nxt_eid = 0
        for v in range(1<<N):
            for k in range(N):
                if v>>k & 1 == 0:
                    edge[v].append((v^(1<<k),nxt_eid))
                    edge[v^(1<<k)].append((v,nxt_eid))
                    nxt_eid += 1
        
        euler_tour = find_euler_tour(1<<N,nxt_eid,edge)
        res = [[0]*(1<<N) for i in range(1<<N)]

        for ce in euler_tour:
            for u,v in zip(ce,ce[1:]):
                res[u][v] = 1
                res[v][u] = -1
        
        return ("Yes",res)
    
    else:
        edge = [[] for v in range(1<<N)]
        nxt_eid = 0
        for v in range(1<<N):
            for k in range(3,N):
                if v>>k & 1 == 0:
                    edge[v].append((v^(1<<k),nxt_eid))
                    edge[v^(1<<k)].append((v,nxt_eid))
                    nxt_eid += 1
        
        euler_tour = find_euler_tour(1<<N,nxt_eid,edge)
        res = [[0]*(1<<N) for i in range(1<<N)]

        for ce in euler_tour:
            for u,v in zip(ce,ce[1:]):
                res[u][v] = 1
                res[v][u] = -1

        for root in range(0,1<<N,8):
            two_in = [0,3,5,6]
            two_out = [1,2,4,7]

            for a,b in zip(two_in,two_out):
                res[a^root][b^root] = -2
                res[b^root][a^root] = 2
            
            for a in two_in:
                for k in [2,4]:
                    res[a^root][a^k^root] = 1
                    res[a^k^root][a^root] = -1
        
        return ("Yes",res)
    
def checker(N,res):
    for i in range(1<<N):
        for k in range(N):
            if res[i][i^(1<<k)] == 0:
                print(i,i^(1<<k))
                return False
        for j in range(1<<N):
            d = i^j
            if (d & (d-1)) != 0 and res[i][j]!=0:
                return (False,"must 0")
            
        if sum(res[i])!=0:
            return (False,"sum")
    return True

N = int(input())
ans,res = solve(N)

if ans == "Yes":
    print(ans)
    for i in range(1<<N):
        print(*res[i])
else:
    print(ans)