OFFSET
1,1
COMMENTS
Suggested in an email by J. M. Bergot.
Primes prime(k) such that prime(k)^2 + prime(k+1)^2 is in A051448.
LINKS
Robert Israel, Table of n, a(n) for n = 1..2000
EXAMPLE
a(3) = 23 is a term because 23 and 29 are consecutive primes with 23^2 + 29^2 = 1370 = 2*5*137, and 2+5+137 = 144 = 12^2.
MAPLE
q:= 2: R:= NULL: count:=0:
while count < 100 do
p:= q; q:= nextprime(p);
s:= p^2 + q^2;
if issqr(add(t[1]*t[2], t = ifactors(s)[2])) then
R:= R, p; count:= count+1;
fi;
od:
R;
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Jan 01 2023
STATUS
approved