A073798 pi(n) is a power of 2, where pi(n) = A000720(n) is the number of primes <= n.
2, 3, 4, 7, 8, 9, 10, 19, 20, 21, 22, 53, 54, 55, 56, 57, 58, 131, 132, 133, 134, 135, 136, 311, 312, 719, 720, 721, 722, 723, 724, 725, 726, 1619, 1620, 3671, 3672, 8161, 8162, 8163, 8164, 8165, 8166, 17863, 17864, 17865, 17866, 17867, 17868, 17869, 17870
Offset: 1
Keywords
Examples
10 is in the sequence since pi(10)=4=2^2.
Links
- Chai Wah Wu, Table of n, a(n) for n = 1..1103 (n = 1..665 from Ivan Neretin)
Programs
-
Mathematica
pow2[n_] := n==1||(n>1&&IntegerQ[n/2]&&pow2[n/2]); Select[Range[20000], pow2[PrimePi[ # ]]&] Flatten@Table[Range[p = Prime[2^k], NextPrime[p] - 1], {k, 0, 11}] (* Ivan Neretin, Jan 21 2017 *)
-
PARI
isok(n) = my(pi = primepi(n)); (pi==1) || (pi==2) || (ispower(primepi(n),,&k) && (k==2)); \\ Michel Marcus, Jan 23 2017
Extensions
Edited by Dean Hickerson, Aug 15 2002
Comments