%I #23 Oct 29 2023 17:55:47
%S 53,89,157,173,211,251,257,263,293,331,337,359,367,373,389,409,449,
%T 479,509,541,547,557,563,577,587,593,607,631,653,683,691,701,709,719,
%U 727,733,751,787,797,839,919,929,947,953,977,983,991,997,1039,1069,1103,1109,1117
%N Primes p with property that (p - previous prime) >= 6 and (next prime - p) >= 6.
%H Robert Israel, <a href="/A137869/b137869.txt">Table of n, a(n) for n = 1..10000</a>
%p M:=1000; t1:=[];
%p for i from 2 to M do
%p p:=ithprime(i); o:=prevprime(p); q:=nextprime(p);
%p if p-o >= 6 and q-p >= 6 then t1:=[op(t1),p]; fi; od:
%p t1; # _N. J. A. Sloane_
%t lst={};Do[p=Prime[n];If[ !PrimeQ[p-2]&&!PrimeQ[p+2]&&!PrimeQ[p-4]&&!PrimeQ[p+4],AppendTo[lst,p]],{n,6!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Jul 20 2009 *)
%t Select[Partition[Prime[Range[200]],3,1],Min[Differences[#]]>5&][[;;,2]] (* _Harvey P. Dale_, Oct 29 2023 *)
%o (PARI) p=q=2;forprime(r=2,999,r-q>4 & q-p>4 & print1(q","); p=q; q=r) \\ _M. F. Hasler_, May 02 2009
%o (Magma) [p:p in PrimesInInterval(3,1200)|p-PreviousPrime(p) ge 6 and NextPrime(p)-p ge 6]; // _Marius A. Burtea_, Aug 11 2019
%Y These are the primes not in A167773.
%Y Cf. A053070.
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Apr 29 2008
%E Edits and more terms from _N. J. A. Sloane_, May 02 2009