A368709 a(n) = hypergeom([-1 - n, -n, 1 - n], [2, 3], +2).

1, 1, -1, -3, 13, 17, -241, 121, 5081, -13327, -106705, 609589, 1850661, -23392159, -6796193, 811545073, -1688514383, -25224774367, 123764707231, 650087614573, -6385330335427, -9591188592399, 279171512779759, -318526766092183, -10665705513959287, 40625771132796817

Offset

0,4

Links

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

Formula

a(n) = (-1/2)*B(n, -2) where B(n, x) are the Baxter polynomials with coefficients A359363, for n > 0. - Peter Luschny, Jan 04 2024

a(0) = 1, a(n) = (-1)^n*2^(n + 1)/(n*(n + 1)^2)*Sum_{k=1..n} (-1/2)^k*binomial(n + 1, k - 1)*binomial(n + 1, k)*binomial(n + 1, k + 1). - Detlef Meya, May 29 2024

Mathematica

Table[HypergeometricPFQ[{-1-n, -n, 1-n}, {2, 3}, 2], {n, 0, 30}] (* Vaclav Kotesovec, Jan 04 2024 *)
a[0] := 1; a[n_] := (-1)^n*2^(n + 1)/(n*(n + 1)^2)*Sum[(-1/2)^k*Binomial[n + 1, k - 1]*Binomial[n + 1, k]*Binomial[n + 1, k + 1], {k , 1, n}]; Table[a[n], {n, 0, 25}] (* Detlef Meya, May 29 2024 *)

Prog

(SageMath)
def A368709(n): return PolyA359363(n, -2) // (-2) if n > 0 else 1
print([A368709(n) for n in range(0, 26)])  # Peter Luschny, Jan 04 2024
(Python)
def A368709(n):
    if n == 0: return 1
    return sum((-2)**k * v for k, v in enumerate(A359363Row(n))) // (-2)
print([A368709(n) for n in range(26)]) # Peter Luschny, Jan 04 2024

Crossrefs

Cf. A368708, A001181, A007724, A217800, A359363.

Sequence in context: A070520 A018579 A006486 * A273946 A070518 A263182

Adjacent sequences: A368706 A368707 A368708 * A368710 A368711 A368712

Keyword

sign

Author

Joerg Arndt, Jan 04 2024

Status

approved