#include <fstream>
#include <cstring>
#include <cmath>
#define maxn 1000000
#define MAX 1000100
using namespace std;
int x[MAX],rad,prime[300002],b[300002],m,divizori[1000002],multime[300002],a[100002],n;
long long prod,sol;
bool neprime[MAX],in[MAX];

void ciur()
{
	int pr=0;
	for(int i=2;i<=1000000;i++)
	if(!divizori[i])
	{
		prime[++pr]=i;
		divizori[i]=1;
		for(int j=i+i;j<=1000000;j+=i) divizori[j]++;
	}
}
int main()
{
	int i,j;
	ifstream fi("pairs.in");
	ofstream fo("pairs.out");
	fi>>n;
	sol=n*(n-1)/2;
	ciur();
	memset(neprime,0,sizeof(neprime));
	//neprime[]=in multime
	for(i=1;i<=n;i++) 
	{
		fi>>a[i];
		in[a[i]]=1;
	}
	for(i=2;i<=maxn;i++)
	for(j=i;j<=maxn;j+=i) 
		if(in[j]) x[i]++;

	for(i=2;i<=maxn;i++)
	if(divizori[i] and x[i]>1)
	{
		if(divizori[i]%2) sol-=x[i]*(x[i]-1)/2; else sol+=x[i]*(x[i]-1)/2;
	}
	fo<<sol<<"\n";
	return 0;
}