%I A087857
%S A087857 41,89,281,601,809,1049,1321,2729,4649,5209,5801,9241,11689,12569,15401,
%T A087857 17449,18521,21929,26921,41641,46681,52009,53849,55721,59561,71849,
%U A087857 80681,87641,94889,97369,99881,107609,112921,118361,126761,135449
%N A087857 Primes of the form 16m^2+25, for m=1,2,3...
%o A087857 (PARI) fourmp1(m,n) = { forstep(x=1,m,1, y=16*(x^2+n^2)+8*n+1; if(isprime(y),
               print1(y",")) ) }
%Y A087857 Sequence in context: A139924 A155572 A107145 this_sequence A139995 A044179 
               A044560
%Y A087857 Adjacent sequences: A087854 A087855 A087856 this_sequence A087858 A087859 
               A087860
%K A087857 nonn
%O A087857 1,1
%A A087857 Cino Hilliard (hillcino368(AT)gmail.com), Oct 09 2003