OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
lst={}; Do[p=Prime[n]; If[Mod[p, 32]==31, AppendTo[lst, p]], {n, 7!}]; lst (* Vladimir Joseph Stephan Orlovsky, Feb 26 2009 *)
Select[Prime[Range[6000]], Mod[#, 32]==31&] (* Vincenzo Librandi, Aug 14 2012 *)
Select[Range[31, 5300, 32], PrimeQ] (* Harvey P. Dale, Oct 03 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(6000) | p mod 32 eq 31 ]; // Vincenzo Librandi, Aug 14 2012
(PARI) is(n)=isprime(n) && n%32==31 \\ Charles R Greathouse IV, Jul 01 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Artur Jasinski, Jan 19 2007
EXTENSIONS
Corrected by N. J. A. Sloane, Jul 11 2008
Edited by N. J. A. Sloane at the suggestion of R. J. Mathar, Jul 20 2008
STATUS
approved