1,2
(PARI) nud(n) = 2^omega(n); pi(n) = s=0; for(x=1, n, if(isprime(x), s++)); s; sopf(n, s, fac, i)=fac=factor(n); for(i=1, matsize(fac)[1], s=s+fac[i, 1]); return(s); for(n=1, 10^6, if(sopf(n)*nud(n)==pi(n), print(n)))
Cf. A008472, A034444, A000720.
Sequence in context: A115693 A112403 A043932 * A044576 A104453 A143945
Adjacent sequences: A064012 A064013 A064014 * A064016 A064017 A064018
more,nonn
Jason Earls (zevi_35711(AT)yahoo.com), Sep 06 2001