STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
E.g.f.: (exp(x) - Sum_{k=0..3} x^k/k!) * exp(exp(x) - 1). - Ilya Gutkovskiy, Jun 25 2022
approved
editing
editing
approved
# second Maple program:
b:= proc(n) option remember; `if`(n=0, [1, 0], add((p-> p+[0,
`if`(j>3, p[1], 0)])(b(n-j)*binomial(n-1, j-1)), j=1..n))
end:
a:= n-> b(n)[2]:
seq(a(n), n=4..30); # Alois P. Heinz, Jan 06 2022
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
b[n_] := b[n] = If[n == 0, 1, Sum[b[n - j]*Binomial[n-1, j-1], {j, 1, n}]];
g[n_, k_] := g[n, k] = If[n < k, 0, g[n, k+1] + Binomial[n, k]*b[n - k]];
a[n_] := g[n, 4];
Table[a[n], {n, 4, 30}] (* Jean-François Alcover, May 28 2018, from Maple *)
approved
editing