OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..941
FORMULA
a(n) = (-2)^n*Sum_{k=0..n} A080247(n,k)/(-2)^k.
a(n) = ((8 - 4*n)*a(n-3) + (30 - 24*n)*a(n-2) + (17 - 37*n)*a(n-1))/(3*n + 3).
a(n) = [x^n] 2/(1 + sqrt(1 + 4*x*(x + 3)).
a(n) = [x^n] reverse((3*x^2 + x)/(1 - x^2))/x.
MAPLE
a := proc(n) option remember; if n < 3 then return [1, -3, 17][n+1] fi;
((8 - 4*n)*a(n-3) + (30 - 24*n)*a(n-2) + (17 - 37*n)*a(n-1))/(3*n + 3) end:
seq(a(n), n=0..20);
# Alternative:
gf := 2/(1 + sqrt(1 + 4*x*(x + 3))):
ser := series(gf, x, 24):
seq(coeff(ser, x, n), n=0..20);
# Or:
series((3*x^2 + x)/(1 - x^2), x, 24):
gfun:-seriestoseries(%, 'revogf'):
convert(%, polynom) / x: seq(coeff(%, x, n), n=0..20);
PROG
(SageMath)
R.<x> = PowerSeriesRing(QQ)
f = (3*x^2 + x)/(1 - x^2)
f.reverse().shift(-1).list()
(PARI) N=20; x='x+O('x^N); Vec(2/(1+sqrt(1+4*x*(x+3)))) \\ Seiichi Manyama, Feb 03 2020
CROSSREFS
KEYWORD
sign
AUTHOR
Peter Luschny, Jan 02 2020
STATUS
approved