N = int(input())
A = list(range(1,N+1))
ans = sum(A)*N
A.sort()
from bisect import *
for a in A:
  for p in range(a,A[-1]+1,a):
    ind = bisect_left(A,p)
    ans -= a*(N-ind)
print(ans)