|
%I #13 Apr 25 2022 08:09:33
%S 1,3,6,13,29,64,137,285,584,1197,2479,5198,10989,23275,49122,102951,
%T 213915,440478,899177,1821577,3668080,7358199,14742653,29582240,
%U 59593613,120742119,246254254,505494077,1043327165,2161750892,4488143133,9319949375,19326922734
%N Binomial transform of A002487.
%C The transform omits the initial two terms 0, 1 of A002487. - _Georg Fischer_, Oct 15 2021
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%p b:= proc(n) option remember; `if`(n<2, n,
%p (q-> b(q)+(n-2*q)*b(n-q))(iquo(n, 2)))
%p end:
%p a:= n-> add(b(n+2-j)*binomial(n-1, j-1), j=1..n):
%p seq(a(n), n=1..33); # _Alois P. Heinz_, Oct 15 2021
%t b[n_] := b[n] = If[n < 2, n,
%t With[{q = Quotient[n, 2]}, b[q] + (n-2*q)*b[n-q]]];
%t a[n_] := Sum[b[n+2-j]*Binomial[n-1, j-1], {j, 1, n}];
%t Table[a[n], {n, 1, 33}] (* _Jean-François Alcover_, Apr 25 2022, after _Alois P. Heinz_ *)
%K nonn
%O 1,2
%A _N. J. A. Sloane_, May 24 2002
|
