A223037 a(n) = largest prime p such that Sum_{primes q = 2, ..., p} 1/q does not exceed n.

3, 271, 5195969, 1801241230056600467

Offset

1,1

Comments

Since Sum_{prime q} 1/q diverges, the sequence is infinite.

In fact, by the Prime Number Theorem Prime(k) ~ k log k as k -> infinity, and by integration Sum_{k <= n} 1/(k log k) ~ log log n, so a(n) ~ Prime(Floor(e^e^n)).

a(4) = A000040(A046024(4)-1) = Prime[43922730588128389], but Mathematica 7.0.0 cannot compute this prime on a Mac computer running OS X.

Instead, using a(4) = largest prime < A016088(4) = 1801241230056600523, Mathematica's PrimeQ function finds that a(4) = 1801241230056600467.

See A016088 for other relevant comments, references, links, and programs.

Links

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

Formula

a(n) = A000040(A046024(n)-1) = largest prime < A016088(n).

a(n) ~ Prime(Floor(e^e^n)) = A000040(A096232(n)) as n -> infinity.

Example

a(1) = 3 because 1/2 + 1/3 < 1 < 1/2 + 1/3 + 1/5.

Crossrefs

Cf. A016088, A046024, A024451, A096232.

Sequence in context: A230373 A003761 A216471 * A171358 A115477 A051365

Adjacent sequences: A223034 A223035 A223036 * A223038 A223039 A223040

Keyword

nonn,hard,more

Author

Jonathan Sondow, Apr 16 2013

Status

approved