OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
MAPLE
select(n -> numtheory:-sigma(n+1)=3*numtheory:-phi(n), [$1..10000]); # Robert Israel, May 29 2018
MATHEMATICA
Select[Range[6000], DivisorSigma[1, #+1]==3*EulerPhi[#]&] (* Harvey P. Dale, Jun 11 2017 *)
PROG
(PARI) isok(n) = sigma(n+1) == 3*eulerphi(n); \\ Michel Marcus, May 29 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 21 2002
STATUS
approved
