login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A223853 a(n) = ceiling(li(2*2^n) - li(2^n)) - (pi(2*2^n) - pi(2^n)) with li(x) the logarithmic integral and pi(x) the prime counting function. 2
1, 1, 2, 1, 2, 2, 2, 1, 3, 4, 1, 7, 1, 13, 10, 4, 25, -5, 49, 17, 38, 82, 103, -55, 245, 290, 105, 621, -107, 1219, 1196, -274, 1749, 5329, 2881, 2451, 6836, 2910, 15905, 28044, -10652, 55758, 18068, 129994, -95925, 52787, 443983, 253331, 151395, 740898, -352415 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

COMMENTS

This is the difference between the estimate for the number of primes in power of two intervals determined by the li approximation, and the actual number of primes in the power of two interval. The MAGMA program gives the ceiling of the difference between the li estimate at the end points of the interval and the actual number of primes in the interval (A036378).

H. J. J. te Riele (1987) using methods developed by Lehman (1966) showed that between 6.62*10^370 and 6.69*10^370 there are more than 10^180 consecutive integers where pi(x) > li(x). It is worth noting that this falls entirely within the power of two interval starting at 2^1231, and while the condition "li underestimates the number of primes in an interval" is not sufficient to imply that pi(x) > li(x), for example in (2^18, 2^19) li(x) underestimates by 5 but li(x) > pi(x) at every point in the interval, it does seem to be necessary for this to occur, assuming runs of consecutive values where pi(x) > li(x) do not cross a power of two.

LINKS

Brad Clardy, Table of n, a(n) for n = 1..74

R. S. Lehman, On the difference pi(x) - li(x), Acta Arithmetica XI (1966), p. 397-410

H. J. J. te Riele, On the sign of the difference pi(x) - li(x), Math. Comp. 48 (1987), p.323-328

FORMULA

a(n) = A223900(n) - A036378(n).

MATHEMATICA

pi = Table[PrimePi[2^n], {n, 1, 30}];

li = Table[LogIntegral[2^n], {n, 1, 30}];

Ceiling[Rest@li - Most@li] - (Rest@pi - Most@pi) (* Peter Luschny, Oct 14 2017 *)

PROG

(MAGMA)

1;

for i := 2 to 29 do

    x := 2^i;

    y := 2^(i+1);

    delta_li := Ceiling(LogIntegral(y) - LogIntegral(x));

    delta_pi := #PrimesInInterval(x, y);

    delta_li - delta_pi;

end for;

CROSSREFS

Cf. A000720, A036378, A052435, A223900.

Sequence in context: A068307 A158946 A303428 * A023645 A167865 A218654

Adjacent sequences:  A223850 A223851 A223852 * A223854 A223855 A223856

KEYWORD

sign

AUTHOR

Brad Clardy, Mar 28 2013

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 1 05:04 EDT 2020. Contains 333155 sequences. (Running on oeis4.)