OFFSET
1,2
COMMENTS
a(22) > 10^10. - _Amiram Eldar_, Jan 27 2019
EXAMPLE
sigma(96) = 252, phi(96) = 32, d(96) = 12 and phi(252 - 32 - 12) = 96.
MAPLE
with(numtheory): P:=proc(q) local n;
for n from 1 to q do
if phi(sigma(n)-phi(n)-tau(n))=n then print(n);
fi; od; end: P(10^10);
MATHEMATICA
Select[Range[10^4], EulerPhi[DivisorSigma[1, #] - EulerPhi[#] - DivisorSigma[0, #]] == # &] (* _Alonso del Arte_, Jun 18 2014 *)
PROG
(PARI) isok(n) = (spd = sigma(n) - eulerphi(n) - numdiv(n)) && (eulerphi(spd) == n); \\ _Michel Marcus_, Jun 18 2014
CROSSREFS
KEYWORD
nonn,more
AUTHOR
_Paolo P. Lava_, Jun 18 2014
EXTENSIONS
a(15)-a(16) from _Jens Kruse Andersen_, Jul 23 2014
a(17)-a(21) from _Amiram Eldar_, Jan 27 2019
STATUS
approved
