Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A223037
a(n) = largest prime p such that Sum_{primes q = 2, ..., p} 1/q does not exceed n.
3
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.
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
KEYWORD
nonn,hard,more
AUTHOR
Jonathan Sondow, Apr 16 2013
STATUS
approved