A361558 Expansion of e.g.f. exp((x + x^2/2 + x^3/6 + x^4/24)/(1-x)).

1, 1, 4, 20, 127, 976, 8776, 90084, 1037555, 13233077, 184956386, 2809098986, 46038214729, 809411443790, 15189361799522, 302932433571356, 6396529241755881, 142523960797017589, 3341115707515530400, 82187749261419720712, 2116421112495023612311

Offset

0,3

Links

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

Formula

a(n) = (2*n-1) * a(n-1) - (n-1)*(n-3) * a(n-2) - binomial(n-1,3) * a(n-4) - 3*binomial(n-1,4) * a(n-5) for n > 4.

Prog

(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(exp((x+x^2/2+x^3/6+x^4/24)/(1-x))))

Crossrefs

Cf. A361532, A361548, A361557.

Sequence in context: A266490 A135886 A007550 * A361557 A355167 A080795

Adjacent sequences: A361555 A361556 A361557 * A361559 A361560 A361561

Keyword

nonn

Author

Seiichi Manyama, Mar 15 2023

Status

approved