OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
FORMULA
a(n) = A327673(n,n).
EXAMPLE
a(3) = 19: 3abc, 3acb, 3bac, 3bca, 3cab, 3cba, 2ab1c, 2ac1b, 2ba1c, 2bc1a, 2ca1b, 2cb1a, 1a2bc, 1a2cb, 1b2ac, 1b2ca, 1c2ab, 1c2ba, 1a1b1c.
MAPLE
b:= proc(n, i, k, p) option remember;
`if`(n=0, p!, `if`(i<1, 0, add(binomial(k^i, j)*
b(n-i*j, min(n-i*j, i-1), k, p+j)/j!, j=0..n/i)))
end:
a:= n-> add(b(n$2, i, 0)*(-1)^(n-i)*binomial(n, i), i=0..n):
seq(a(n), n=0..21);
MATHEMATICA
b[n_, i_, k_, p_] := b[n, i, k, p] =
If[n == 0, p!, If[i < 1, 0, Sum[Binomial[k^i, j]*
b[n - i j, Min[n - i j, i - 1], k, p + j]/j!, {j, 0, n/i}]]];
a[n_] := Sum[b[n, n, i, 0] (-1)^(n-i) Binomial[n, i], {i, 0, n}];
Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Aug 01 2021, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Sep 21 2019
STATUS
approved