OFFSET
1,4
MAPLE
a[1] := 1: for n from 2 to 100 do a[n] := n*a[n-1]/ifactors(n)[2][1][1]: od:seq(a[j], j=1..100);
MATHEMATICA
nxt[{n_, a_}]:={n+1, (a(n+1))/FactorInteger[n+1][[1, 1]]}; Transpose[ NestList[ nxt, {1, 1}, 30]][[2]] (* Harvey P. Dale, Jan 09 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Oct 18 2002
EXTENSIONS
More terms from Sascha Kurz, Jan 21 2003
STATUS
approved