OFFSET
1,1
COMMENTS
The initial primes 3 and 7 give essentially A000945.
EXAMPLE
5*2*11*3 + 1 = 331, which is prime; the least prime factor of 330*331 + 1 = 109231 = 19*5749 is 19, so a(6) = 19.
MATHEMATICA
a[1]=5; a[n_] := First[ Flatten[ FactorInteger[ 1+Product[ a[ j ], {j, 1, n-1} ] ] ] ]; Array[a, 10]
PROG
(PARI) spf(n)=my(f=factor(n)[1, 1]); f;
first(m)=my(v=vector(m)); v[1]=5; for(i=2, m, v[i]=spf(1+prod(j=1, i-1, v[j]))); v; \\ Anders Hellström, Aug 15 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
a(38)-a(39) from Robert Price, Jul 19 2015
STATUS
approved