OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: 1 / (1 - x / (1 - x / (1 - 2*x / (1 - x / (2 - 3*x / (1 - 2*x / (3 - 4*x / ... ))))))).
D-finite with recurrence: 0 = (4*n + 6) * a(n) - (17*n + 27) * a(n+1) + (24*n + 42) * a(n+2) - (9*n + 21) * a(n+3) + (n + 3) * a(n+4). - Sign flipped by R. J. Mathar, Feb 16 2020
0 = a(n) * (16*a(n+1) - 74*a(n+2) + 120*a(n+3) - 66*a(n+4) + 10*a(n+5))+ a(n+1) * (-62*a(n+1) + 361*a(n+2) - 480*a(n+3) + 265*a(n+4) - 41*a(n+5)) + a(n+2) * (-342*a(n+2) + 615*a(n+3) - 335*a(n+4) + 54*a(n+5)) + a(n+3) * (-90*a(n+3) + 75*a(n+4) - 15*a(n+5)) + a(n+4) * (-3*a(n+4) + a(n+5)).
a(n) = A129775(n) if n>0.
HANKEL transform is A000012.
INVERT transform is A073525.
EXAMPLE
G.f. = 1 + x + 2*x^2 + 6*x^3 + 21*x^4 + 78*x^5 + 298*x^6 + 1157*x^7 + ...
MATHEMATICA
a[ n_] := SeriesCoefficient[ 2 / (2 - x - x / Sqrt[1 - 4 x]), {x, 0, n}]
PROG
(PARI) {a(n) = if( n<0, 0, polcoeff( 2 / (2 - x - x / sqrt(1 - 4*x + x * O(x^n))), n))}
(Magma) m:=25; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R!(2/( 2-x - x/Sqrt(1-4*x)))); // G. C. Greubel, Aug 07 2018
CROSSREFS
KEYWORD
dead
AUTHOR
Michael Somos, Jan 09 2014
STATUS
approved