%I #20 Sep 08 2022 08:45:31
%S 13,43,73,103,163,193,223,283,313,373,433,463,523,613,643,673,733,823,
%T 853,883,1033,1063,1093,1123,1153,1213,1303,1423,1453,1483,1543,1663,
%U 1693,1723,1753,1783,1873,1933,1993,2053,2083,2113,2143,2203,2293,2383
%N Primes congruent to 13 (mod 30).
%C Primes ending in 3 with (SOD-1)/3 integer where SOD is sum of digits. - _Ki Punches_
%C Subsequence of primes of A082369. - _Michel Marcus_, Jan 23 2016
%H Vincenzo Librandi, <a href="/A132233/b132233.txt">Table of n, a(n) for n = 1..1000</a>
%H C. K. Caldwell, <a href="http://primes.utm.edu">The Prime Pages</a>.
%H Omar E. Pol, <a href="http://www.polprimos.com">Determinacion geometrica de los numeros primos y perfectos</a>.
%F a(n) = A158746(n)*30 + 13. - _Ray Chandler_, Apr 07 2009
%F Intersection of A030431 and A002476. - _Ray Chandler_, Apr 07 2009
%p select(isprime, [seq(30*i+13,i=0..1000)]); # _Robert Israel_, Jan 24 2016
%t Select[Prime[Range[1000]],MemberQ[{13},Mod[#,30]]&] (* _Vincenzo Librandi_, Aug 14 2012 *)
%o (Magma) [p: p in PrimesUpTo(3000) | p mod 30 eq 13 ]; // _Vincenzo Librandi_, Aug 14 2012
%o (PARI) lista(nn) = forprime(p=2, nn, if(p % 30 == 13, print1(p, ", "))); \\ _Altug Alkan_, Jan 23 2016
%Y Cf. A000040, A039949, A082369, A132230-A132236.
%K nonn,easy
%O 1,1
%A _Omar E. Pol_, Aug 15 2007
%E Extended by _Ray Chandler_, Apr 07 2009