OFFSET
0,2
LINKS
FORMULA
a(n) = (1/(n+1)) * Sum_{k=0..n} (-1)^k * binomial(n+k,k) * binomial(4*n+2*k+2,n-k).
D-finite with recurrence: (-768*n^3 - 2688*n^2 - 3024*n - 1080)*a(n) + (3073*n^3 + 17925*n^2 + 35318*n + 23400)*a(n + 1) + (-3374*n^3 - 28761*n^2 - 81889*n - 77916)*a(n + 2) + (925*n^3 + 9540*n^2 + 32555*n + 36780)*a(n + 3) + (-92*n^3 - 1104*n^2 - 4324*n - 5520)*a(n + 4) = 0. - Robert Israel, Mar 11 2026
MAPLE
f:= gfun:-rectoproc({(-768*n^3 - 2688*n^2 - 3024*n - 1080)*a(n) + (3073*n^3 + 17925*n^2 + 35318*n + 23400)*a(n + 1) + (-3374*n^3 - 28761*n^2 - 81889*n - 77916)*a(n + 2) + (925*n^3 + 9540*n^2 + 32555*n + 36780)*a(n + 3) + (-92*n^3 - 1104*n^2 - 4324*n - 5520)*a(n + 4), a(0) = 1, a(1) = 2, a(2) = 5, a(3) = 11}, 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))/x)
(PARI) a(n) = sum(k=0, n, (-1)^k*binomial(n+k, k)*binomial(4*n+2*k+2, n-k))/(n+1);
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Mar 23 2024
STATUS
approved
