A061137 Number of degree-n odd permutations of order dividing 6.

0, 0, 1, 3, 6, 30, 270, 1386, 6048, 46656, 387180, 2469060, 17204616, 158065128, 1903506696, 18887563800, 163657221120, 2095170230016, 30792968596368, 346564643468976, 3905503235814240, 58609511127871200, 866032039742528736

Offset: 0

Author

Vladeta Jovovic, Apr 14 2001

Keywords

• easy
• nonn

Links

• Robert Israel, Table of n, a(n) for n = 0..522
• Lev Glebsky, Melany Licón, Luis Manuel Rivera, On the number of even roots of permutations, arXiv:1907.00548 [math.CO], 2019.
• Robert Israel, Recurrence for this sequence

Crossrefs

Cf. A000085, A001470, A001472, A052501, A053496-A053505, A001189, A001471, A001473, A061121-A061128, A000704, A061129-A061132, A048099, A051695, A061133-A061135, A001465, A061136-A061140.

Programs

Magma

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x + x^3/3)*Sinh(x^2/2 + x^6/6) )); [0,0] cat [Factorial(n+1)*b[n]: n in [1..m-2]]; // G. C. Greubel, Jul 02 2019

Maple

Egf:= exp(x + x^3/3)*sinh(x^2/2 + x^6/6):
S:= series(Egf,x,31):
seq(coeff(S,x,j)*j!,j=0..30); # Robert Israel, Jul 13 2018

Mathematica

With[{m=30}, CoefficientList[Series[Exp[x + x^3/3]*Sinh[x^2/2 + x^6/6], {x, 0, m}], x]*Range[0,m]!] (* Vincenzo Librandi, Jul 02 2019 *)

PARI

my(x='x+O('x^30)); concat([0,0], Vec(serlaplace( exp(x + x^3/3)*sinh(x^2/2 + x^6/6) ))) \\ G. C. Greubel, Jul 02 2019

Sage

m = 30; T = taylor(exp(x + x^3/3)*sinh(x^2/2 + x^6/6), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Jul 02 2019

Formula

E.g.f.: exp(x + x^3/3)*sinh(x^2/2 + x^6/6).

Linear recurrence of order 12 whose coefficients are polynomials in n of degree up to 15: see link. - Robert Israel, Jul 13 2018