OFFSET
1,1
COMMENTS
Odd primes p such that -1 is a 32nd power mod p. - Eric M. Schmidt, Mar 27 2014
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
MATHEMATICA
lst={}; Do[p=Prime[n]; If[Mod[p, 64]==1, AppendTo[lst, p]], {n, 7!}]; lst (* Vladimir Joseph Stephan Orlovsky, Feb 26 2009 *)
Select[Prime[Range[1800]], MemberQ[{1}, Mod[#, 64]] &] (* Vincenzo Librandi, Sep 06 2012 *)
Select[Range[1, 12000, 64], PrimeQ] (* Harvey P. Dale, Dec 25 2023 *)
PROG
(Haskell)
a142925 n = a142925_list !! (n-1)
a142925_list = filter ((== 1) . a010051) [1, 65..]
-- Reinhard Zumkeller, Mar 06 2012
(Magma) [p: p in PrimesUpTo(12000) | p mod 64 eq 1 ] ; // Vincenzo Librandi, Sep 06 2012
(PARI) is(n)=isprime(n) && n%64==1 \\ Charles R Greathouse IV, Jul 01 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Jul 11 2008
STATUS
approved