import sys
input = sys.stdin.readline
N = int(input())

res = [[0] * N for _ in range(2 * N - 1)]
def solve(n, l):
  global res
  if n == 1:
    res[0][0] = l
    return
  for i in range(2 * n): res[2 * n - 2 - (i & 1)][i >> 1] = l + i
  for i in range(2 * n - 1): res[i][n - 1] = l + 2 * n - 2
  #print(res)
  solve(n - 1, l + 1)

solve(N, 1)
print(2 * N - 1)
for r in res: print(*r)