login
A051322
Euclid-Mullin sequence (A000945) with initial value a(1)=61 instead of a(1)=2.
1
61, 2, 3, 367, 7, 5, 1217, 17, 13, 59, 3271, 19, 11, 76938833, 337, 10711, 1021, 31, 83, 1290503, 15608141, 195648590992627, 109, 15983112328343713807564263439978033717550587578630760604057976854157331
OFFSET
1,1
LINKS
Tyler Busby, Table of n, a(n) for n = 1..42 (terms 1..33 from Robert Price)
MATHEMATICA
a[1]=61; 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]=61; for(i=2, m, v[i]=spf(1+prod(j=1, i-1, v[j]))); v \\ Anders Hellström, Dec 04 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved