Pagini recente » Cod sursa (job #1793866) | Cod sursa (job #662159) | Cod sursa (job #457797) | Cod sursa (job #1424236) | Cod sursa (job #2295386)
# include <fstream>
# define DIM 1000010
using namespace std;
ifstream fin("pairs.in");
ofstream fout("pairs.out");
int d[DIM],f[DIM],nr[DIM],p[DIM],v[DIM/10],n,x,i,j,val,maxim;
long long sol;
int main () {
fin>>n;
for(i=1;i<=n;i++){
fin>>x;
f[x]=1;
v[i]=x;
maxim=max(maxim,x);
}
/*for(i=2;i*i<=maxim;i++)
if(p[i]==0){
val=0;
for(j=i;j<=maxim;j+=i){
p[j]=1;
if(f[j])
val++;
}
if(val>0)
p[i]=0;
}*/
for(i=1;i<=maxim;i++)
d[i]=1;
for(i=2;i<=maxim;i++)
if(d[i]==1&&p[i]==0)
for(j=i;j<=maxim;j+=i){
d[j]*=i*(1-p[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;
}