N = int(input())

A = list(map(int, input().split()))

A.sort()

ans = sum(A) * N

from bisect import bisect_left
X = A[-1]
note = [0] * (10 ** 5 + 1)

for a in A:
    note[a] += 1
    
S = note.copy()

for i in range(10 ** 5, 0, -1):
    S[i - 1] += S[i]

for i in range(1, 10 ** 5 + 1):
    for j in range(i, 10 ** 5 + 1, i):
        ans -= note[i] * i * S[j]
    
print(ans)