A087872 Primes of the form (4*k + 3)^2 + (4*k + 2)^2 where k=0,1,2,3,...

13, 421, 1013, 1861, 2381, 3613, 5101, 9941, 13613, 16381, 21013, 24421, 30013, 34061, 47741, 52813, 60901, 82013, 99013, 106261, 168781, 218461, 223781, 234613, 304981, 337021, 370661, 377581, 391613, 405901, 413141, 427813, 450301, 481181

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

Cf. A088759.

Sequence in context: A258178 A266486 A142484 * A308341 A098890 A012023

Adjacent sequences: A087869 A087870 A087871 * A087873 A087874 A087875

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