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