A362166 Expansion of e.g.f. exp(-x * (1-3*x)^(1/3)).

1, -1, 3, -1, 41, 299, 4531, 74507, 1474481, 33540119, 864507491, 24891022199, 791755864153, 27571976573699, 1043247441846611, 42615848603499779, 1869129393654945761, 87605345727468933167, 4369604246576366377411, 231091472431638655755119

Offset

0,3

Links

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

Formula

a(n) = (-1)^n * n! * Sum_{k=0..n} 3^k * binomial((n-k)/3,k)/(n-k)!.

D-finite with recurrence +(-9*n+11)*a(n) +2*(27*n^2-121*n+72)*a(n-1) +3*(-27*n^3+304*n^2-1053*n+1056)*a(n-2) +(-612*n^3+6984*n^2-23677*n+21227) *a(n-3) +4*(27*n-23)*(n-3)*a(n-4) -48*(9*n-10) *(n-3)*(n-4) *a(n-5) +64*(n-5)*(n-4)*(9*n^2-62*n+78)*a(n-6) +256*(n-5) *(n-6)*(17*n-24)*(n-4)*a(n-7)=0. - R. J. Mathar, Dec 04 2023

Maple

A362166 := proc(n)
    (-1)^n*n!*add(3^k * binomial((n-k)/3, k)/(n-k)!, k=0..n) ;
end proc:
seq(A362166(n), n=0..70) ; # R. J. Mathar, Dec 04 2023

Prog

(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(-x*(1-3*x)^(1/3))))

Crossrefs

Cf. A362162, A362164.

Sequence in context: A125082 A307803 A356819 * A136517 A366479 A104097

Adjacent sequences: A362163 A362164 A362165 * A362167 A362168 A362169

Keyword

sign,easy

Author

Seiichi Manyama, Apr 10 2023

Status

approved