OFFSET
0,2
COMMENTS
The sequence is obtained by including all those numbers between primorial(n) and primorial(n+1) which have the largest prime divisor < = prime(n).
LINKS
Matthieu Pluntz, Table of n, a(n) for n = 0..34184
MATHEMATICA
{{1}}~Join~Array[Function[{p, q}, Select[Range[p, p NextPrime[q] - 1], FactorInteger[#][[-1, 1]] <= q &]] @@ {Product[Prime@ i, {i, #}], Prime[#]} &, 4] // Flatten (* Michael De Vlieger, Nov 18 2017 *)
PROG
(PARI) v = vector(1000); pr = 1; forprime(p = 2, 1000, pr *= p; v[p] = pr); for (n = 2, 1000, f = factor(n); p = f[matsize(f)[1], 1]; if (n >= v[p], print(n))); \\ David Wasserman, Feb 22 2006
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Dec 31 2003
EXTENSIONS
More terms from David Wasserman, Feb 22 2006
STATUS
approved