#include <bits/stdc++.h>
  
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/pb_ds/detail/standard_policies.hpp>
using namespace __gnu_pbds;
 
using namespace std;
 
#if !ONLINE_JUDGE
#define debug
#endif
 
typedef long long ll;
typedef vector<int> vi;
typedef pair<int,int> ii;
typedef vector<pair<int,int> > vii;
typedef pair<ll,ll> pll;
typedef vector<vector<int> > vvi;
typedef vector<vector<ii> > vvii;
 
#define pb push_back
#define mp make_pair
#define all(x) x.begin(),x.end()
#define sz(x) (int)x.size()
#define fill(a,x) memset(a,x,sizeof(a))
#define ff first
#define ss second
#define trav(a, x) for(auto& a : x) 
#define FOR(i,a,b) for(int i = a; i<=b; ++i)
#define NFOR(i,a,b) for(int i = a; i>=b; --i)

#ifdef debug
#define dbg(...) { cerr<<"[ "; dump(#__VA_ARGS__, __VA_ARGS__);}
#undef endl
template <typename Arg1,typename Arg2>
ostream& operator << (ostream& out, const pair<Arg1,Arg2> &x) {
    return out<<"("<<x.ff<<","<<x.ss<<")";
}
 
template <typename Arg1>
ostream& operator << (ostream& out, const vector<Arg1> &a) {
    out<<"[";for(const auto &x:a)out<<x<<",";return out<<"]";
}
 
template <typename Arg1>
ostream& operator << (ostream& out, const set<Arg1> &a) {
    out<<"[";for(const auto &x:a)out<<x<<",";return out<<"]";
}
 
template <typename Arg1,typename Arg2>
ostream& operator << (ostream& out, const map<Arg1,Arg2> &a) {
    out<<"[";for(const auto &x:a)out<<x<<",";return out<<"]";
}
 
template <typename Arg1,typename Arg2>
ostream& operator << (ostream& out, const unordered_map<Arg1,Arg2> &a) {
    out<<"[";for(const auto &x:a)out<<x<<",";return out<<"]";
} 
 
template <typename Arg1>
void dump(const string name, Arg1&& arg1){
    cerr << name << " : " << arg1 << " ] " << endl;
}
 
template <typename Arg1, typename... Args>
void dump(const string names, Arg1&& arg1, Args&&... args){
    const string name = names.substr(0,names.find(','));
    cerr<<name<<" : "<<arg1<<" | ";
    dump(names.substr(1+(int)name.size()), args...);
}
#else
    #define dbg(args...)
#endif
 
ll gcd(ll a,ll b){if (a==0) return b;return gcd(b%a,a);}
ll lcm(ll a,ll b){return a*(b/gcd(a,b));}
ll Abs(ll a){if(a>0)return a;return -a;}
ll Ceil(ll a,ll b){if(a%b==0)return a/b;else return a/b+1;}
double Abs(double a){if(a>0)return a;return -a;}
 
//*X.find_by_order(2) element at index=2
//X.order_of_key(1) how many elements strictly less than 1
 
#define ordered_set tree<int, null_type,less_equal<int>, rb_tree_tag,tree_order_statistics_node_update> 

const ll INF = 1e18;
const int mod = 1e9 + 7;
 
inline int pow_(ll x, ll y,ll p){
  int r = 1; while(y){if(y & 1) r = r * x % p; y >>= 1; x = x * x % p;} return r;
}
inline int inv_(int x){return pow_(x, mod - 2, mod);}
inline int add(ll a,ll b){a += b; if(a >= mod) a -= mod; return a;}
inline int mul(int a, int b){return a * 1LL * b % mod;}
inline int sub(int a, int b){a -= b; if(a < 0) a += mod; return a;}
  
mt19937_64 mrand(chrono::steady_clock::now().time_since_epoch().count());
 
const int N = 1e4 + 5;

int arr[N];
int arr1[N];

const int m = 1e6 + 5;

int p[m] = {0};

map<int,int> a[N];
 
int main(int argc, char const *argv[])
{
    #ifdef LOCAL_TEST
    freopen("in.txt","r",stdin);
    //freopen("out.txt", "w", stdout);
    #else
    ios_base::sync_with_stdio(false);cin.tie(NULL);
    #endif

    auto start = std::chrono::high_resolution_clock::now();

    int n;
    cin>>n;

    for(int i=1;i<=n;i++){
        cin>>arr[i];
        arr1[i] = arr[i];
    }

    for(int i=1;i<=n;i++){
        vector<ii> v;
        for(ll j=2;j*j<=arr[i];j++){
            int d = 0;
            while(arr[i]%j == 0 && arr[i] != 1){
                arr[i] /= j;
                d++;
            }

            if(d){
                a[i].insert({j,d});
                v.pb({j,d});
            }
        }

        if(arr[i] != 1){
            a[i].insert({arr[i],1});
            v.pb({arr[i],1});
        }

        for(int i=0;i<v.size();i++){
            if(v[i].ss > p[v[i].ff]){
                p[v[i].ff] = v[i].ss;
            }
        }
    }

    ll ans_ = 0;

    vector<ii> l;

    for(int i=1;i<=1e6;i++){
        if(p[i]){
            // dbg(i,p[i]);
            l.pb({i,p[i]});
        }
    }

    for(int i=1;i<=n;i++){
        ll ans = 1;
        ll ans1 = 1;
        for(int j=0;j<l.size();j++){
            int val = l[j].ff;
            int val2 = l[j].ss;
            if(a[i].find(val) != a[i].end()){
                int val1 = a[i][val];
                // dbg(val,val1,val2,i);
                if(val1 < val2){
                    ans = mul(ans,pow_(val,val2-val1,mod));
                    ans1 = mul(ans1,pow_(val,val1,mod));
                } else {
                    ans1 = mul(ans1,pow_(val,val2,mod));
                }
            } else {
                ans = mul(ans,pow_(val,val2,mod));
            }
        }


        ans1 = mul(arr1[i],inv_(ans1));


        ans = mul(ans,inv_(ans1));


        ans_ = add(ans_,ans);
    }

    cout << ans_ << "\n";

    auto finish = std::chrono::high_resolution_clock::now();  
    std::chrono::duration<double> elapsed = finish - start;
    cerr << "Elapsed time: " << elapsed.count() << " s\n";
    
    return 0;    
}

Measured memory_bytes: 5632