# 最小公倍数(ユーグリッド互除法)
def gcd(a,b):
    if b == 0:
        return a
    else:
        return gcd( b, a % b)

K = int(input()) +1

# データパターンを求める式
ans = 0
for h in range(1, K):
  for i in range(1, K):
    for j in range(1, K):
      # 最小公倍数算出
      ans += gcd(h,gcd(i,j))

print(ans)