A167569 -id:A167569 - cp's OEIS frontend

A047053 a(n) = 4^n * n!.

1, 4, 32, 384, 6144, 122880, 2949120, 82575360, 2642411520, 95126814720, 3805072588800, 167423193907200, 8036313307545600, 417888291992371200, 23401744351572787200, 1404104661094367232000, 89862698310039502848000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

Original name was "Quadruple factorial numbers".
For n >= 1, a(n) is the order of the wreath product of the cyclic group C_4 and the symmetric group S_n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
Number of n X n monomial matrices with entries 0, +/-1, +/-i.
a(n) is the product of the positive integers <= 4*n that are multiples of 4. - Peter Luschny, Jun 23 2011
Also, a(n) is the number of signed permutations of length 2*n that are equal to their reverse-complements. (See the Hardt and Troyka reference.) - Justin M. Troyka, Aug 13 2011.
Pi^n/a(n) is the volume of a 2*n-dimensional ball with radius 1/2. - Peter Luschny, Jul 24 2012
Equals the first right hand column of A167557, and also equals the first right hand column of A167569. - Johannes W. Meijer, Nov 12 2009
a(n) is the order of the group U_n(Z[i]) = {A in M_n(Z[i]): A*A^H = I_n}, the group of n X n unitary matrices over the Gaussian integers. Here A^H is the conjugate transpose of A. - Jianing Song, Mar 29 2021

			G.f. = 1 + 4*x + 32*x^2 + 384*x^3 + 6144*x^4 + 122880*x^5 + 2949120*x^6 + ...

Vincenzo Librandi, Table of n, a(n) for n = 0..100
R. B. Brent, Generalizing Tuenter's Binomial Sums, J. Int. Seq. 18 (2015), # 15.3.2.
CombOS - Combinatorial Object Server, Generate colored permutations
R. Coquereaux and J.-B. Zuber, Maps, immersions and permutations, arXiv preprint arXiv:1507.03163 [math.CO], 2015-2016. Also J. Knot Theory Ramifications 25, 1650047 (2016), DOI.
Sylvie Corteel and Lauren Williams, Tableaux Combinatorics for the Asymmetric Exclusion Process II, arXiv:0810.2916 [math.CO], 2008-2009.
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).
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 492.
Alexsandar Petojevic, The Function vM_m(s; a; z) and Some Well-Known Sequences, Journal of Integer Sequences, 5 (2002), Article 02.1.7.
M. D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications, J. Int. Seq. 13 (2010), Article 10.6.7, p. 39.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, 9 (2006), Article 06.1.1.

Cf. A000142, A000165, A007696, A008545, A032031, A047058, A087299, A092042, A092616.

a(n)= A051142(n+1, 0) (first column of triangle).

