#include <bits/stdc++.h>
#define DIM 1000005
using namespace std;
int pr[DIM];
bool ap[DIM];
long long x,n,ans,nr,i,j,ma;
int main()
{
    ifstream f ("pairs.in");
    ofstream g ("pairs.out");
    f>>n;
    for(i=1; i<=n; ++i)
    {
        f>>x;
        ap[x]=1;
        ma=max(ma,x);
    }
    for(i=2; i<=ma; ++i)
    {
        if(!pr[i])
            for(j=i; j<=ma; j+=i)
                pr[j]++;
        if(1LL*i*i>ma)continue;
        for(j=i*i; j<=ma; j+=i*i)
            pr[j]=-1;
    }
    for(i=2; i<=ma; ++i)
        if(pr[i]>-1)
        {
            nr=0;
            for(j=i; j<=ma; j+=i)
                nr+=ap[j];
            if(pr[i]%2)ans+=(nr)*(nr-1);
            else ans-=nr*(nr-1);
        }
    g<<n*(n-1)/2-ans/2;
    return 0;
}