|
%I #19 Jan 24 2023 06:56:51
%S 13,17,43,47,73,103,107,137,163,167,193,197,223,227,257,283,313,317,
%T 347,373,433,463,467,523,557,587,613,617,643,647,673,677,733,797,823,
%U 827,853,857,883,887,947,977,1033,1063,1093,1097,1123,1153,1187,1213,1217
%N Primes congruent to 13 or 17 mod 30.
%H Vincenzo Librandi, <a href="/A140533/b140533.txt">Table of n, a(n) for n = 1..1000</a>
%p A140533 := proc(n) local a; if n = 1 then 13; else a := nextprime(procname(n-1)) ; while not a mod 30 in {13,17} do a := nextprime(a) ; end do: return a; end if; end: seq(A140533(n),n=1..80) ; # _R. J. Mathar_, Oct 22 2009
%t Select[Prime[Range[300]],MemberQ[{13,17},Mod[#,30]]&] (* _Vincenzo Librandi_, Aug 15 2012 *)
%o (Magma) [p: p in PrimesUpTo(3000)|p mod 30 in {13,17}]; // _Vincenzo Librandi_, Dec 18 2010
%Y Cf. A132233, A039949.
%Y Cf. A000040.
%K nonn,easy
%O 1,1
%A _Juri-Stepan Gerasimov_, Jul 28 2008
|
