OFFSET
1,3
LINKS
Paul Barry, Centered polygon numbers, heptagons and nonagons, and the Robbins numbers, arXiv:2104.01644 [math.CO], 2021.
FORMULA
Given g.f. A(x), then 1 = (1/A(x) + 1/A(-x)) / 2.
a(n) = -(-1)^n * binomial(n + m, n - m) / (2*m + 1) where m = floor(n/2) if n>0.
From Michael Somos, Apr 13 2012 (Start)
Reversion of A080956 with offset 1.
Hankel transform is A005161 omitting first 1.
n * a(n) = -(-1)^n * A099576(n-1). (End)
D-finite with recurrence +8*n*(n+1)*a(n) -36*n*(n-2)*a(n-1) +6*(-9*n^2+18*n-14)*a(n-2) +27*(3*n-7)*(3*n-8)*a(n-3)=0. - R. J. Mathar, Sep 24 2021
a(n) = (-1)^(n-1)*binomial(2*n, n-1)*hypergeom([-n+1, n], [-2*n], -1) / n. - Detlef Meya, Dec 26 2023
EXAMPLE
G.f. = x - x^2 + 2*x^3 - 3*x^4 + 7*x^5 - 12*x^6 + 30*x^7 - 55*x^8 + 143*x^9 + ...
MATHEMATICA
a[ n_] := With[ {m = Quotient[n, 2]}, If[n < 1, 0, -(-1)^n Binomial[n + m, n - m] / (2 m + 1)]]; (* Michael Somos, Oct 16 2015 *)
a[ n_] := If[n < 1, 0, SeriesCoefficient[ InverseSeries[ Series[(x - 2 x^2) / (1 - x)^3, {x, 0, n}]], {x, 0, n}]]; (* Michael Somos, Oct 16 2015 *)
a[n_] := (-1)^(n-1)*Binomial[2*n, n-1]*Hypergeometric2F1[-n+1, n, -2*n, -1] / n; Flatten[Table[a[n], {n, 1, 32}]] (* Detlef Meya, Dec 26 2023 *)
PROG
(PARI) {a(n) = my( m = n\2); if( n<1, 0, -(-1)^n * binomial( n + m, n - m) / (2 * m + 1))};
(PARI) {a(n) = if( n<1, 0, polcoeff( serreverse( (x - 2 * x^2) / (1 - x)^3 + x * O(x^n) ), n))};
(PARI) {a(n) = if( n<1, 0, polcoeff( 1 / ( 1 + 1 / serreverse( x - x^3 + x * O(x^n) )), n))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Nov 01 2007
STATUS
approved