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

data = [0 for i in range(2*10**5+1)]
for a in A:
    data[a]+=1

cum = [data[i] for i in range(2*10**5+1)]
cum_val = [i*data[i] for i in range(2*10**5+1)]
for i in range(1,2*10**5+1):
    cum[i] += cum[i-1]
    cum_val[i] += cum_val[i-1]

ans = 0
Max = 2*10**5
for i in range(1,Max+1):
    ans += data[i] * cum_val[i-1]
    for j in range(1,Max//i + 1):
        L = i*j
        R = min(i*(j+1) - 1,Max)
        ans -= i * j * (cum[R]-cum[L-1]) * data[i]
    ans += data[i] * (cum_val[Max] - cum_val[i-1])

print(ans)