OFFSET
1,1
COMMENTS
Subsequence of A062324.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
For prime p = 3, p^2+4 = 13 and p^2+4p+2 = 23 are prime; for p = 67, p^2+4 = 4493 and p^2+4p+2 = 4759 are prime.
MAPLE
a := proc (n) if isprime(n) = true and isprime(n^2+4) = true and isprime(n^2+4*n+2) = true then n else end if end proc: seq(a(n), n = 1 .. 7000); # Emeric Deutsch, Jan 02 2009
MATHEMATICA
Select[Prime[Range[10000]], PrimeQ[#^2+4]&&PrimeQ[#^2 +4#+2]&] (* Vincenzo Librandi, Jul 27 2012 *)
PROG
(Magma) [ p: p in PrimesUpTo(7000) | IsPrime(p^2+4) and IsPrime(p^2+4*p+2) ];
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Dec 30 2008
EXTENSIONS
Edited, corrected (three terms deleted) and extended beyond a(10) by Klaus Brockhaus, Jan 02 2009
Corrected and extended by Emeric Deutsch, Jan 02 2009
STATUS
approved