login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A016088 a(n) = smallest prime p such that Sum_{primes q = 2, ..., p} 1/q exceeds n. 24
2, 5, 277, 5195977, 1801241230056600523 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,1
COMMENTS
The indices of these primes are in A046024: a(n) = A000040(A046024(n)).
Comment from Eric Bach, Jun 03 2005: Schoenfeld (Math. Comp. 1976) has explicit estimates that imply, assuming the Riemann hypothesis, that the sum first exceeds 4 for some x in the range (1.80124093... * 10^18, 1.80124152... * 10^18).
Eric Bach (Sep 14 2005) comments that the next element in the sequence is about 4.2 * 10^49, so it may remain unknown for all eternity.
Comment from Richard C. Schroeppel, Nov 09 2006: The Bach-Sorenson algorithm takes around x^(1/3) space and x^(2/3) time. When x = 4 * 10^49, these are roughly 10^16 and 10^33. We have facilities today that handle 10^16 storage. Current world computing capability is ~10^25 instructions/year (10^8.5 machines, 10^9.5 inst/sec, 10^7.5 sec/year). The algorithm seems very parallelizable. So with current resources, we could have the value for a(5) in a mere 10^8 years. This might seem like a long time, but it's far short of eternity.
The sequence is less than 2^3^n for all n >= 1. Moreover, the limit a of a_n^e^(-n) seems to exist and is approximately 2.16 and thus a^e^n is an estimate for the sequence which is not completely wrong. - Wolfgang Burmeister (Wolfgang-Burmeister(AT)t-online.de), May 05 2007
Sequence can be approximated by the simple expression a(n) = exp(exp(n-0.2615)) due to the behavior of the sum of reciprocals of primes. This gives: a(4)=1.801..*10^18; a(5)=4.2..*10^49 and a(6)=7.7..*10^134. - Carmine Suriano, Mar 25 2014
REFERENCES
Calvin C. Clawson, Mathematical Mysteries, The Beauty and Magic of Numbers, Plenum Press, NY and London, 1996, page 64.
LINKS
E. Bach, D. Klyve, and J. P. Sorenson, Computing prime harmonic sums, Math. Comp. 78 (268) (2009) 2283-2305.
Eric Bach and Jonathan Sorenson, Computing Prime Harmonic Sums [Wayback Machine cached version]
J. Barkley Rosser and Lowell Schoenfeld, Sharper bounds for the Chebyshev functions theta(x) and psi(x), Collection of articles dedicated to Derrick Henry Lehmer on the occasion of his seventieth birthday. Math. Comp. 29 (1975), 243-269.
Lowell Schoenfeld, Sharper bounds for the Chebyshev functions theta (x) and psi (x). II. Math. Comp. 30 (1976), number 134, 337-360.
FORMULA
From Jonathan Sondow, Apr 17 2013: (Start)
a(n) = A000040(A000720(A223037(n))+1).
a(n) ~ prime(floor(e^e^n)) = A000040(A096232(n)) as n -> oo (see comment in A223037). (End)
MATHEMATICA
s = 0; k = 1; Do[ While[ s = N[ s + 1/Prime[ k ], 36 ]; s <= n, k++ ]; Print[ Prime[ k ] ]; k++, {n, 1, 3} ]
s = 0; n = 0; For[k = 1, k > 0, k++, If[(s = N[s + 1/(p = Prime[k]), 40]) > n, Print[p|s]; n++ ]] (* Wolfgang Burmeister (Wolfgang-Burmeister(AT)t-online.de), May 05 2007 *)
PROG
(PARI) a(n)=my(t); forprime(p=2, , t+=1./p; if(t>n, return(p))) \\ Charles R Greathouse IV, Apr 29 2015
CROSSREFS
Sequence in context: A309675 A042341 A300900 * A248961 A353050 A308270
KEYWORD
nonn,nice,hard,more
AUTHOR
EXTENSIONS
a(0) from Wolfgang Burmeister (Wolfgang-Burmeister(AT)t-online.de), May 05 2007
a(3) corrected by Ulrich Schimke (UlrSchimke(AT)aol.com)
a(4) computed by Eric Bach and Jon Sorenson, Sep 14 2005. They used a variant of the Lagarias-Miller-Odlyzko algorithm for pi(x) and found that sum_{p <= 1801241230056600467} 1/p = 3.99999999999999999966 and sum_{p <= 1801241230056600523} 1/p = 4.00000000000000000021. There are no primes between 1801241230056600467 and 1801241230056600523. Total computing time was about two weeks, divided between two workstations (i.e., about a week on each).
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:13 EDT 2024. Contains 371235 sequences. (Running on oeis4.)