OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..446
FORMULA
a(n) = n! * Sum_{k=0..floor(n/3)} binomial(n-2*k,n-3*k)/(6^k * k!).
D-finite with recurrence: 2*(n + 1)*(n + 2)*(n + 3)*a(n) - 3*(n + 2)*(n + 3)*a(n+1) + 6*(n + 3)^2*a(n + 2) - 6*(7 + 2*n)*a(n + 3) + 6*a(n + 4). - Robert Israel, May 06 2026
MAPLE
f:= gfun:-rectoproc({{2*a(n)*(n + 1)*(n + 2)*(n + 3) - 3*a(n + 1)*(n + 2)*(n + 3) + 6*(n + 3)^2*a(n + 2) - 6*(7 + 2*n)*a(n + 3) + 6*a(n + 4), a(0) = 1, a(1) = 1, a(2) = 2, a(3) = 7}, a(n), remember):
map(f, [$0..30]); # Robert Israel, May 06 2026
PROG
(PARI) a(n) = n!*sum(k=0, n\3, binomial(n-2*k, n-3*k)/(6^k*k!));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jun 18 2024
STATUS
approved
