Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
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
OFFSET
1,1
LINKS
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
KEYWORD
nonn,easy
AUTHOR
Cino Hilliard, Oct 09 2003
STATUS
approved