login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A330803 Evaluation of the Big-Schröder polynomials at -1/2 and normalized with (-2)^n. 3
1, -3, 17, -123, 1001, -8739, 79969, -756939, 7349657, -72798003, 732681489, -7471545435, 77031538377, -801616570947, 8408819677377, -88821190791915, 943928491520249, -10085451034660947, 108275140773938545, -1167408859459660923, 12635538801834255401 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
LINKS
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
Sequence in context: A305307 A074543 A216314 * A267736 A129115 A246963
KEYWORD
sign
AUTHOR
Peter Luschny, Jan 02 2020
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 03:39 EDT 2024. Contains 371264 sequences. (Running on oeis4.)