OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
610 is in A249667: the least prime>610 is 613, and 610+613=1223 is prime; the largest prime<610 is 607, and 610+607=1217 is prime. Also, 613-610=610-607=3, then 610 is in the current sequence.
MAPLE
P:= select(isprime, [seq(i, i=3..12000, 2)]):
A:= (P[1..-2] + P[2..-1])/2:
A[select(i -> isprime(A[i]+P[i]) and isprime(A[i]+P[i+1]), [$1..nops(A)])]; # Robert Israel, Oct 03 2025
PROG
(PARI) {for(i=3, 2*10^4, m=nextprime(i+1); k=i+m; n=precprime(i-1); q=i+n; if(isprime(k)&&isprime(q)&&m-i==i-n, print1(i, ", ")))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Antonio Roldán, Dec 03 2014
STATUS
approved