[4^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Jul 20 2011

A047053:= n -> mul(k, k = select(k-> k mod 4 = 0, [$1..4*n])): seq(A047053(n), n = 0.. 16); # Peter Luschny, Jun 23 2011

a[n_]:= With[{m=2n}, If[ m<0, 0, m!*SeriesCoefficient[1 +Sqrt[Pi]*x*Exp[x^2]*Erf[x], {x, 0, m}]]]; (* Michael Somos, Jan 03 2015 *)
Table[4^n n!,{n,0,20}] (* Harvey P. Dale, Sep 19 2021 *)

PARI
```
a(n)=4^n*n!;
```

a(n) = 4^n * n!.
E.g.f.: 1/(1 - 4*x).
Integral representation as the n-th moment of a positive function on a positive half-axis: a(n) = Integral_{x=0..oo} x^n*exp(-x/4)/4, n >= 0. This representation is unique. - Karol A. Penson, Jan 28 2002 [corrected by Jason Yuen, May 04 2025]
Sum_{k>=0} (-1)^k/(2*k + 1)^n = (-1)^n * n * (PolyGamma[n-1, 1/4] - PolyGamma[n-1, 3/4]) / a(n) for n > 0. - Joseph Biberstine (jrbibers(AT)indiana.edu), Jul 27 2006
a(n) = Sum_{k=0..n} C(n,k)*(2k)!*(2(n-k))!/(k!(n-k)!) = Sum_{k=0..n} C(n,k)*A001813(k)*A001813(n-k). - Paul Barry, May 04 2007
E.g.f.: With interpolated zeros, 1 + sqrt(Pi)*x*exp(x^2)*erf(x). - Paul Barry, Apr 10 2010
From Gary W. Adamson, Jul 19 2011: (Start)
a(n) = sum of top row terms of M^n, M = an infinite square production matrix as follows:
  2, 2, 0, 0, 0, 0, ...
  4, 4, 4, 0, 0, 0, ...
  6, 6, 6, 6, 0, 0, ...
  8, 8, 8, 8, 8, 0, ...
  ... (End)
G.f.: 1/(1 - 4*x/(1 - 4*x/(1 - 8*x/(1 - 8*x/(1 - 12*x/(1 - 12*x/(1 - 16*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 08 2012
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - 8*x*(k + 1)/(8*x*(k + 1) - 1 + 8*x*(k + 1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
G.f.: 1/Q(0), where Q(k) = 1 - 4*x*(2*k + 1) - 16*x^2*(k + 1)^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 28 2013
a(n) = A000142(n) * A000302(n). - Michel Marcus, Nov 28 2013
a(n) = A087299(2*n). - Michael Somos, Jan 03 2015
D-finite with recurrence: a(n) - 4*n*a(n-1) = 0. - R. J. Mathar, Jan 27 2020
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/4) (A092042).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/4) (A092616). (End)

Edited by Karol A. Penson, Jan 22 2002

A034177 a(n) is the n-th quartic factorial number divided by 4.

Original entry on oeis.org

1, 8, 96, 1536, 30720, 737280, 20643840, 660602880, 23781703680, 951268147200, 41855798476800, 2009078326886400, 104472072998092800, 5850436087893196800, 351026165273591808000, 22465674577509875712000, 1527665871270671548416000, 109991942731488351485952000

Offset: 1

Wolfdieter Lang

			G.f. = x + 8*x^2 + 96*x^3 + 1536*x^4 + 30720*x^5 + 737820*x^6 + ...

Michael De Vlieger, Table of n, a(n) for n = 1..365
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 513.
Norihiro Nakashima and Shuhei Tsujie, Enumeration of Flats of the Extended Catalan and Shi Arrangements with Species, arXiv:1904.09748 [math.CO], 2019.

Cf. A007696, A000407, A034176. First column of triangle A048786.
A052570 is an essentially identical sequence. - Philippe Deléham, Sep 18 2008
Equals the second right hand column of A167569 divided by 2. - Johannes W. Meijer, Nov 12 2009

List([1..20], n-> 4^(n-1)*Factorial(n) ); # G. C. Greubel, Aug 15 2019

[4^(n-1)*Factorial(n): n in [1..20]]; // G. C. Greubel, Aug 15 2019

[seq(n!*4^(n-1), n=1..16)]; # Zerinvary Lajos, Sep 23 2006

Array[4^(# - 1) #! &, 16] (* Michael De Vlieger, May 30 2019 *)

vector(20, n, 4^(n-1)*n!) \\ G. C. Greubel, Aug 15 2019

[4^(n-1)*factorial(n) for n in (1..20)] # G. C. Greubel, Aug 15 2019

4*a(n) = (4*n)(!^4) = Product_{j=1..n} 4*j = 4^n * n!.
E.g.f.: (-1 + 1/(1-4*x))/4.
D-finite with recurrence: a(n) -4*n*a(n-1)=0. - R. J. Mathar, Feb 24 2020
From Amiram Eldar, Jan 08 2022: (Start)
Sum_{n>=1} 1/a(n) = 4*(exp(1/4)-1).
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*(1-exp(-1/4)). (End)

A167560 The ED2 array read by ascending antidiagonals.

Original entry on oeis.org

1, 2, 1, 6, 4, 1, 24, 16, 6, 1, 120, 80, 32, 8, 1, 720, 480, 192, 54, 10, 1, 5040, 3360, 1344, 384, 82, 12, 1, 40320, 26880, 10752, 3072, 680, 116, 14, 1, 362880, 241920, 96768, 27648, 6144, 1104, 156, 16, 1

Offset: 1

Johannes W. Meijer, Nov 10 2009

The coefficients in the upper right triangle of the ED2 array (m>n) were found with the a(n,m) formula while the coefficients in the lower left triangle of the ED2 array (m<=n) were found with the recurrence relation, see below. We use for the array rows the letter n (>=1) and for the array columns the letter m (>=1).
The ED2 array is related to the EG1 matrix, see A162005, because sum(EG1(2*m-1,n) * z^(2*m-1), m=1..infinity) = ((2*n-1)!/(4^(n-1)*(n-1)!^2))*int(sinh(y*(2*z))/cosh(y)^(2*n), y=0..infinity).
For the ED1, ED3 and ED4 arrays see A167546, A167572 and A167584.

			The ED2 array begins with:
1, 1, 1, 1, 1, 1, 1, 1, 1, 1
2, 4, 6, 8, 10, 12, 14, 16, 18, 20
6, 16, 32, 54, 82, 116, 156, 202, 254, 312
24, 80, 192, 384, 680, 1104, 1680, 2432, 3384, 4560
120, 480, 1344, 3072, 6144, 11160, 18840, 30024, 45672, 66864
720, 3360, 10752, 27648, 61440, 122880, 226800, 392832, 646128, 1018080

Johannes W. Meijer, The four Escher-Droste arrays, Mar 08 2013.

A000012, A005843 (n>=1), 2*A104249 (n>=1), A167561, A167562 and A167563 equal the first sixth rows of the array.
A000142 equals the first column of the array.
A047053 equals the a(n, n) diagonal of the array.
2*A034177 equals the a(n+1, n) diagonal of the array.
A167570 equals the a(n+2, n) diagonal of the array,
A167564 equals the row sums of the ED2 array read by antidiagonals.
A167565 is a triangle related to the a(n) formulas of the rows of the ED2 array.
A167568 is a triangle related to the GF(z) formulas of the rows of the ED2 array.
A167569 is the lower left triangle of the ED2 array.
Cf. A162005 (EG1 triangle).
Cf. A167546 (ED1 array), A167572 (ED3 array), A167584 (ED4 array).

nmax:=10; mmax:=10; for n from 1 to nmax do for m from 1 to n do a(n,m) := 4^(m-1)*(m-1)!*(n+m-1)!/(2*m-1)! od; for m from n+1 to mmax do a(n,m):= n! + sum((-1)^(k-1)*binomial(n-1,k)*a(n,m-k),k=1..n-1) od; od: for n from 1 to nmax do for m from 1 to n do d(n,m):=a(n-m+1,m) od: od: T:=1: for n from 1 to nmax do for m from 1 to n do a(T):= d(n,m): T:=T+1: od: od: seq(a(n),n=1..T-1);
# alternative
A167560 := proc(n,m)
    option remember ;
    if m > n then
        n!+add( (-1)^(k-1)*binomial(n-1,k)*procname(n,m-k),k=1..n-1) ;
    else
        4^(m-1)*(m-1)!*(n+m-1)!/(2*m-1)! ;
    end if;
end proc:
seq( seq(A167560(d-m,m),m=1..d-1),d=2..12) ; # R. J. Mathar, Jun 28 2024

nmax = 10; mmax = 10; For[n = 1, n <= nmax, n++, For[m = 1, m <= n, m++, a[n, m] = 4^(m - 1)*(m - 1)!*((n + m - 1)!/(2*m - 1)!)]; For[m = n + 1, m <= mmax, m++, a[n, m] = n! + Sum[(-1)^(k - 1)*Binomial[n - 1, k]*a[n, m - k], {k, 1, n - 1}]]; ]; For[n = 1, n <= nmax, n++, For[m = 1, m <= n, m++, d[n, m] = a[n - m + 1, m]]; ]; t = 1; For[n = 1, n <= nmax, n++, For[m = 1, m <= n, m++, a[t] = d[n, m]; t = t + 1]]; Table[a[n], {n, 1, t - 1}] (* Jean-François Alcover, Dec 20 2011, translated from Maple *)

a(n,m) = ((m-1)!/((m-n-1)!))*int(sinh(y*(2*n))/(cosh(y))^(2*m),y=0..infinity) for m>n.
The (n-1)-differences of the n-th array row lead to the recurrence relation
sum((-1)^k*binomial(n-1,k)*a(n-1,m-k),k=0..n-1) = n!
which in its turn leads to, see A167569,
a(n,m) = 4^(m-1)*(m-1)!*(n+m-1)!/(2*m-1)! if m<=n.

cp's OEIS Frontend

A167570 The third right hand column of triangle A167569.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A167571 The row sums of triangle A167569.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A047053 a(n) = 4^n * n!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A034177 a(n) is the n-th quartic factorial number divided by 4.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A167560 The ED2 array read by ascending antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula