A087857 Primes of the form 16*m^2 + 25 for m=1,2,3,...

41, 89, 281, 601, 809, 1049, 1321, 2729, 4649, 5209, 5801, 9241, 11689, 12569, 15401, 17449, 18521, 21929, 26921, 41641, 46681, 52009, 53849, 55721, 59561, 71849, 80681, 87641, 94889, 97369, 99881, 107609, 112921, 118361, 126761, 135449

Offset

1,1

Links

Vincenzo Librandi, Table of n, a(n) for n = 1..5000

Mathematica

Select[Table[16n^2+25, {n, 0, 4000}], PrimeQ] (* Vincenzo Librandi, Dec 11 2011 *)

Prog

(PARI) fourmp1(m, n) = { forstep(x=1, m, 1, y=16*(x^2+n^2)+8*n+1; if(isprime(y), print1(y", ")) ) }
(Magma) [a: n in [0..100] | IsPrime(a) where a is 16*n^2+25]; // Vincenzo Librandi, Dec 11 2011

Crossrefs

Sequence in context: A155572 A356246 A107145 * A139995 A044179 A044560

Adjacent sequences: A087854 A087855 A087856 * A087858 A087859 A087860

Keyword

nonn,easy

Author

Cino Hilliard, Oct 09 2003

Status

approved