|
%I #12 Sep 08 2022 08:45:35
%S 37,103,367,433,499,631,829,1093,1291,1423,1489,1621,1753,1951,2017,
%T 2083,2281,2347,2677,3271,3469,3733,3931,4129,4261,4327,4591,4657,
%U 4723,4789,4987,5119,5449,5581,5647,5779,6043,6373,6571,6637,6703,6967,7297,7561
%N Primes congruent to 4 mod 33.
%H Vincenzo Librandi, <a href="/A142051/b142051.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) ~ 20n log n. - _Charles R Greathouse IV_, Jul 03 2016
%t Select[Range[4, 20000, 33], PrimeQ] (* _Vladimir Joseph Stephan Orlovsky_, Jun 25 2011 *)
%t Select[Prime[Range[3000]],MemberQ[{4},Mod[#,33]]&] (* _Vincenzo Librandi_, Aug 18 2012 *)
%o (Magma) [p: p in PrimesUpTo(8000) | p mod 33 eq 4 ]; // _Vincenzo Librandi_, Aug 18 2012
%o (PARI) is(n)=isprime(n) && n%33==4 \\ _Charles R Greathouse IV_, Jul 03 2016
%Y Cf. A000040.
%K nonn,easy
%O 1,1
%A _N. J. A. Sloane_, Jul 11 2008
|
