%I #10 Sep 08 2022 08:45:11

%S 233,313,569,953,1193,2473,3769,4793,6569,10169,11833,12713,16553,

%T 22073,28393,29753,32569,40169,52153,53993,61673,71993,78569,115769,

%U 129769,135593,144569,160169,169913,186793,222953,238313,250169,258233,262313

%N Primes of the form 16*m^2 + 169, m=1,2,3,...

%H Vincenzo Librandi, <a href="/A087862/b087862.txt">Table of n, a(n) for n = 1..10000</a>

%t Select[Table[16n^2+169,{n,0,700}],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..200] | IsPrime(a) where a is 16*n^2+169]; // _Vincenzo Librandi_, Dec 11 2011

%K nonn,easy

%O 1,1

%A _Cino Hilliard_, Oct 09 2003