|
%I #23 Dec 16 2023 20:11:18
%S 131,151,229,233,311,571,643,727,941,1013,1051,1153,1373,1531,1667,
%T 1669,1723,1783,1787,1831,1951,1979,2029,2131,2213,2239,2311,2441,
%U 2593,2621,2633,2659,2663,2887,3001,3011,3019,3121,3169,3209,3253,3347,3413,3457
%N Primes prime(k) such that all integers in [(prime(k-1)+1)/2,(prime(k)-1)/2] are composite, excluding those primes in A080359.
%C The primes of A080359 larger than 3 all have the property that the integers in the interval selected by halving the value of the preceding prime and halving their own value are all composite. This sequence here collects the primes that are not in A080359 but still share this property of the prime-free subinterval.
%H Jean-François Alcover, <a href="/A164294/b164294.txt">Table of n, a(n) for n = 1..1106</a>
%H V. Shevelev, <a href="http://arXiv.org/abs/0908.2319">On critical small intervals containing primes</a>, arXiv:0908.2319 [math.NT], 2009. [From _Vladimir Shevelev_, Aug 20 2009]
%H V. Shevelev, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL15/Shevelev/shevelev19.html">Ramanujan and Labos primes, their generalizations, and classifications of primes</a>, J. Integer Seq. 15 (2012) Article 12.5.4.
%F A164333 \ A080359.
%e For the prime 1531=A000040(242), the preceding prime is A000040(241)=1523, and the integers from (1523+1)/2 = 762 up to (1531-1)/2 = 765 are all composite, as they fall in the gap between A000040(135) and A000040(136). In addition, 1531 is not in A080359, which adds 1531 to this sequence here.
%t maxPrime = 3500;
%t kmax = PrimePi[maxPrime];
%t A164333 = Select[Table[{(Prime[k - 1] + 1)/2, (Prime[k] - 1)/2}, {k, 3, kmax}], AllTrue[Range[#[[1]], #[[2]]], CompositeQ] &][[All, 2]]*2 + 1;
%t b[1] = 2; b[n_] := b[n] = Module[{k = b[n - 1]}, While[(PrimePi[k] - PrimePi[Quotient[k, 2]]) != n, k++]; k];
%t A080359 = Reap[For[n = 1, b[n] <= maxPrime, n++, Sow[b[n]]]][[2, 1]];
%t Complement[A164333, A080359] (* _Jean-François Alcover_, Sep 14 2018 *)
%o (PARI) okprime(p) = { my(k = primepi(p)); for (i = (prime(k-1)+1)/2, (prime(k)-1)/2, if (isprime(i), return (0));); return (1);}
%o lista(nn) = {vlp = readvec("b080359.txt"); forprime (p=2, nn, if (! vecsearch(vlp, p) && okprime(p), print1(p, ", ")););} \\ _Michel Marcus_, Jan 15 2014
%Y Cf. A080359, A104272, A164288, A001262, A001567, A062568, A141232, A164368.
%K nonn
%O 1,1
%A _Vladimir Shevelev_, Aug 12 2009
%E Extended beyond 571 by _R. J. Mathar_, Oct 02 2009
|
