OFFSET
1,1
COMMENTS
A prime of the form 16n+9 is represented by exactly one of x^2 + 32y^2 and x^2 + 64y^2 (see Kaplanski link). - Michel Marcus, Dec 23 2012
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
Irving Kaplansky, The forms x+32y^2 and x+64y^2, Proc. Amer. Math. Soc. 131 (2003), 2299-2300
MAPLE
M:=500; f:=proc(n) local t1, t2; t1:=[]; for k from 0 to M do t2:=2^n*k+2^(n-1)+1; if isprime(t2) then t1:=[op(t1), t2]; fi; od; t1; end; f(4);
MATHEMATICA
lst={}; Do[p=16*n+9; If[PrimeQ[p], AppendTo[lst, p]], {n, 0, 3*5!, 1}]; lst (* Vladimir Joseph Stephan Orlovsky, Jan 27 2009 *)
PROG
(PARI) select(n->n%16==9, primes(500)) \\ Charles R Greathouse IV, Apr 29 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, based on correspondence from Marco Matosic, Apr 11 2005
STATUS
approved