OFFSET
0,1
COMMENTS
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.
Lowell Schoenfeld, Corrigendum: "Sharper bounds for the Chebyshev functions theta (x) and psi (x). II" (Math. Comput. 30 (1976), number 134, 337-360), Math. Comp. 30 (1976), number 136, 900.
FORMULA
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
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