STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
b[n_] := b[n] = If[n < 2, n,
With[{q = Quotient[n, 2]}, b[q] + (n-2*q)*b[n-q]]];
a[n_] := Sum[b[n+2-j]*Binomial[n-1, j-1], {j, 1, n}];
Table[a[n], {n, 1, 33}] (* Jean-François Alcover, Apr 25 2022, after Alois P. Heinz *)
approved
editing
proposed
approved
editing
proposed
b:= proc(n) option remember; `if`(n<2, n,
(q-> b(q)+(n-2*q)*b(n-q))(iquo(n, 2)))
end:
a:= n-> add(b(n+2-j)*binomial(n-1, j-1), j=1..n):
seq(a(n), n=1..33); # Alois P. Heinz, Oct 15 2021
proposed
editing
editing
proposed
The transform omits the initial two terms 0, 1 of A002487. - Georg Fischer, Oct 15 2021
approved
editing
_N. J. A. Sloane (njas(AT)research.att.com), _, May 24 2002