OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MAPLE
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
MATHEMATICA
Select[Prime[Range[300]], MemberQ[{13, 17}, Mod[#, 30]]&] (* Vincenzo Librandi, Aug 15 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(3000)|p mod 30 in {13, 17}]; // Vincenzo Librandi, Dec 18 2010
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Juri-Stepan Gerasimov, Jul 28 2008
STATUS
approved