OFFSET
0,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..100
FORMULA
a(n) = 2*((2*n)!*Sum_{k=1..(2*n)} binomial(2*n+k,2*n)*Sum_{j=1..k} binomial(k,j)*(Sum_{l=0..(j-1)} (binomial(j,l)*Sum_{i=0..(j-l)/2} binomial(j-l,i)*(l-j+2*i)^(2*n-l+j)*(-1)^(n-i)))/(2*n-l+j)!))))), a(0)=1.
a(n) ~ 2^(2*n+1) * n^(2*n) / (3^(1/4) * exp(2*n) * (sqrt(3) - Pi/3)^(2*n+1/2)). - Vaclav Kotesovec, Jan 26 2014
MATHEMATICA
Table[(CoefficientList[InverseSeries[Series[-x+2*Sin[x], {x, 0, 31}], x], x]*Range[0, 31]!)[[n]], {n, 2, 30, 2}] (* Vaclav Kotesovec, Jan 26 2014 *)
PROG
(Maxima) a(n):=if n=0 then 1 else 2*((2*n)!*sum(binomial(2*n+k, 2*n)*sum(binomial(k, j)*(sum((binomial(j, l)*sum(binomial(j-l, i)*(l-j+2*i)^(2*n-l+j)*(-1)^(n-i), i, 0, (j-l)/2))/(2*n-l+j)!, l, 0, j-1)), j, 1, k), k, 1, 2*n));
(PARI) seq(n)={my(p=serlaplace(serreverse(2*sin(x + O(x^(2*n+2))) - x))); vector(n+1, i, polcoef(p, 2*i-1))} \\ Andrew Howroyd, Jan 04 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Kruchinin, Feb 04 2012
EXTENSIONS
Terms a(13) and beyond from Andrew Howroyd, Jan 04 2020
STATUS
approved