from atcoder.fenwicktree import FenwickTree

MOD = 10 ** 9 + 7
N = int(input())
A = list(map(int,input().split()))
A.sort()
Amax = A[-1]
ft = FenwickTree(Amax + 1)
t1 = 0
t2 = 0
for i in range(N):
    t1 += A[i] * (2 * i - N + 1)
for i in reversed(range(N)):
    a = A[i]
    tmp = 0
    for j in range(a,Amax+1,a):
        tmp += j * ft.sum(j,min(j+a,Amax+1))
        tmp %= MOD
    ft.add(a,1)
    t2 += tmp
    t2 %= MOD
print((t2-t1)%MOD)