#include <fstream>
#include <vector>
#define N 100001
#define M 1000001
using namespace std;

ifstream fin("pairs.in");
ofstream fout("pairs.out");

bool verif[M];
int prime[N], pi, cont[M];
vector <int> divi;

void ciur()
{
    long long i, j;
    for(i=2;i<M;i++)
    {
        if(!verif[i])
        {
            prime[++pi]=i;
            for(j=i*i;j<M;j+=i)
            {
                verif[j]=true;
            }
        }
    }
}

int main()
{
    int n, x, i, j, ind, sz, semn, nr;
    long long sol;
    fin>>n;
    sol=1LL*n*(n-1)/2;
    ciur();
    for(ind=0;ind<n;ind++)
    {
        fin>>x;
        divi.clear();
        for(i=1;x>1&&prime[i]*prime[i]<=x;i++)
        {
            if(x%prime[i]==0)
            {
                divi.push_back(prime[i]);
                while(x%prime[i]==0) x/=prime[i];
            }
        }
        if(x>1) divi.push_back(x);
        sz=divi.size();
        for(i=1;i<(1<<sz);i++)
        {
            nr=semn=1;
            for(j=0;j<sz;j++)
            {
                if(i&(1<<j))
                {
                    semn=-semn;
                    nr*=divi[j];
                }
            }
            sol+=cont[nr]*semn;
            cont[nr]++;
        }
    }
    fout<<sol;
}