|
|
EXAMPLE
|
a(15) = 8: [3,2,1,2,1,3,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,3,2,1,2,1,2,1], [3,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,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, 6), polynom))
end:
a:= n-> coeff(b(n, 0), x, 5):
seq(a(n), n=15..50);
|
|
|
MATHEMATICA
|
k = 5;
b[n_, i_] := b[n, i] = If[n == 0, 1, Sum[b[n - j, j]*
If[j < i, x, 1], {j, n}] + O[x]^(k + 1)];
a[n_] := SeriesCoefficient[b[n, 0], {x, 0, k}];
a /@ Range[15, 50] (* Jean-François Alcover, Aug 27 2021, after Maple code *)
|
