cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

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.

This page as a plain text file.
%I A339344 #13 Oct 12 2023 09:51:27
%S A339344 3,5,17,257,65537,4294967311,1229782942255939601,
%T A339344 88962710886098567818446141338419231,
%U A339344 255302062200114858892457591448999891874349780170241684791167583265041
%N 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.
%C A339344 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.
%C A339344 The first 5 terms are the Fermat primes (A019434).
%C A339344 a(10) = 7.455916... * 10^135 is too large to be included in the data section.
%H A339344 Amiram Eldar, <a href="/A339344/b339344.txt">Table of n, a(n) for n = 1..12</a>
%F A339344 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).
%F A339344 Product_{n=>1} (1 - 1/a(n)) = 1/2.
%t A339344 s = {}; r = 1; p = 3; Do[AppendTo[s, p]; r *= 1 - 1/p; p = NextPrime[r/(r - 1/2)], {9}]; s
%Y A339344 Cf. A000215, A019434, A080307, A125045, A262228, A339345.
%K A339344 nonn
%O A339344 1,1
%A A339344 _Amiram Eldar_, Nov 30 2020