|
PROGRAM
|
(PARI) a=vector(10^5); a[1]=1; for(n=2, 10^5, a[n]=a[n-1]+n*abs(moebius(n))); forprime(p=2, 10^5, forprime(q=p+1, 10^5, n=a[q]-a[p-1]; if((n%q==0)&&(n%p==0), k=n/(p*q); u=factor(k); w=matsize(u); if((w[1]==0)||(u[1, 1]>=p)&&(u[w[1], 1]<=q), print(n))))) - Robert Gerbicz (robert.gerbicz(AT)gmail.com), Jun 09 2007
|