A360026 a(n) = Sum_{k=0..floor(n/4)} (-1)^k * binomial(n-3*k,k) * Catalan(k).

1, 1, 1, 1, 0, -1, -2, -3, -2, 1, 6, 13, 17, 13, -4, -39, -83, -113, -92, 31, 279, 605, 850, 701, -219, -2129, -4736, -6749, -5690, 1569, 17114, 38713, 55957, 48249, -11498, -142163, -326860, -478957, -421262, 84015, 1210831, 2829363, 4197670, 3762583, -601732

Offset

0,7

Links

Seiichi Manyama, Table of n, a(n) for n = 0..1000

Formula

a(n) = 1 - Sum_{k=0..n-4} a(k) * a(n-k-4).

G.f. A(x) satisfies: A(x) = 1/(1-x) - x^4 * A(x)^2.

G.f.: 2 / ( 1-x + sqrt((1-x)^2 + 4*x^4*(1-x)) ).

D-finite with recurrence +(n+4)*a(n) +2*(-n-3)*a(n-1) +(n+2)*a(n-2) +4*(n-2)*a(n-4) +2*(-2*n+5)*a(n-5)=0. - R. J. Mathar, Jan 25 2023

Prog

(PARI) a(n) = sum(k=0, n\4, (-1)^k*binomial(n-3*k, k)*binomial(2*k, k)/(k+1));
(PARI) my(N=50, x='x+O('x^N)); Vec(2/(1-x+sqrt((1-x)^2+4*x^4*(1-x))))

Crossrefs

Cf. A360024, A360025, A360027.

Cf. A000108, A346073, A349048.

Sequence in context: A077870 A294408 A358605 * A062323 A291087 A020858

Adjacent sequences: A360023 A360024 A360025 * A360027 A360028 A360029

Keyword

sign

Author

Seiichi Manyama, Jan 22 2023

Status

approved