#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

#pragma GCC target ("avx2")
#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")
#pragma comment(linker, "/stack:200000000")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

using namespace std;
using namespace __gnu_pbds;

typedef long long ll;
typedef string str;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector <pii> vpii;
typedef vector <pll> vpll;
typedef map <str,int> mapsi;
typedef map <str,int> :: iterator mapsitr;
typedef map <int , int> mint;
typedef map <ll , ll> mll;
typedef set <int> si;
typedef set <ll> sll;
typedef si :: iterator sitr;
typedef si :: reverse_iterator rsitr;
typedef sll :: iterator sltr;
typedef sll :: reverse_iterator rsltr;
#define mset multiset

typedef mset <int> msi;
typedef mset <ll> msll;
typedef msi :: iterator msitr;
typedef msi :: reverse_iterator msritr;
typedef msll :: iterator msltr;
typedef msll :: reverse_iterator mslritr;

#define ordered_set tree<int, null_type,less<int>, rb_tree_tag,tree_order_statistics_node_update>
#define mp make_pair
#define pb push_back
#define pob pop_back
#define pf push_front
#define pof pop_front
#define fi first
#define se second
#define fs first.second
#define ss second.second
#define ff first.first
#define sf second.first
#define newl '\n'
#define fbo find_by_order
#define ook order_of_key

char to_upper (char x){
    if( 97 <= int(x) && int(x) <= 122)return char(x-32);
    else if( 65 <= int(x) && int(x) <= 90)return x;
}
char to_lower (char x){
    if( 97 <= int(x) && int(x) <= 122)return x;
    else if( 65 <= int(x) && int(x) <= 90)return char(x+32);
}
int numerize (char x){
    if(48 <= int(x) && int(x) <= 57)return int(x-'0');
    else if( 97 <= int(x) && int(x) <= 122)return int(x-96);
    else if( 65 <= int(x) && int(x) <= 90)return int(x-64);
}
bool isect (int l1, int r1, int l2, int r2){
    pii p1,p2;
    p1 = mp(l1,r1); p2 = mp(l2,r2);
    if(p1>p2)swap(p1,p2);
    if(p2.fi <= p1.se)return true;
    else return false;
}
ll quickpow (ll num1, ll num2, ll MOD){
    if(num2==0)return 1;
    else if(num2==1)return num1;
    else{
        ll temp = quickpow (num1,num2/2,MOD); ll res = ((temp%MOD) * (temp%MOD))%MOD;
        if(num2%2==1) res = ((res%MOD)*(num1%MOD))%MOD; return res;
    }
}
ll invmod (ll num, ll MOD){return quickpow (num,MOD-2,MOD);}
ll gcd (ll num1, ll num2){
    if(num1 < num2) swap(num1,num2); ll num3 = num1 % num2 ;
    while(num3 > 0){ num1 = num2; num2 = num3; num3 = num1 % num2;}
    return num2;
}
ll lcm (ll num1 , ll num2){return (ll) (num1/__gcd(num1,num2))*num2;}
// end of Template

vi prime;
bool isprime[1000010];
ll mx[1000010],a[10010],ans;
ll md = 1e9+7;
int n;
void sieve(int x){
    for(int i=2;i<=x;i++) isprime[i] = true;
    for(int i=2;i<=x;i++){
        if(isprime[i]){
            prime.pb(i);
            for(int j=i;j<=x;j+=i) isprime[j] = false;
        }
    }
}
vpii factor[10010];
vpii L;
void factorize(int idx){
    map <int,int> f;
    int x = a[idx];
    while(x>1){
        for(int i=0;i<prime.size();i++){
            int temp = prime[i];
            if(temp > sqrt(x)){
                f[x] ++;
                x = 1;
                break;
            }

            if(x % temp == 0){
                f[temp] ++;
                x/=temp;
                break;
            }
        }
    }

    for(map<int,int> :: iterator itr = f.begin(); itr!=f.end(); ++itr){
        factor[idx].pb(mp(itr->fi,itr->se));
        mx[itr->fi] = max(mx[itr->fi],(ll)itr->se);
    }
}
int main(){

//    freopen("input.txt", "r", stdin);
//    freopen("output.txt", "w", stdout);

    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    cout.tie(NULL);

    sieve(1e6);

    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        factorize(i);
//        cout<<i<<" : ";
//        for(int j=0;j<factor[i].size();j++) cout<<factor[i][j].fi<<"^"<<factor[i][j].se<<" ";cout<<newl;
    }

    for(int i=1;i<=1e6;i++){
        if(mx[i]){
            L.pb(mp(i,mx[i]));
//            cout<<i<<" "<<mx[i]<<newl;
        }
    }

    for(int idx=1;idx<=n;idx++){
        ll temp = 1;
        for(int i=0;i<L.size();i++){
            int a = L[i].fi;
            int b = L[i].se;
            int x = 0;
            for(int j=0;j<factor[idx].size();j++){
                if(factor[idx][j].fi == a){
                    x = factor[idx][j].se; break;
                }
            }
            temp = (temp * quickpow(a,b-x,md)) % md;
        }
        ans = (ans + temp) % md;
    }
    cout<<ans<<newl;

}


Measured memory_bytes: 9976