OFFSET
1,2
MAPLE
with(numtheory): A056013:=n->`if`(sigma(n)-tau(n)+phi(n) mod n = 0, n, NULL): seq(A056013(n), n=1..10^5); # Wesley Ivan Hurt, Apr 26 2017
MATHEMATICA
Do[If[Mod[DivisorSigma[1, n]-DivisorSigma[0, n]+EulerPhi[n], n]==0, Print[n]], {n, 1, 10^7}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jul 24 2000
EXTENSIONS
a(6)-a(8) from Donovan Johnson, Nov 12 2012
STATUS
approved