%I #23 Jul 05 2022 23:14:48
%S 3,5,7,9,11,13,17,19,23,25,29,31,37,41,43,47,49,53,59,61,67,71,73,79,
%T 81,83,89,97,101,103,107,109,113,121,127,131,137,139,149,151,157,163,
%U 167,169,173,179,181,191,193,197,199,211,223,227,229,233,239,241
%N Odd numbers with prime number of divisors.
%C Here but not in A062090: [729, 15625, 59049, 117649, 531441]; in A062090 but not here: [1, 6561, 390625]. - _Klaus Brockhaus_, Nov 01 2001
%H T. D. Noe, <a href="/A050150/b050150.txt">Table of n, a(n) for n = 1..1000</a>
%F Numbers of the form p^e where p is an odd prime and e+1 is a prime.
%F A010051(a100995(a(n)) + 1) = 1. - _Reinhard Zumkeller_, Aug 16 2013
%F a(n) ~ n log n. - _Charles R Greathouse IV_, Aug 28 2013
%e Numbers of the form p^6 for example (such as 3^6 = 729) are here but not in A062090.
%t Select[ Range[1, 250, 2], PrimeQ[ Length[ Divisors[ # ]]] & ]
%t Select[Range[1,799,2],PrimeQ[DivisorSigma[0,#]]&] (* _Harvey P. Dale_, Jun 22 2011 *)
%o (PARI) forstep(n=1,1000,2, if(isprime(numdiv(n)),print1(n,",")))
%o (PARI) is(n)=n%2 && isprime(isprimepower(n)+1) \\ _Charles R Greathouse IV_, Aug 28 2013
%o (Haskell)
%o a050150 n = a050150_list !! (n-1)
%o a050150_list = filter ((== 1) . a010051 . (+ 1) . a100995) [1, 3 ..]
%o -- _Reinhard Zumkeller_, Aug 16 2013
%o (Python)
%o from sympy import divisor_count, isprime
%o def ok(n): return n and n%2 and isprime(divisor_count(n))
%o print([k for k in range(250) if ok(k)]) # _Michael S. Branicky_, Jul 05 2022
%Y Cf. A062090 (a different sequence).
%K easy,nonn,nice
%O 1,1
%A _Jason Earls_, Jul 04 2001
%E More terms from _Jud McCranie_, Oct 31 2001
|