OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..2000
EXAMPLE
a(3) = 51071 is in the sequence because p=51071, q=51109, r=51131 are consecutive primes such that (3*p+q)/2=102161, (p+3*q)/2=102199, (3*q+r)/2=102229, (q+3*r)/2=102251 are all prime.
MAPLE
q:= 3: r:= 5:
count:= 0: Res:= NULL:
while count < 100 do
p:= q; q:= r; r:= nextprime(q);
if isprime((3*p+q)/2) and isprime((p+3*q)/2) and isprime((3*q+r)/2)
and isprime((q+3*r)/2) then
count:= count+1; Res:= Res, p;
fi
od:
Res;
MATHEMATICA
Select[Partition[Prime[Range[175000]], 3, 1], AllTrue[{(3#[[1]]+#[[2]])/2, (#[[1]]+ 3#[[2]])/2, (3#[[2]]+#[[3]])/2, (#[[2]]+3#[[3]])/2}, PrimeQ]&][[;; , 1]] (* Harvey P. Dale, Mar 19 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 26 2020
STATUS
approved