N = int(input())
A = list(map(int, input().split()))
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)