OFFSET
1,1
COMMENTS
Given a set of prime numbers P, finite or infinite, the set of numbers which are divisible by at least one of the primes in P has an asymptotic density Product_{p in P} (1 - 1/p). If P is finite, then this density is equal to 1/2 only when P = {2}. Otherwise, the density is 1/2 for infinitely many sets P. This sequence is the lexicographically earliest infinite sequence of such primes.
The first 5 terms are the Fermat primes (A019434).
a(10) = 7.455916... * 10^135 is too large to be included in the data section.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..12
FORMULA
a(1) = 3, a(n) = nextprime(r(n-1)/(r(n-1) - 1/2)), where r(n) = Product_{k=1..n-1} 1 - 1/a(n).
Product_{n=>1} (1 - 1/a(n)) = 1/2.
MATHEMATICA
s = {}; r = 1; p = 3; Do[AppendTo[s, p]; r *= 1 - 1/p; p = NextPrime[r/(r - 1/2)], {9}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 30 2020
STATUS
approved