VOOZH
about
URL: https://oeis.org/A049126
⇱ A049126 - OEIS
login
A049126
Revert transform of ((x - 1)(3x - 1))/(1 - 3x + x^3).
0
1, 1, 2, 6, 22, 90, 393, 1789, 8378, 40050, 194492, 956448, 4752519, 23822571, 120315345, 611644737, 3127389106, 16072642890, 82981119252, 430187414196, 2238469102212, 11687227631892, 61208286479382, 321465732705594
(
list
;
graph
;
refs
;
listen
;
history
;
text
;
internal format
)
OFFSET
1,3
LINKS
Table of n, a(n) for n=1..24.
Paul Barry,
Centered polygon numbers, heptagons and nonagons, and the Robbins numbers
, arXiv:2104.01644 [math.CO], 2021.
Index entries for reversions of series
FORMULA
Revert transform of
A122100
offset 1. -
Michael Somos
, May 11 2012
Recurrence: 12*(n-1)*n*(117*n - 604)*a(n) = 4*(n-1)*(4329*n^2 - 29719*n + 39384)*a(n-1) - 3*(23985*n^3 - 230056*n^2 + 679659*n - 630308)*a(n-2) + 3*(38961*n^3 - 437004*n^2 + 1609875*n - 1955812)*a(n-3) - 27*(n-4)*(2223*n^2 - 18496*n + 38133)*a(n-4) + 81*(n-5)*(n-4)*(117*n - 487)*a(n-5). -
Vaclav Kotesovec
, Jan 02 2021
a(n) ~ 3^(n - 1/4) * (2 + sqrt(3))^(n - 3/2) / (sqrt(Pi) * n^(3/2) * 2^(n + 1/2)). -
Vaclav Kotesovec
, Jan 02 2021
EXAMPLE
G.f. = x + x^2 + 2*x^3 + 6*x^4 + 22*x^5 + 90*x^6 + 393*x^7 + 1789*x^8 + ...
MATHEMATICA
Rest[CoefficientList[InverseSeries[Series[x*((x - 1)(3x - 1))/(1 - 3x + x^3), {x, 0, 40}], x], x]] (*
Vaclav Kotesovec
, Jan 02 2021 *)
PROG
(PARI) {a(n) = if( n<1, 0, polcoeff( serreverse(x * (1 - 4*x + 3*x^2) / (1 - 3*x + x^3) + x * O(x^n)), n))}; /*
Michael Somos
, May 11 2012 */
CROSSREFS
Cf.
A122100
.
Sequence in context:
A150270
A049136
A165523
*
A049134
A086456
A006318
Adjacent sequences:
A049123
A049124
A049125
*
A049127
A049128
A049129
KEYWORD
nonn
AUTHOR
Olivier Gérard
STATUS
approved