Pagini recente » Cod sursa (job #2225095) | Cod sursa (job #1826444) | Cod sursa (job #2108123) | Cod sursa (job #2038181) | Cod sursa (job #162216)
Cod sursa(job #162216)
#include<stdio.h>
#include<math.h>
int n,max,prim;
int a[100002],i,j,ok,p,aux,nr,fact,cs;
//const int sirprim[20] = {0,2,3,5,7,11,13,17,19};
int v[1000002],contor;
int frcv[1000002];
int main()
{
freopen("pairs.in","r",stdin);
freopen("pairs.out","w",stdout);
scanf("%ld",&n);
for(i=1;i<=n;i++)
{
scanf("%ld",&a[i]);
frcv[a[i]]=1;
if (a[i]>max)
max=a[i];
}
for(i=2;i<=max;i++)
{
for(j=i;j<=max;j+=i)
if (frcv[j]) v[i]++;
}
for(p=2;p<=max;p++)
{
aux=p;
nr=0;
fact=1;
while(aux!=1)
{
fact++;
contor=0;
while (aux%fact==0) {aux/=fact;contor++;}
if (contor==1) nr++;
else
if (contor>1) {break;nr=-1;}
if (fact>sqrt(aux)&&aux!=1) {nr++;break;}
}
if (nr!=-1)
if (nr%2) cs+=(v[p]*(v[p]-1))/2;
else cs-=(v[p]*(v[p]-1))/2;
}
printf("%ld",n*(n-1)/2-cs);
return 0;
}