OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
MATHEMATICA
t = {}; Do[If[DivisorSigma[1, n] == EulerPhi[n]*PrimeNu[n], AppendTo[t, n]], {n, 700000}]; t (* Jayanta Basu, Jul 01 2013 *)
Select[Range[15*10^5], DivisorSigma[1, #]==EulerPhi[#]PrimeNu[#]&] (* Harvey P. Dale, Mar 04 2023 *)
PROG
(PARI) for(n=1, 10000000, if(sigma(n)==eulerphi(n)*omega(n), print1(n, ", "))) (Klasen)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Aug 31 2002
EXTENSIONS
More terms from Lambert Klasen (lambert.klasen(AT)gmx.de), Dec 23 2004
STATUS
approved