OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..150
MAPLE
b:= proc(n, m) option remember; `if`(n=0, m, add((j-1)!*
b(n-j, ilcm(m, j))*binomial(n-1, j-1), j=1..n))
end:
a:= n-> add(b(j, 1)*n^(n-j)*binomial(n-1, j-1), j=0..n):
seq(a(n), n=0..25);
MATHEMATICA
T[n_, k_] := T[n, k] = If[n == 0, k, Sum[(j - 1)! * T[n - j, LCM[k, j]]*Binomial[n - 1, j - 1], {j, n}]]; {1}~Join~Table[Sum[T[j, 1]*n^(n - j)*Binomial[n - 1, j - 1], {j, 0, n}], {n, 19}] (* Michael De Vlieger, Aug 17 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 13 2017
STATUS
approved