OFFSET
1,1
EXAMPLE
For n = 87380 = 4*5*17*257, gpf(n) = 257, phi(n) = 65536, commutator[87380] = phi(257) - gpf(65536) = 256 - 2 = 254.
MATHEMATICA
pf[x_] := Part[Reverse[Flatten[FactorInteger[x]]], 2] Do[s=EulerPhi[pf[n]]-pf[EulerPhi[n]]; If[Equal[s, 254], Print[{n, n/257, pf[n/257]}]], {n, 3, 1000000}] (* Terms of sequence are n *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, May 09 2002
STATUS
approved