A368575 a(n) = n! * Sum_{k=0..n} binomial(k+3,4) / k!.

0, 1, 7, 36, 179, 965, 5916, 41622, 333306, 3000249, 30003205, 330036256, 3960436437, 51485675501, 720799459394, 10811991893970, 172991870307396, 2940861795230577, 52935512314156371, 1005774733968978364, 20115494679379576135, 422425388266971109461

Offset

0,3

Links

Table of n, a(n) for n=0..21.

Formula

a(0) = 0; a(n) = n*a(n-1) + binomial(n+3,4).

E.g.f.: x * (1+3*x/2+x^2/2+x^3/24) * exp(x) / (1-x).

Prog

(PARI) my(N=30, x='x+O('x^N)); concat(0, Vec(serlaplace(x*sum(k=0, 3, binomial(3, k)*x^k/(k+1)!)*exp(x)/(1-x))))

Crossrefs

Cf. A007526, A103519, A368574, A368576.

Cf. A000332, A337002.

Sequence in context: A246417 A110310 A054493 * A037538 A037482 A147546

Adjacent sequences: A368572 A368573 A368574 * A368576 A368577 A368578

Keyword

nonn,easy

Author

Seiichi Manyama, Dec 31 2023

Status

approved