OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3)=397 is in the sequence because it is prime, the next prime is 401, and 2*397^2-401^2 = 154417 and 2*401^2-397^2 = 163993 are prime.
MAPLE
N:= 100: # for the first N terms
q:= 2: count:= 0: R:= NULL:
for iter from 1 while count < N do
p:= q; q:= nextprime(q);
if isprime(2*p^2-q^2) and isprime(2*q^2-p^2) then
count:= count+1; R:= R, p;
fi
od:
R;
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Nov 11 2020
STATUS
approved