OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
For n=9: sigma(9)=13, phi(9)=6, gcd(13,6)=1, core(9)=3, sigma(3)=4, phi(3)=2, gcd(4,2)=2.
MATHEMATICA
ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] cor[x_] := Apply[Times, ba[x]] g1[x_] := GCD[DivisorSigma[1, x], EulerPhi[x]] g2[x_] := GCD[DivisorSigma[1, cor[x]], EulerPhi[cor[x]]] Do[s1=g1[n]; s2=g2[n]; If[Greater[s2, s1], Print[n]], {n, 1, 256}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Oct 17 2002
STATUS
approved