%I #10 Sep 08 2022 08:45:11
%S 41,89,281,601,809,1049,1321,2729,4649,5209,5801,9241,11689,12569,
%T 15401,17449,18521,21929,26921,41641,46681,52009,53849,55721,59561,
%U 71849,80681,87641,94889,97369,99881,107609,112921,118361,126761,135449
%N Primes of the form 16*m^2 + 25 for m=1,2,3,...
%H Vincenzo Librandi, <a href="/A087857/b087857.txt">Table of n, a(n) for n = 1..5000</a>
%t Select[Table[16n^2+25,{n,0,4000}],PrimeQ] (* _Vincenzo Librandi_, Dec 11 2011 *)
%o (PARI) fourmp1(m,n) = { forstep(x=1,m,1, y=16*(x^2+n^2)+8*n+1; if(isprime(y),print1(y",")) ) }
%o (Magma) [a: n in [0..100] | IsPrime(a) where a is 16*n^2+25]; // _Vincenzo Librandi_, Dec 11 2011
%K nonn,easy
%O 1,1
%A _Cino Hilliard_, Oct 09 2003
|