OFFSET
1,1
COMMENTS
Subsequence of A052291.
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..10000 (first 1000 terms from Vincenzo Librandi)
MAPLE
A122424:=n->`if`(isprime(n) and isprime(4*n^2+1) and isprime(4*(4*n^2+1)^2+1), n, NULL): seq(A122424(n), n=1..10^5); # Wesley Ivan Hurt, Aug 04 2014
MATHEMATICA
Select[Prime[Range[3500]], PrimeQ[4 #^2 + 1] && PrimeQ[64 #^4 + 32 #^2 + 5]&] (* Vincenzo Librandi, Apr 09 2013 *)
PROG
(Magma) [p: p in PrimesUpTo(30000) | IsPrime(q) and IsPrime(4*q^2+1) where q is 4*p^2+1]; // Vincenzo Librandi, Apr 09 2013
(PARI)
f(x)=4*x^2+1;
forprime(p=1, 10^5, if(isprime(f(p))&&isprime(f(f(p))), print1(p, ", "))) \\ Derek Orr, Jul 31 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Zak Seidov, Oct 20 2006
STATUS
approved