OFFSET
1,1
EXAMPLE
For n = 13954, phi(phi(n)) + sigma(sigma(n)) - phi(sigma(n)) - sigma(phi(n)) = phi(6976) + sigma(20934) - phi(20934) - sigma(6976) = 3456 + 45396 - 6972 - 13970 = 27910 = 6976 + 20934 = phi(n) + sigma(n), so 13954 is in the sequence.
MATHEMATICA
g[x_] := Module[{a, b, c, d, e, f}, a=EulerPhi[x]; b=DivisorSigma[1, x]; c=EulerPhi[a]; d=DivisorSigma[1, b]; e=EulerPhi[b]; f=DivisorSigma[1, a]; a+b==c+d-e-f]; Do[If[g[n]==True, Print[n]], {n, 1, 10^6}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Jan 24 2002
EXTENSIONS
Edited by Dean Hickerson, Jan 26 2002
a(8)-a(22) from Donovan Johnson, Jan 02 2009
STATUS
approved