OFFSET
0,2
LINKS
FORMULA
a(n) = (1/(n+1)) * Sum_{k=0..floor(n/2)} (-1)^k * binomial(n+k,k) * binomial(4*n+k+2,n-2*k).
D-finite with recurrence: (-2048*n^3 - 7168*n^2 - 8064*n - 2880)*a(n) + (16640*n^3 + 86400*n^2 + 149328*n + 86040)*a(n + 1) + (-13440*n^3 - 111360*n^2 - 308668*n - 286644)*a(n + 2) + (2648*n^3 + 27184*n^2 + 92312*n + 103776)*a(n + 3) + (-155*n^3 - 1860*n^2 - 7285*n - 9300)*a(n + 4) = 0. - Robert Israel, Mar 11 2026
MAPLE
f:= gfun:-rectoproc({(-2048*n^3 - 7168*n^2 - 8064*n - 2880)*a(n) + (16640*n^3 + 86400*n^2 + 149328*n + 86040)*a(n + 1) + (-13440*n^3 - 111360*n^2 - 308668*n - 286644)*a(n + 2) + (2648*n^3 + 27184*n^2 + 92312*n + 103776)*a(n + 3) + (-155*n^3 - 1860*n^2 - 7285*n - 9300)*a(n + 4), a(0) = 1, a(1) = 3, a(2) = 14, a(3) = 76}, a(n), remember):
map(f, [$0..30]); # Robert Israel, Mar 11 2026
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serreverse(x*((1-x)^3+x^2))/x)
(PARI) a(n) = sum(k=0, n\2, (-1)^k*binomial(n+k, k)*binomial(4*n+k+2, n-2*k))/(n+1);
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Mar 23 2024
STATUS
approved
