#include <iostream>
#include <vector>
#include <algorithm>

using namespace std;

long long mod = 1000000007;

int gcd (int a,int b){
    return b ? gcd(b,a%b):a;
}

int main(){
    int N;
    cin>>N;
    int tmp;
    vector<int>G;
    vector<int>M;
    int res=0;
    for(int i=0;i<N;++i){
        cin>>tmp;
        G.emplace_back(tmp);
        M.emplace_back(tmp);
        res+=G[i]*M[i];
        res %= mod;
    }
    for(int i=N-1;i>=0;i--){
        for(int j=0;j<i;j++){
            G[j]=gcd(G[j],G[j+1]);
            M[j]=max(M[j],M[j+1]);
            res+= G[j]*M[j];
            res %= mod;
        }
    }


    cout<<res<<endl;

    return 0;
}