N = int(input())
A = list(map(int, input().split()))
ans = sum(A)*N
M = max(A)
lis = [0]*(M+1)
cum = [0]*(M+1)
for a in A:
  lis[a] += 1
for i in range(1,M+1):
  cum[i] = cum[i-1]+lis[i]

for i in range(1,M+1):
  if lis[i]==0:
    continue
  for p in range(i,M+1,i):
    ans -= lis[i]*i*(N-cum[p-1])
print(ans)