login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A052503 Number of permutations sigma of [2n] without fixed points such that sigma^4 = Id. 4

%I

%S 1,1,9,105,2625,76545,3440745,176080905,12034447425,922995698625,

%T 87505195602825,9203114782686825,1141501848477415425,

%U 155540530213013570625,24232951756530007115625,4112826185329479728735625,781060320618828163499210625

%N Number of permutations sigma of [2n] without fixed points such that sigma^4 = Id.

%H G. C. Greubel, <a href="/A052503/b052503.txt">Table of n, a(n) for n = 0..280</a>

%H INRIA Algorithms Project, <a href="http://ecs.inria.fr/services/structure?nbr=28">Encyclopedia of Combinatorial Structures 28</a>

%F a(n) = (2n)! * [x^(2n)] exp(x^2/2 + x^4/4).

%F D-finite with recurrence a(n) +(-2*n+1)*a(n-1) -2*(n-1)*(2*n-1)*(2*n-3)*a(n-2)=0, with a(0)=1, a(1)=1, a(2)=9. - Corrected by _R. J. Mathar_, Feb 20 2020 to skip zeros.

%F a(n) = 2^n*Gamma(n+1/2)*A047974(n)/Pi^(1/2). - _Mark van Hoeij_, Oct 30 2011.

%p spec := [S,{S=Set(Union(Cycle(Z,card=2),Cycle(Z,card=4)))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);

%t With[{nmax = 40}, CoefficientList[Series[Exp[x^2*(2 + x^2)/4], {x, 0, nmax}], x]*(Range[0, nmax])!][[1 ;; -1 ;; 2]] (* _G. C. Greubel_, May 14 2019 *)

%o (PARI) x='x+O('x^40); v=Vec(serlaplace( exp(x^2/2 + x^4/4) )); vector(#v\2, n, v[2*n-1]) \\ _G. C. Greubel_, May 14 2019

%o (MAGMA) m:=40; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x^2/2 + x^4/4) )); [Factorial(2*n-2)*b[2*n-1]: n in [1..Floor((m-2)/2)]]; // _G. C. Greubel_, May 14 2019

%o (Sage) m = 40; T = taylor(exp(x^2/2 + x^4/4), x, 0, 2*m+2); [factorial(2*n)*T.coefficient(x, 2*n) for n in (0..m)] # _G. C. Greubel_, May 14 2019

%Y Cf. A001472, A261317, A261381.

%Y Bisection of column k=4 of A261430.

%K easy,nonn,changed

%O 0,3

%A encyclopedia(AT)pommard.inria.fr, Jan 25 2000

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 25 04:32 EST 2020. Contains 332217 sequences. (Running on oeis4.)