#include<bits/stdc++.h>
using namespace std;using ll=long long;using vi=vector<int>;using vvi=vector<vi>;using vl=vector<ll>;using vvl=vector<vl>;using P=pair<int,int>;using PL=pair<ll,ll>;using vp=vector<P>;using vpl=vector<PL>;
#define pqueue priority_queue
template<typename T>constexpr auto inf=numeric_limits<T>::max()/2;constexpr int INF=inf<int>,MOD=1000000007;constexpr ll LINF=inf<ll>;
#define _ol3(_1,_2,_3,name,...)name
#define _rep(i,n)_repi(i,0,n)
#define _repi(i,a,b)for(int i=a,i##_l=(b);i<i##_l;++i)
#define REP(...)_ol3(__VA_ARGS__,_repi,_rep,)(__VA_ARGS__)
#define REPR(i,n)for(int i=n-1;i>=0;--i)
#define REPA(i,v)REP(i,(v).size())
#define all(v)(v).begin(),(v).end()
#define rall(v)(v).rbegin(),(v).rend()
#define bit(n)(1ll<<(n))
#define F first
#define S second
#define endl '\n'
#define cho(n,a,b)cout<<((n)?a:b)<<endl
void YES(int n){cho(n,"YES","NO");}void Yes(int n){cho(n,"Yes","No");}void Poss(int n){cho(n,"Possible","Impossible");}
void _print(ostream&){}template<class T,class...U>void _print(ostream&s,const T&t,const U&...u){s<<t<<(sizeof...(u)?' ':'\n');_print(s,u...);}
template<class...T>void print(const T&...t){_print(cout,t...);}template<class...T>void dprint(const T&...t){_print(cerr,t...);}
#ifndef LOCAL
struct osd{template<class T>osd&operator<<(const T&t){return*this;}};osd cer_;
#define dprint(...)
#define cerr cer_
#endif
#define dbg(x) dprint("@l",__LINE__,':',#x,'=',x)
template<class T>bool chmax(T&a,const T&b){bool x=a<b;x?a=b:b;return x;}template<class T>bool chmin(T&a,const T&b){bool x=a>b;x?a=b:b;return x;}
template<class T,class U>ostream&operator<<(ostream&o,const pair<T,U>&p){return o<<p.F<<' '<<p.S;}template<class T,class U>istream&operator>>(istream&i,pair<T,U>&p){return i>>p.F>>p.S;}
template<class T>class iterable{static false_type c(string v);template<class U>static auto c(U v)->decltype(all(v),true_type());static false_type c(...);public:const static bool value=decltype(c(declval<T>()))::value;};
template<class T,enable_if_t<iterable<T>::value,int> =0>ostream&operator<<(ostream&o,const T&v){for(auto&&i:v)o<<i<<' ';return o;}
template<class T>istream&operator>>(istream&i,vector<T>&v){for(T&j:v)i>>j;return i;}template<class T>vector<T>&operator<<(vector<T>&v,const T&t){v.push_back(t);return v;}
template<class T,class U>queue<T,U>&operator<<(queue<T,U>&v, const T&t){v.push(t);return v;}template<class T,class U>queue<T,U>&operator>>(queue<T,U>&v, T&t){t=v.front();v.pop();return v;}
template<class T,class U,class V>pqueue<T,U,V>&operator<<(pqueue<T,U,V>&v, const T&t){v.push(t);return v;}template<class T,class U,class V>pqueue<T,U,V>&operator>>(pqueue<T,U,V>&v, T&t){t=v.front();v.pop();return v;}
template<class T,class U>stack<T,U>&operator<<(stack<T,U>&v, const T&t){v.push(t);return v;}template<class T,class U>stack<T,U>&operator>>(stack<T,U>&v, T&t){t=v.front();v.pop();return v;}
template<class T,class U,class V>set<T,U,V>&operator<<(set<T,U,V>&v,const T&t){v.insert(t);return v;}template<class T,class U,class V>multiset<T,U,V>&operator<<(multiset<T,U,V>&v,const T&t){v.insert(t);return v;}
template<class T,enable_if_t<is_integral<T>::value,int> =0>T mod(T a, T b){if(a>0)return a%b;return b-(-a%b);}template<class T,enable_if_t<is_floating_point<T>::value,int> =0>T mod(T a, T b){if(a>0)return fmod(a,b);return b-fmod(-a,b);}
template<class T,class U>pair<T,U>&operator+=(pair<T,U>&a, pair<T,U>&b){a.F+=b.F;a.S+=b.S;return a;}template<class T,class U>pair<T,U>operator+(const pair<T,U>&a, const pair<T,U>&b){T r=a;return r+=b;}
template<class T,class U>pair<T,U>&operator-=(pair<T,U>&a, pair<T,U>&b){a.F-=b.F;a.S-=b.S;return a;}template<class T,class U>pair<T,U>operator-(const pair<T,U>&a, const pair<T,U>&b){T r=a;return r-=b;}
template<class T,class U>pair<T,U>operator-(const pair<T,U>&a){return mp(-a.F, -a.S);}
template<class T,class U>double std::norm(const pair<T,U>&p){return norm(p.F)+norm(p.S);}template<class T,class U>double std::abs(const pair<T,U>&p){return sqrt(norm(p));}
inline ll square(ll a){a%=MOD;return (a*a)%MOD;}
ll lpow(ll a,ll b){
    if(b<0)assert(0);
    a%=MOD;
    b%=MOD-1;
    if(b==0)return 1;
    if(b==1)return a;
    if(b%2){
        return (square(lpow(a,b/2)%MOD)*a)%MOD;
    }else{
        return square(lpow(a,b/2)%MOD);
    }
}

ll inv(ll a){
    return lpow(a,MOD-2);
}

ll fact(ll n){
    ll ans=1;
    while(n){
        ans=(ans*n--)%MOD;
    }
    return ans;
}

ll factto(ll a,ll b){
    ll ans=1;
    for(ll i=a;i<=b;i++){
        ans=(ans*i)%MOD;
    }
    return ans;
}

ll perm(ll n,ll r){ // nPr
    return factto(n-r+1,n)%MOD;
}

ll comb(ll n,ll r){ // nCr
    return (perm(n,r)*lpow(fact(r),MOD-2))%MOD;
}

ll homo(ll n,ll r){ // nHr
    return comb(n+r-1,r);
}
class yurisu{
public:
    map<ll,int> f;
    yurisu(int n){
        if(n == 1) { // n=1 の素因数分解は n^1
            //f[n] = 1;
            return ;
        }
        for(int i = 2, _n = n; i*i <= _n; ++i) {
            while(n % i == 0) {
                ++f[i]; // 素数i^{f[i]}
                n /= i;
            }
        }
        if(n != 1) f[n] = 1;
    }
    void lcm(const yurisu&o){
        for(const auto& x:o.f)if(f.find(x.F)==f.end())f[x.F]=x.S;
        for(auto& x:f)if(o.f.find(x.F)!=o.f.end())chmax(f[x.F],o.f.at(x.F));
    }
    ll div(const yurisu&o){
        ll ret=1;
        for(auto&x:f){
            if(o.f.find(x.F)==o.f.end())
                (ret*=lpow(x.F, x.S))%=MOD;
            else
                (ret*=lpow(x.F, x.S-o.f.at(x.F)))%=MOD;
        }
        return ret;
    }
    void show(){dprint(f);}
    ll get(){
        ll ret=1;
        for(const auto&x:f){
            ret*=lpow(x.F,x.S);
            ret%=MOD;
        }
        return ret;
    }
};

int main(){
    cin.tie(0);ios::sync_with_stdio(0);
    cout<<setprecision(10);
    int n;
    cin>>n;
    vi a(n);
    cin>>a;
    vector<yurisu> b;
    REP(i,n)b.emplace_back(a[i]);
    yurisu l(1);
    //yurisu one(1);
    REP(i,n)l.lcm(b[i]);
    
    ll lc=l.get();
    
    ll ans=0;
    REP(i,n){
        dprint(l.div(b[i]));
        //ans+=l.div(b[i]);
        ans+=lc*inv(b[i].get());
        ans%=MOD;
    }
    cout<<ans<<endl;
}


Measured memory_bytes: 3896