

A087857


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


2



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
(list;
graph;
refs;
listen;
history;
text;
internal format)



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: A139924 A155572 A107145 * A139995 A044179 A044560
Adjacent sequences: A087854 A087855 A087856 * A087858 A087859 A087860


KEYWORD

nonn,easy


AUTHOR

Cino Hilliard, Oct 09 2003


STATUS

approved



