OFFSET
1,1
LINKS
Zak Seidov, Table of n, a(n) for n = 1..503 (all terms up to 10^10)
MATHEMATICA
Select[Prime[Range[132*10^5]], AllTrue[(#^Range[2, 8, 2]+5)/6, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, May 11 2018 *)
PROG
(PARI) forprime(p=1, 10^7, k=0; for(i=1, 4, P=(p^(2*i)+5)/6; if(P\1==P, if(ispseudoprime(P), k++); if(!ispseudoprime(P), k=0; break)); if(P\1!=P, k=0; break)); if(k, print1(p, ", "))) \\ Derek Orr, Jan 21 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jan 21 2015
STATUS
approved