OFFSET
18,1
LINKS
Joerg Arndt and Alois P. Heinz, Table of n, a(n) for n = 18..1000
EXAMPLE
a(18) = 13: [3,2,1,3,2,1,3,2,1], [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,3,2,1,2,1,3,2,1], [3,2,1,2,1,2,1,3,2,1], [3,2,1,2,1,3,2,1,2,1], [2,1,2,1,2,1,3,2,1,2,1], [2,1,2,1,3,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,3,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, 7), polynom))
end:
a:= n-> coeff(b(n, 0), x, 6):
seq(a(n), n=18..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, 6];
Table[a[n], {n, 18, 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