OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
EXAMPLE
a(4) = 1 + 2 + 2 + 2 + 1 + 2 + 2 + 1 = 13: (1)111, (1)1(2), (1)(2)1, (2)(1)1, (2)2, (1)(3), (3)(1), (4).
MAPLE
b:= proc(n, i, p) option remember; `if`(n=0, [p!, 0],
`if`(i<1, 0, add((p-> [0, `if`(j=0, 0, p[1])]+p)(
b(n-i*j, i-1, p+j)/j!), j=0..n/i)))
end:
a:= n-> b(n$2, 0)[2]:
seq(a(n), n=0..38);
MATHEMATICA
b[n_, i_, p_] := b[n, i, p] = If[n == 0, {p!, 0},
If[i<1, {0, 0}, Sum[{0, If[j == 0, 0, #[[1]]]}+#&[
b[n-i*j, i-1, p+j]/j!], {j, 0, n/i}]]];
a[n_] := b[n, n, 0][[2]];
a /@ Range[0, 38] (* Jean-François Alcover, Jun 13 2021, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 06 2020
STATUS
approved