STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
h[n_] := h[n] = DigitCount[n, 2, 1];
a[n_] := a[n] = If[n == 0, 1,
Sum[a[n - i]*h[i], {i, 1, n}]];
Table[a[n], {n, 0, 36}] (* Jean-François Alcover, Apr 14 2022, after Alois P. Heinz *)
approved
editing
editing
approved
Alois P. Heinz, <a href="/A341020/b341020.txt">Table of n, a(n) for n = 0..2984</a>
Number of compositions of n into parts where there are A000120(k) sorts of part k.
G.f.: 1 / (1 - Sum_{k>=1} A000120(k) * x^k).
0
INVERT transform of the binary weight.
h:= proc(n) option remember; add(j, j=Bits[Split](n)) end:
a:= proc(n) option remember; `if`(n=0, 1,
add(a(n-i)*h(i), i=1..n))
end:
seq(a(n), n=0..36);