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

65537, 2, 5, 3, 7, 11, 61, 13, 120052642711, 157, 8837, 71, 3271, 31, 37, 19, 8597, 8868542017, 23, 149, 26443389827415776296437975011182059459845296976025400711, 97, 751, 17, 280121, 4999, 8797543, 12241, 29, 53, 197, 5659, 3109, 1063, 727915549

Offset: 1

Author

Labos Elemer

Keywords

• nonn

Links

• Tyler Busby, Table of n, a(n) for n = 1..70

Crossrefs

Cf. A020639, A000945, A000946, A005265, A005266.

Programs

Mathematica

a[1]=65537; 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 \\ A020639
first(m)=my(v=vector(m)); v[1]=65537; for(i=2, m, v[i]=spf(1+prod(j=1, i-1, v[j]))); v \\ Anders Hellström, Dec 07 2015

Extensions

a(25)-a(35) from Robert Price, Jul 09 2015