OFFSET
1,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..389
FORMULA
a(n) = A222029(n,n).
MAPLE
b:= proc(n, m) option remember; `if`(n=0, x^m, add((j-1)!*
b(n-j, ilcm(m, j))*binomial(n-1, j-1), j=1..n))
end:
a:= n-> add(coeff(b(j, 1), x, n)*n^(n-j)*binomial(n-1, j-1), j=0..n):
seq(a(n), n=1..25);
MATHEMATICA
b[n_, m_] := b[n, m] = If[n == 0, x^m, Sum[(j - 1)!*
b[n - j, LCM[m, j]]*Binomial[n - 1, j - 1], {j, 1, n}]];
a[n_] := Sum[Coefficient[b[j, 1], x, n]*n^(n-j)*Binomial[n-1, j-1], {j, 0, n}];
Table[a[n], {n, 1, 25}] (* Jean-François Alcover, Mar 07 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 15 2017
STATUS
approved