N = int(input())
A = list(map(int,input().split()))
M = 200001
D = [0] * M
DS = [0] * M
S = [0] * M
for a in A:
    D[a] += 1

for i in range(1,M):
    S[i] = S[i-1] + i * D[i]
    DS[i] = DS[i-1]+D[i]

ans = 0
for i in range(1,M):
    if D[i] > 0:
        ans+= S[i-1] * D[i]
        for j in range(1,M//i):
            ans += D[i] * (S[i*(j+1)-1]-S[i*j] - i * j * (DS[i*(j+1)-1] - DS[i*j]))

print(ans)