A334598 a(n) is the largest nonnegative integer m such that m >= pi(m)^(1 + 1/n).

4, 28, 1860, 149052, 12771496, 1221908916, 132662942122, 16354869261256, 2272946910544740, 353076161059625536, 60799066209732571716, 11518836088596729968092

Offset

1,1

Comments

For a nonnegative integer m, pi(m) = A000720(m). It is well-known that if

m >= 17, then m/log(m) < pi(m). [Rosser and Schoenfeld]

Fix a real exponent d > 0. If m is big enough, then m < (m/log(m))^(1 + d). In particular, choosing d = 1/n, with n >= 1, we deduce that a(n) exists.

Note that different choices of the exponent d will produce analogous sequences.

The estimates of pi(m) in [Dusart, Thm. 5.1] and [Axler, Thm. 2] allow us to obtain upper and lower bounds for a(n). In particular, we can conclude that in base 10:

a(13) has 25 digits, starting with 238;

a(14) has 27 digits, starting with 536;

a(15) has 30 digits, starting with 1304;

a(16) has 32 digits, starting with 3409.

The tool primecount [Walisch], used to compute pi(10^28) in A006880, can handle pi(m) for m <= 10^31, and since (a(n)) is monotonically increasing, it seems that the computation of a(n) for n >= 16 will be challenging.

It is easy to see that for every n >= 1, a(n) is even and a(n)+1 is prime. - Eduard Roure Perdices, Nov 07 2021

Links

Table of n, a(n) for n=1..12.

Christian Axler, Estimates for pi(x) for large values of x and Ramanujan's prime counting inequality, Integers 18 (2018), Paper No. A61, 14 pp.

Pierre Dusart, Explicit estimates of some functions over primes, The Ramanujan Journal 45 (2018), no. 1, 227-251.

J. Barkley Rosser and Lowell Schoenfeld, Approximate formulas for some functions of prime numbers, Illinois J. Math. 6 (1962), no. 1, 64-94.

Kim Walisch, primecount, Github, Aug 14 2021.

Crossrefs

Cf. A334599, A087235, A038625.

Sequence in context: A336622 A354842 A084594 * A000838 A218174 A220756

Adjacent sequences: A334595 A334596 A334597 * A334599 A334600 A334601

Keyword

nonn,more,hard

Author

Eduard Roure Perdices, May 07 2020

Extensions

a(8) from Giovanni Resta, May 07 2020

a(9)-a(10) from Daniel Suteu, May 20 2020

a(11)-a(12) from Eduard Roure Perdices, Nov 07 2021

Status

approved