# include <fstream>
# define DIM 1000010
using namespace std;
ifstream fin("pairs.in");
ofstream fout("pairs.out");
int d[DIM],f[DIM],nr[DIM],n,x,i,j,val,maxim;
long long sol;
int main () {
    fin>>n;
    for(i=1;i<=n;i++){
        fin>>x;
        f[x]=1;
        maxim=max(maxim,x);
    }
    for(i=1;i<=maxim;i++)
            d[i]=1;
    for(i=2;i<=maxim;i++)
        if(d[i]==1)
            for(j=i;j<=maxim;j+=i){
                d[j]*=i;
                nr[j]++;
            }
    for(i=1;i<=maxim;i++){
        if(d[i]==i){
            val=0;
            for(j=i;j<=maxim;j+=i)
                if(f[j])
                    val++;
            if(nr[i]%2==0)
                sol+=(1LL*val*(val-1))/2;
            else
                sol-=(1LL*val*(val-1))/2;
        }
    }
    fout<<sol<<"\n";
    return 0;
}