OFFSET
1,1
COMMENTS
Primes of the form 32*k^2 + 40*k + 13. - Vincenzo Librandi, Dec 12 2011
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..5000
MATHEMATICA
Select[Table[32n^2+40n+13, {n, 0, 800}], PrimeQ] (* Vincenzo Librandi, Dec 12 2011 *)
PROG
(PARI) fourmp1b(n) = { for(x=1, n, z=(4*x+3)^2+(4*x+2)^2; if(isprime(z), print1(z", ")) ) }
(Magma) [a: n in [0..200] | IsPrime(a) where a is 32*n^2+40*n+13]; // Vincenzo Librandi, Dec 12 2011
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Cino Hilliard, Oct 10 2003
EXTENSIONS
More terms from Ray Chandler, Oct 19 2003
a(1)=13 inserted by Vincenzo Librandi, Dec 12 2011
Description revised to agree with insertion of new term by Ray Chandler, May 21 2017
STATUS
approved