login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A339344
Lexicographically earliest sequence of odd primes such that the asymptotic density of the numbers which are divisible by at least one of these primes is 1/2.
2
3, 5, 17, 257, 65537, 4294967311, 1229782942255939601, 88962710886098567818446141338419231, 255302062200114858892457591448999891874349780170241684791167583265041
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
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