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

41, 2, 83, 3, 7, 47, 71, 29, 653, 5, 173, 23, 103058819, 11, 389, 73161901, 168593357, 13, 45613, 347, 211, 53, 400947612985987, 28837, 35111, 913011302795748880783905085999338914209329333652950191830525020998365540320068611

Offset

1,1

Links

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

Mathematica

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

Prog

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

Crossrefs

Cf. A000945, A000946, A005265, A005266.

Sequence in context: A053342 A297205 A283895 * A036189 A198190 A106424

Adjacent sequences: A051314 A051315 A051316 * A051318 A051319 A051320

Keyword

easy,nonn

Author

Labos Elemer

Status

approved