import core.bitop, std.algorithm, std.ascii, std.bigint, std.conv, std.math,
    std.functional, std.numeric, std.range, std.stdio, std.string, std.random,
    std.typecons, std.container, std.format;
static import std.ascii;

// dfmt off
T lread(T = long)(){return readln.chomp.to!T();}
T[] aryread(T = long)(){return readln.split.to!(T[])();}
void scan(TList...)(ref TList Args){auto line = readln.split();
foreach (i, T; TList){T val = line[i].to!(T);Args[i] = val;}}
alias sread = () => readln.chomp();enum MOD = 10 ^^ 9 + 7;
alias PQueue(T, alias less = "a<b") = BinaryHeap!(Array!T, less);
// dfmt on

void main()
{
    long N = lread();
    auto A = aryread();
    // auto B = A.map!(x => BigInt(x)).array;
    auto l = lcm(A);
    long ans;
    foreach (a; A)
    {
        long x = (l * powmod(a, MOD - 2, MOD)) % MOD;
        ans = (ans + x) % MOD;
    }
    writeln(ans);
}

long lcm(long[] A)
{
    long[long] f;
    foreach (a; A)
    {
        auto d = factorize(a);
        foreach (key, value; d)
        {
            f[key] = max(f.get(key, 0), value);
        }
    }
    // writeln(f);
    long ans = 1;
    foreach (key, value; f)
    {
        ans = powmod(key, value, MOD) * ans % MOD;
    }
    return ans;
}

/// x^^n % m
T powmod(T = long)(T x, T n, T m)
{
    if (n < 1)
        return 1;
    if (n & 1)
    {
        return x * powmod(x, n - 1, m) % m;
    }
    T tmp = powmod(x, n / 2, m);
    return tmp * tmp % m;
}

/// 素因数分解
long[long] factorize(long x)
{
    assert(0 < x, "x is negative");
    long[long] ps;
    while ((x & 1) == 0)
    {
        x /= 2;
        ps[2] = (2 in ps) ? ps[2] + 1 : 1;
    }
    for (long i = 3; i * i <= x; i += 2)
        while (x % i == 0)
        {
            x /= i;
            ps[i] = (i in ps) ? ps[i] + 1 : 1;
        }
    if (x != 1)
        ps[x] = (x in ps) ? ps[x] + 1 : 1;
    return ps;
}
