#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=1e6+5;
int a[maxn];
vector<int> t[maxn];//质数出现的位置
signed main(){
    //保证只被最小质因子筛去
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    vector<bool> vis(1000000+1);
    vector<int> prm;//1e6以内的质数
    for(int i=2;i<=1000000;i++){
        if(!vis[i])prm.push_back(i);
        for(int j=0;j<prm.size();j++){
            if(1ll*prm[j]*i>1000000)break;
            vis[i*prm[j]]=true;
            if(i%prm[j]==0)break;
        }
    }
    for(int i=1;i<=n;i++){
        int x=a[i];
        for(int j=0;prm[j]*prm[j]<=a[i];j++){
            if(x%prm[j]==0){
                t[prm[j]].push_back(i);
                while(x%prm[j]==0)x/=prm[j];
            }
        }
        if(x>1)t[x].push_back(i);
    }
    // for(int i=0;i<30;i++){
    //     for(int j=0;j<t[prm[i]].size();j++)
    //         cout<<prm[i]<<":"<<t[prm[i]][j]<<":"<<a[t[prm[i]][j]]<<":"<<n-t[prm[i]][j]+1<<endl;
    //     cout<<"--------------"<<endl;
    // }
    int ans=0;
    for(int i=0;i<prm.size();i++){
        int x=prm[i];
        for(int j=0;j<t[x].size();j++){
            int l=t[x][j];
            if(j!=0)l-=t[x][j-1];
            ans+=l*(n-t[x][j]+1);
        }
    }
    cout<<ans<<endl;
    return 0;
}