OFFSET
1,1
COMMENTS
Apparently a(n) = A025584(n) for all n>2. Verified for indices n <= 100000. - R. J. Mathar, Mar 01 2010
Assume a(n-1) is in A025584 and greater than 3. Then, if a(n) = a(n-1) + 4, a(n-1) + 2 cannot be prime, and a(n) is in A025584. Otherwise, a(n) must still be in A025584, because if it were not, a(n) would not be the smallest valid prime. - Charlie Neder, Dec 16 2018
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..10000
MAPLE
p := 2: for n from 1 to 100 do printf("%d, ", p) ; p := nextprime(p+3) ; od: # R. J. Mathar, Feb 21 2009
MATHEMATICA
a[1] = 2; a[n_] := a[n] = NextPrime[a[n - 1] + 3]; Array[a, 60] (* Amiram Eldar, Dec 16 2018 *)
NestList[NextPrime[#+3]&, 2, 60] (* Harvey P. Dale, Jan 05 2019 *)
PROG
(Magma) [n le 1 select 2 else NextPrime(Self(n-1)+3): n in [1..60]]; // Vincenzo Librandi, Dec 21 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Feb 16 2009
EXTENSIONS
More terms from R. J. Mathar, Feb 21 2009
STATUS
approved