OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..386
FORMULA
a(n) = [x^n] 1/(1+Sum_{j=1..n} (n+1-j)*(-x)^j).
EXAMPLE
a(3) = 16: 000, 002, 020, 021, 022, 100, 102, 110, 111, 200, 202, 210, 211, 220, 221, 222 (using ternary alphabet {0, 1, 2}).
MAPLE
b:= proc(n, k) option remember; `if`(n<0, 0, `if`(n=0, 1,
-add((-1)^j*(k+1-j)*b(n-j, k), j=1..k)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..25);
MATHEMATICA
b[n_, k_] := b[n, k] = If[n < 0, 0, If[n == 0, 1,
-Sum[(-1)^j (k+1-j) b[n-j, k], {j, 1, k}]]];
a[n_] := b[n, n];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Nov 02 2021, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 26 2016
STATUS
approved