login
This site is supported by donations to The OEIS Foundation.

 

Logo

Annual Appeal: Today, Nov 11 2014, is the 4th anniversary of the launch of the new OEIS web site. 70,000 sequences have been added in these four years, all edited by volunteers. Please make a donation (tax deductible in the US) to help keep the OEIS running.

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. 18
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 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 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

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

E. Bach, D. Klyve, J. P. Sorenson, Computing prime harmonic sums, Math. Comp. 78 (268) (2009) 2283-2305.

Eric Bach and Jonathan Sorenson, Computing Prime Harmonic Sums

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.

Index entries for sequences related to decimal expansion of 1/n

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 -> infinity (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 *)

CROSSREFS

Cf. A046024, A223037.

Sequence in context: A121270 A085603 A042341 * A248961 A042909 A183129

Adjacent sequences:  A016085 A016086 A016087 * A016089 A016090 A016091

KEYWORD

nonn,nice,hard,more

AUTHOR

Robert G. Wilson v

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 | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified November 24 04:08 EST 2014. Contains 249867 sequences.