OFFSET
1,1
COMMENTS
Corresponding values of k: 625678,704613,717718,1182168,3147353,4869813,5339178,5363578,5411562,5846777,7848283,7970403,8152962,9220303,9727978.
MATHEMATICA
apQ[p_]:=Module[{q=4p^2+1, r}, r=4q^2+1; AllTrue[{p, q, r, 4r^2+1}, PrimeQ]]; Select[ 4*Range[ 10^7]^2+1, apQ] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Aug 28 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Mar 05 2015
STATUS
approved