|
%I #35 Jan 01 2019 07:11:02
%S 1,6,8,60,180,1680,8064,90720,604800,7983360,68428800,1037836800,
%T 10897286400,186810624000,2324754432000,44460928512000,
%U 640237370572800,13516122267648000,221172909834240000,5109094217170944000,93666727314800640000,2350183339898634240000,47726800133326110720000
%N a(n) is the reciprocal of the coefficient of x^n in the power series of the function defined by ((1+2x)*exp(x) + 3*exp(-x) - 4)/ (4x^2).
%H Muniru A Asiru, <a href="/A322544/b322544.txt">Table of n, a(n) for n = 0..445</a>
%F a(n) = (n+2)!/(3*floor(n/2)-n+2).
%F a(n) = (4*(n+2)!)/(2n+5+3*(-1)^n).
%F a(n) = 4/([x^n]((exp(x)*(1+2x)+3*exp(-x)-4)/x^2)).
%F a(n) = (n+2)!/(A028242(n)+1)
%F a(n) = (n+2)!/A030451(n+1)
%p a:=n->factorial(n+2)/(3*floor(n/2)-n+2): seq(a(n),n=0..25); # _Muniru A Asiru_, Dec 20 2018
%t Table[4*Factorial[n + 2]/(2*n + 5 + 3*(-1)^n), {n, 0, 25}]
%t (* or *)
%t Function[x, 1/x] /@
%t CoefficientList[Series[(Exp[x]/4 + 3/4*Exp[-x] + x/2*Exp[x] - 1)/x^2, {x, 0, 20}], x]
%o (PARI) a(n)={(4*(n+2)!)/(5 + 3*(-1)^n + 2*n)} \\ _Andrew Howroyd_, Dec 14 2018
%o (PARI) my(x='x + O('x^30)); Vec(apply(x->1/x, ((1+2*x)*exp(x) + 3*exp(-x) - 4)/ (4*x^2))) \\ _Michel Marcus_, Dec 19 2018
%o (GAP) List([0..25],n->(4*Factorial(n+2))/(2*n+5+3*(-1)^n)); # _Muniru A Asiru_, Dec 20 2018
%Y Cf. A060593 (even bisection, shifted), A028242 (denominator minus 1), A030451 (denominator, shifted), A107991 (Expansion of a similar function)
%K nonn
%O 0,2
%A _Pierre-Alain Sallard_, Dec 14 2018
| 