OFFSET
1,1
COMMENTS
A subsequence of the primes of the form 2k^2+4k+1 = 2*(k+1)^2-1, A066436. - R. J. Mathar, Aug 10 2009
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
lst={}; Do[p=Prime@n; a=2*p^2+4*p+1; If[PrimeQ@a, AppendTo[lst, a]], {n, 7!}]; lst (* Vladimir Joseph Stephan Orlovsky, Aug 12 2009 *)
Select[Table[2p^2+4p+1, {p, Prime[Range[200]]}], PrimeQ] (* Harvey P. Dale, Aug 25 2019 *)
PROG
(Magma) [a: p in PrimesUpTo(700)|IsPrime(a) where a is 2*p^2+4*p+1 ] // Vincenzo Librandi, Sep 01 2012
(PARI) lista(nn) = {forprime(p=2, nn, if(isprime(q=2*p^2+4*p+1), print1(q, ", "))); } \\ Altug Alkan, Mar 29 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Aug 08 2009
EXTENSIONS
a(29) corrected by R. J. Mathar, Aug 11 2009
Edited by N. J. A. Sloane, Aug 11 2009
STATUS
approved