OFFSET
1,2
EXAMPLE
MAPLE
aklist := [] ;
A129645 := proc(n)
option remember;
global aklist ;
local prev;
if n = 1 then
return 1;
else
prev := procname(n-1) ;
sort([op(numtheory[divisors](prev))]) ;
aklist := [op(aklist), op(%)] ;
return prev*op(n-1, aklist)+1 ;
end if;
end proc:
seq(A129645(n), n=1..26) ; # R. J. Mathar, Aug 24 2013
MATHEMATICA
aklist = {};
a[n_] := a[n] = If[n == 1, 1, With[{prev = a[n-1]}, aklist = Join[aklist, Divisors[prev]]; prev*aklist[[n-1]]+1]];
Table[a[n], {n, 1, 26}] (* Jean-François Alcover, Jan 27 2024, after R. J. Mathar *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Apr 25 2007
EXTENSIONS
More terms from Sean A. Irvine, May 19 2010; corrected (based on the ordered lists of divisors) by Franklin T. Adams-Watters, Aug 24 2013
STATUS
approved