OFFSET
1,2
COMMENTS
Does the ratio of adjacent terms converge?
It appears that lim_{n->infinity} a(n+1)/a(n) = 1.7532... - Jon E. Schoenfield, Feb 21 2019
For n > 1, a(n) is smallest prime p = prime(k) such that no fewer than (n-1)/n of any p# consecutive integers are divisible by a prime not greater than p. Cf. A053144(k)/A002110(k). - Peter Munn, Apr 29 2017
Also, the smallest prime p such that the sum of the reciprocals of the p-smooth numbers converges to at least n. - Keith F. Lynch, Apr 29 2023
Also, if m is a random integer much larger than the square of a(n), and m is not divisible by any prime less than or equal to a(n), the probability that m is prime is n/log(m). - Keith F. Lynch, Dec 17 2023
LINKS
Jon E. Schoenfield, Table of n, a(n) for n = 1..44
Eric Weisstein's World of Mathematics, Totient Function
Eric Weisstein's World of Mathematics, Primorial
EXAMPLE
7#/phi(7#) = (2*3*5*7)/(1*2*4*6) = 4.375 >= 4, 5#/phi(5#) = 3.75. Hence a(4) = 7.
MATHEMATICA
prod=1; i=0; Table[While[prod<n, i++; prod=prod/(1-1/Prime[i])]; Prime[i], {n, 1, 20}]
PROG
(PARI) al(lim) = local(mm, n, m); mm=3; n=2; m=1; forprime(x=3, lim, n*=x; m*= (x-1); if (n\m >= mm, print1(x", "); mm++)); /* This will generate all terms of this sequence from the 3rd onward, up to lim. The computation slows down for large values because of the size of the internal values. */ \\ Fred Schneider, Aug 13 2009, modified by Franklin T. Adams-Watters, Aug 29 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Jan 09 2004
EXTENSIONS
More terms from David W. Wilson, Sep 28 2005
Sequence reference in name corrected by Peter Munn, Apr 29 2017
STATUS
approved