n = int(input())
*a, = map(int,input().split())
N = 2*10**5+3
r = [0]*N
for ai in a:
    r[ai] += 1
acc = r[:] + [0]
for i in range(N)[::-1]:
    acc[i] += acc[i+1]

ans = sum(a)*(n-r[1])
for i in range(2,N):
    v = 0
    for j in range(1,N//i+1):
        v += j*(acc[j*i] - acc[min(N,(j+1)*i)])
    ans -= r[i]*i*v
print(ans)