A098388 a(n) = floor(log_2(prime(n))).

1, 1, 2, 2, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9

Offset

1,3

Comments

a(n) is the greatest k such that 2^k does not exceed prime(n). - David James Sycamore, Sep 14 2021

Links

T. D. Noe, Table of n, a(n) for n = 1..10000

Formula

a(n) = A000523(A000040(n)); A098391(n) = A000523(a(n)).

a(n) = A035100(n) - 1. - Michel Marcus, Sep 17 2017

Maple

map(ilog2, select(isprime, [2, seq(2*i+1, i=1..1000)])); # Robert Israel, Jun 08 2015

Mathematica

Floor[Log[2, Prime[Range[105]]]] (* data *) (* parameter changed by Hartmut F. W. Hoft, Jun 02 2015 *)

Prog

(PARI) a(n) = logint(prime(n), 2); \\ Michel Marcus, Sep 17 2017
(Python)
from sympy import prime
def A098388(n): return prime(n).bit_length()-1 # Chai Wah Wu, Nov 19 2024

Crossrefs

Cf. A000040, A000523, A023196, A035100, A098391.

Sequence in context: A070319 A057142 A320837 * A371672 A157792 A171896

Adjacent sequences: A098385 A098386 A098387 * A098389 A098390 A098391

Keyword

nonn,easy

Author

Reinhard Zumkeller, Sep 06 2004

Status

approved