login
a(n) is the smallest prime p such that, for m >= nextprime(p), there are more composites than primes in the range [2, m], where multiples of primes prime(1) through prime(n) are excluded.
0

%I #11 Jun 16 2023 13:43:39

%S 13,113,1069,5051,18553,44417,99439,190921,356351,603149,933073,

%T 1416223,2044201,2856559,3957883,5379287,7093217,9113263,11693687,

%U 14701529,18345209,22758829,27879563,33938257,40808759,48364003,57099061,67292237,78919781,92417891

%N a(n) is the smallest prime p such that, for m >= nextprime(p), there are more composites than primes in the range [2, m], where multiples of primes prime(1) through prime(n) are excluded.

%e The number of primes, N_p, and the number of composite, N_c, in the range [2, m] are listed in the table below, where N_p = N_c occurs at m = 9, 11 and 13. For m >= nextprime(13) = 17, N_c > N_p. So, a(0) = 13 is the case for n = 0, in which none of the multiples of primes is excluded from the integer list.

%e m: 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, ...

%e N_p: 1, 2, 2, 3, 3, 4, 4, 4, 4, 5, 5, 6, 6, 6, 6, 7, ...

%e N_c: 0, 0, 1, 1, 2, 2, 3, 4, 5, 5, 6, 6, 7, 8, 9, 9, ...

%e If the multiples of prime(1) are excluded from the list, 113 is the smallest prime such that N_c > N_p for m >= nextprime(113) = 127 and, thus, a(1) = 113 (see below).

%e m: 3, 5, 7, ..., 93, 95, 97, 99, 101, 103, 105, 107, 109, 111, 113, 115, ...

%e N_p: 1, 2, 3, ..., 23, 23, 24, 24, 25, 26, 26, 27, 28, 28, 29, 29, ...

%e N_c: 0, 0, 0, ..., 23, 24, 24, 25, 25, 25, 26, 26, 26, 27, 27, 28, ...

%e If multiples of prime(1) and prime(2) are excluded, a(2) = 1069. If multiples of prime(1), prime(2) and prime(3) are excluded, a(3) = 5051.

%o (Python)

%o from sympy import isprime, prime

%o R = []; L = [x for x in range(2, 100000001)]

%o for n in range(30):

%o np = 0; nc = 0; found = 0

%o if n > 0: q = prime(n); L = [x for x in L if x%q != 0]

%o for m in L:

%o if isprime(m): np += 1; p = m

%o else: nc += 1

%o if np == nc: Lp = p; found = 1

%o if found: R.append(Lp)

%o print(*R, sep = ", ")

%Y Cf. A000040, A002808, A072731, A097454.

%K nonn

%O 0,1

%A _Ya-Ping Lu_, Mar 29 2023