#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+5,mod=1e9+7;
int n,a[N],vis[N];
vector<int>v;
bool vv[N];
ll qpow(ll a,ll n)
{
    ll ans=1;
    while(n)
    {
        if(n&1) ans=ans*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return ans;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)
    {
        int x=a[i];
        for(int j=2;j*j<=x;j++)
        {
            int tot=0;
            while(x%j==0)
            {
                x/=j;tot++;
            }
            if(tot)
            {
                if(!vv[j]) vv[j]=true,v.push_back(j);
            }
            vis[j]=max(vis[j],tot);
        }
        if(x>1)
        {
            vis[x]=max(vis[x],1);
            if(!vv[x]) vv[x]=true,v.push_back(x);
        }
    }
    ll ans=0;
    for(int i=1;i<=n;i++)
    {
        ll res=1;
        for(int j=0;j<v.size();j++)
        {
            int sum=vis[v[j]];
            while(a[i]%v[j]==0)
            {
                a[i]/=v[j];
                sum--;
            }
            res=res*qpow(v[j],sum)%mod;
        }
        ans=(ans+res)%mod;
    }
    printf("%lld\n",ans);
}
