A051321 Euclid-Mullin sequence (A000945) with initial value a(1)=59 instead of a(1)=2.

59, 2, 7, 827, 3, 13, 4583, 5, 610478010871, 61, 292658543, 4483, 47, 11, 31, 43, 16453, 41, 3671, 1982639, 628319, 841476613, 449, 997793, 73, 983, 53, 28475917, 19, 673

Offset: 1

Author

Labos Elemer

Keywords

• easy
• nonn

Links

• Robert Price, Table of n, a(n) for n = 1..48

Crossrefs

Cf. A000945, A000946, A005265, A005266.

Programs

Mathematica

a[1]=59; a[n_] := First[ Flatten[ FactorInteger[ 1+Product[ a[ j ], {j, 1, n-1} ] ] ] ]; Array[a, 10]

PARI

spf(n)=my(f=factor(n)[1, 1]); f
first(m)=my(v=vector(m)); v[1]=59; for(i=2, m, v[i]=spf(1+prod(j=1, i-1, v[j]))); v \\ Anders Hellström, Dec 04 2015

Extensions

a(28)-a(48) from Robert Price, Jul 13 2015