login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A144828 Partial products of successive terms of A017113; a(0)=1. 15
1, 4, 48, 960, 26880, 967680, 42577920, 2214051840, 132843110400, 9033331507200, 686533194547200, 57668788341964800, 5305528527460761600, 530552852746076160000, 57299708096576225280000, 6646766139202842132480000, 824199001261152424427520000, 108794268166472120024432640000 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
a(n) is the number of signed permutations of length 4n that are equal to their reverse-inverses. Note that the reverse-inverse of a permutation is equivalent to a 90-degree rotation of the permutation's diagram (see the Hardt and Troyka link). - Justin M. Troyka, Aug 11 2011
Define the bar operation as an operation on signed permutation that flips the sign of each entry. Then a(n) is the number of signed permutations of length 2n that are equal to the bar of their inverses and equal to their reverse-complements (see the Hardt and Troyka link). - Justin M. Troyka, Aug 11 2011
LINKS
A. Hardt and J. M. Troyka, Restricted symmetric signed permutations, Pure Mathematics and Applications, Vol. 23 (No. 3, 2012), pp. 179-217.
A. Hardt and J. M. Troyka, Slides (associated with the Hardt and Troyka reference above).
FORMULA
a(n) = Sum_{k=0..n} A132393(n,k)*4^k*8^(n-k).
a(n) = A052714(n+1). - R. J. Mathar, Oct 01 2008
a(n) = 2^n *(2*n)! / n!. - Justin M. Troyka, Aug 11 2011
G.f.: 1/(1-4x/(1-8x/(1-12x/(1-16x/(1-20x/(1-24x/(1-28x/(1-32x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
a(n) = (-4)^n*Sum_{k=0..n} 2^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
E.g.f.: 1/sqrt(1-8*x). - Philippe Deléham, May 14 2015
a(n) = 4^n * A001147(n). - Philippe Deléham, May 14 2015
a(n) = 8^n * Gamma(n + 1/2) / sqrt(Pi). - Daniel Suteu, Jan 06 2017
0 = a(n)*(8*a(n+1) - a(n+2)) + a(n+1)*(+a(n+1)) and a(n) = (-1)^n / a(-n) for all n in Z. - Michael Somos, Jan 06 2017
a(n) = 2^n * (n+1)! * Catalan(n). - G. C. Greubel, Apr 02 2021
Sum_{n>=0} 1/a(n) = 1 + e^(1/8)*sqrt(Pi)*erf(1/(2*sqrt(2)))/(2*sqrt(2)), where erf is the error function. - Amiram Eldar, Dec 20 2022
EXAMPLE
a(0)=1, a(1)=4, a(2)=4*12=48, a(3)=4*12*20=960, a(4)=4*12*20*28=26880, ...
Since a(1) = 4, there are 4 signed permutations of 4 that are equal to their reverse-inverses. These are: (+2,+4,+1,+3), (+3,+1,+4,+2), (-2,-4,-1,-3), (-3,-1,-4,-2). - Justin M. Troyka, Aug 11 2011
G.f. = 1 + 4*x + 48*x^2 + 960*x^3 + 26880*x^4 + 967680*x^5 + 42577920*x^6 + ...
MAPLE
A144828:= n-> 2^n*n!*binomial(2*n, n); seq(A144828(n), n=0..30); # G. C. Greubel, Apr 02 2021
MATHEMATICA
Table[4^n (2 n - 1)!!, {n, 0, 15}] (* Vincenzo Librandi, May 14 2015 *)
Join[{1}, FoldList[Times, (8*Range[0, 20]+4)]] (* Harvey P. Dale, Dec 01 2015 *)
PROG
(Magma) [2^k *Factorial(2*k) / Factorial(k): k in [0..20]]; // Vincenzo Librandi, Aug 11 2011
(PARI) a(n)=binomial(2*n, n)*n!<<n \\ Charles R Greathouse IV, Jan 17 2012
(PARI) {a(n) = if( n<0, (-1)^n / a(-n), 2^n *(2*n)! / n!)}; /* Michael Somos, Jan 06 2017 */
(Sage) [2^n*factorial(n+1)*catalan_number(n) for n in (0..30)] # G. C. Greubel, Apr 02 2021
CROSSREFS
Essentially the same as A052714. - N. J. A. Sloane, Feb 03 2013
Sequences of the form m^(n-1)*n!*Catalan(n-1): A001813 (m=1), A052714 (or this sequence) (m=2), A221954 (m=3), A052734 (m=4), A221953 (m=5), A221955 (m=6).
Sequence in context: A162676 A141119 A052714 * A138448 A071221 A198038
KEYWORD
nonn,easy
AUTHOR
Philippe Deléham, Sep 21 2008
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 08:13 EDT 2024. Contains 371265 sequences. (Running on oeis4.)