import math

N = int(input())
A = [int(x) for x in input().split()]

sum = int(0)

def P(x, max, dev):
    if A[x] > max:
        max = A[x]
    tmp = math.gcd(dev, A[x])
    if tmp < dev:
        dev = tmp
    return max, dev

for i in range(0, N):
    max = int(0)
    dev = A[i]
    for j in range(i, N):
        max, dev = P(j, max, dev)
        sum += max*dev

sum %= 1000000007
print(int(sum))