A385310
Expansion of e.g.f. 1/(1 - 2 * x * cos(x))^(1/2).
Original entry on oeis.org
1, 1, 3, 12, 69, 500, 4455, 46928, 571977, 7914384, 122585355, 2100940864, 39470867469, 806555184448, 17808628411119, 422498774818560, 10717948285126545, 289501146405400832, 8295124400250875667, 251300745071590317056, 8025654235707259740885, 269482309052945201181696
Offset: 0
-
a185951(n, k) = binomial(n, k)/2^k*sum(j=0, k, (2*j-k)^(n-k)*binomial(k, j));
a001147(n) = prod(k=0, n-1, 2*k+1);
a(n) = sum(k=0, n, a001147(k)*I^(n-k)*a185951(n, k));
A352647
Expansion of e.g.f. 1/(1 - 3 * x * cos(x)).
Original entry on oeis.org
1, 3, 18, 153, 1728, 24315, 410400, 8079729, 181786752, 4601232243, 129402385920, 4003157532297, 135098815002624, 4939266681129963, 194472450526169088, 8203835046344538465, 369151362125290045440, 17649035213360472293091, 893431062200523039178752
Offset: 0
-
With[{m = 18}, Range[0, m]! * CoefficientList[Series[1/(1 - 3*x*Cos[x]), {x, 0, m}], x]] (* Amiram Eldar, Mar 26 2022 *)
-
my(N=20, x='x+O('x^N)); Vec(serlaplace(1/(1-3*x*cos(x))))
-
a(n) = if(n==0, 1, 3*sum(k=0, (n-1)\2, (-1)^k*(2*k+1)*binomial(n, 2*k+1)*a(n-2*k-1)));
Showing 1-2 of 2 results.