OFFSET
1,2
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
a(9)=31 because 9*a(8) + 1 = 9*179 + 1 = 1612 = 2*2*13*31.
MAPLE
with(numtheory): a:=proc(n) local f, L: if n = 1 then 1 else f:=factorset(1+n*a(n-1)): L:=convert(f, list): L[nops(L)] end if end proc: seq(a(n), n=1..35); # Emeric Deutsch, Nov 12 2007
MATHEMATICA
a = {1}; Do[AppendTo[a, FactorInteger[(Length[a] + 1)*a[[ -1]] + 1][[ -1, 1]]], {70}]; a (* Stefan Steinerberger, Oct 30 2007 *)
nxt[{a_, n_}]:={FactorInteger[a(n+1)+1][[-1, 1]], n+1}; Transpose[ NestList[ nxt, {1, 1}, 50]][[1]] (* Harvey P. Dale, Jun 02 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Oct 28 2007
EXTENSIONS
More terms from Stefan Steinerberger, Oct 30 2007
STATUS
approved