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

83, 2, 167, 3, 7, 67, 5, 13, 719, 37, 11, 31, 1367, 31440216015620321911, 988487183108868589955299792587646370011, 19, 499, 937, 23, 29, 863, 1852812869, 1157393039341097158113740816072403984079, 445807, 13833892959043, 359, 109, 331

Offset: 1

Author

Labos Elemer

Keywords

• nonn

Links

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

Crossrefs

Cf. A000945, A000946, A005265, A005266.

Programs

Mathematica

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

PARI

gpf(n)=my(f=factor(n)[, 1]); f[#f];
first(m)=my(v=vector(m)); v[1]=83; for(i=2, m, v[i]=gpf(1+prod(j=1, i-1, v[j]))); v \\ Anders Hellström, Nov 25 2015

Extensions

More terms from Sean A. Irvine, Apr 12 2011