OFFSET
21,1
LINKS
Joerg Arndt and Alois P. Heinz, Table of n, a(n) for n = 21..1000
EXAMPLE
a(21) = 21: [3,2,1,2,1,3,2,1,3,2,1], [3,2,1,3,2,1,2,1,3,2,1], [3,2,1,3,2,1,3,2,1,2,1], [2,1,3,2,1,3,2,1,3,2,1], [2,1,2,1,2,1,3,2,1,3,2,1], [3,2,1,3,2,1,2,1,2,1,2,1], [2,1,2,1,3,2,1,2,1,3,2,1], [2,1,3,2,1,2,1,3,2,1,2,1], [3,2,1,2,1,3,2,1,2,1,2,1], [3,2,1,2,1,2,1,2,1,3,2,1], [2,1,3,2,1,2,1,2,1,3,2,1], [3,2,1,2,1,2,1,3,2,1,2,1], [2,1,2,1,3,2,1,3,2,1,2,1], [2,1,3,2,1,3,2,1,2,1,2,1], [3,2,1,2,1,2,1,2,1,2,1,2,1], [2,1,2,1,2,1,2,1,2,1,3,2,1], [2,1,3,2,1,2,1,2,1,2,1,2,1], [2,1,2,1,2,1,2,1,3,2,1,2,1], [2,1,2,1,2,1,3,2,1,2,1,2,1], [2,1,2,1,3,2,1,2,1,2,1,2,1], [2,1,2,1,2,1,2,1,2,1,2,1,2,1].
MAPLE
b:= proc(n, i) option remember;
`if`(n=0, 1, convert(series(add(b(n-j, j)*
`if`(j<i, x, 1), j=1..n), x, 8), polynom))
end:
a:= n-> coeff(b(n, 0), x, 7):
seq(a(n), n=21..50);
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0, 1, Sum[b[n - j, j]*
If[j < i, x, 1], {j, 1, n}] // Expand];
a[n_] := Coefficient[b[n, 0], x, 7];
Table[a[n], {n, 21, 50}] (* Jean-François Alcover, Nov 28 2023, from Maple code *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Joerg Arndt and Alois P. Heinz, Apr 26 2014
STATUS
approved