OFFSET
1,1
COMMENTS
From Dean Hickerson, Jan 13 2003: (Start)
Suppose you have a list of the first n prime numbers p_1, ..., p_n and you want to estimate the next one. The probability that a random integer is not divisible by any of p_1, ..., p_n is (1-1/p_1) * ... * (1-1/p_n). In other words, 1 out of every 1/((1-1/p_1) * ... * (1-1/p_n)) integers is relatively prime to p_1, ..., p_n.
So we might expect the next prime to be roughly this much larger than p_n; i.e. p_(n+1) may be about p_n + 1/((1-1/p_1) * ... * (1-1/p_n)). This sequence and A003067, A003068 are obtained by replacing this approximation by an exact equation, using 3 different ways of making the results integers. (End)
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
M. D. Hirschhorn, How unexpected is the prime number theorem?, Amer. Math. Monthly, 80 (1973), 675-677.
M. D. Hirschhorn, How unexpected is the prime number theorem?, Amer. Math. Monthly, 80 (1973), 675-677. [Annotated scanned copy]
R. C. Vaughan, The problime number theorem, Bull. London Math. Soc., 6 (1974), 337-340.
MAPLE
a[1] := 2: for i from 1 to 150 do a[i+1] := floor(a[i]+1/product((1-1/a[j]), j=1..i)): od: # James A. Sellers, Mar 07 2000
MATHEMATICA
a[1] = 2; a[n_] := a[n] = Floor[a[n-1] + 1/Product[1-1/a[j], {j, 1, n-1}]]; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Mar 09 2012, after James A. Sellers *)
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
EXTENSIONS
More terms from James A. Sellers, Mar 07 2000
STATUS
approved