|
%I #31 Jan 01 2018 21:20:11
%S 1,1,1,1,3,3,15,15,105,105,945,945,10395,10395,135135,135135,2027025,
%T 2027025,34459425,34459425,654729075,654729075,13749310575,
%U 13749310575,316234143225,316234143225,7905853580625,7905853580625,213458046676875,213458046676875
%N A001147 with each term repeated.
%C Normally such sequences are excluded from the OEIS, but I have made an exception for this one because so many variants of it have occurred in recent submissions.
%C For n>=2, a(n) = product of odd positive integers <=(n-1). - _Jaroslav Krizek_, Mar 21 2009
%C a(n) is, for n>=3, the number of way to choose floor((n-1)/2) disjoint pairs of items from n-1 items. It is then a fortiori the size of the conjugacy class of the reversal permutation [n-1,n-2,n-3,...,3,2,1]=(1 n-1)(2 n-2)(3 n-3)... in the symmetric group on n-1 elements. - _Karl-Dieter Crisman_, Nov 03 2009
%H G. C. Greubel, <a href="/A133221/b133221.txt">Table of n, a(n) for n = 0..800</a>
%F E.g.f.: x*U(0) where U(k)= 1 + (2*k+1)/(x - x^4/(x^3 + (2*k+2)*(2*k+3)/U(k+1))) ; (continued fraction, 3rd kind, 3-step). - _Sergei N. Gladkovskii_, Sep 25 2012
%F G.f.: 1+x*G(0), where G(k)= 1 + x*(2*k+1)/(1 - x/(x + 1/G(k+1))); (continued fraction). - _Sergei N. Gladkovskii_, Jun 07 2013
%t f[x_] := E^(x^2/2) + Sqrt[Pi/2]*Erfi[x/Sqrt[2]]; CoefficientList[ Series[f[x], {x, 0, 29}], x]*Range[0, 29]! (* _Jean-François Alcover_, Sep 25 2012, after _Sergei N. Gladkovskii_ *)
%t Table[(n - 1 - Mod[n, 2])!!, {n, 0, 20}] (* _Eric W. Weisstein_, Dec 31 2017 *)
%t Table[((2 n + (-1)^n - 3)/2)!!, {n, 0, 20}] (* _Eric W. Weisstein_, Dec 31 2017 *)
%o (Sage)
%o def Gauss_factorial(N, n): return mul(j for j in (1..N) if gcd(j, n) == 1)
%o def A133221(n): return Gauss_factorial(n-1, 2)
%o [A133221(n) for n in (0..29)] # _Peter Luschny_, Oct 01 2012
%o (PARI) a(n) = my(k = (2*n + (-1)^n - 3)/2); prod(i=0, (k-1)\2, k - 2*i) \\ _Iain Fox_, Dec 31 2017
%Y Cf. A055634.
%Y Appears in A161736. - _Johannes W. Meijer_, Jun 18 2009
%K nonn
%O 0,5
%A _N. J. A. Sloane_, Oct 13 2007
|
