OFFSET
1,1
COMMENTS
Subsequence of A007519 (primes of form 8n+1). - Zak Seidov, May 16 2012
Primes p such that p XOR 14 = p + 14. - Brad Clardy, Jul 23 2012
A prime of the form 16n+1 is represented either by both x^2+32y^2 and x^2+64y^2 or by neither (see Kaplansky link). - Michel Marcus, Dec 23 2012
Odd primes p such that -1 is an 8th power mod p. - Eric M. Schmidt, Mar 27 2014
LINKS
T. D. Noe, Table of n, a(n) for n=1..1000
C, Caldwell, Prime test.
Irving Kaplansky, The forms x+32y^2 and x+64y^2, Proc. Amer. Math. Soc. 131 (2003), 2299-2300
MAPLE
p:=proc(n) if isprime(16*n+1)=true then 16*n+1 else fi end:seq(p(n), n=1..200); # Emeric Deutsch, Dec 23 2004
MATHEMATICA
lst={}; Do[p=16*n+1; If[PrimeQ[p], AppendTo[lst, p]], {n, 6!}]; lst (* Vladimir Joseph Stephan Orlovsky, Feb 26 2009 *)
Select[16*Range[200]+1, PrimeQ] (* Harvey P. Dale, Nov 04 2017 *)
PROG
(Haskell)
a094407 n = a094407_list !! (n-1)
a094407_list = filter ((== 1) . a010051) [1, 17..]
-- Reinhard Zumkeller, Mar 06 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jun Mizuki (suzuki32(AT)sanken.osaka-u.ac.jp), Jun 03 2004
EXTENSIONS
More terms from Emeric Deutsch, Dec 23 2004
STATUS
approved