cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A167558 The second right hand column of triangle A167557.

Original entry on oeis.org

1, 12, 160, 2688, 55296, 1351680, 38338560, 1238630400, 44920995840, 1807409479680, 79906524364800, 3850733459865600, 200907832688640000, 11282983883794022400, 678650586195610828800, 43527244493925384192000
Offset: 2

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Comments

Triangle A167557 is the lower left triangle of the ED1 array A167546.

Crossrefs

A167546 is the ED1 array.
Equals the second right hand column of triangle A167557.

Programs

  • Mathematica
    Table[4^(n-2) (n-2)! (2n-3)!/(2n-4)!,{n,2,20}] (* Harvey P. Dale, Aug 21 2013 *)

Formula

a(n) = 4^(n-2)*(n-2)!*(2*n-3)!/(2*n-4)!.
(-2*n+5)*a(n) +4*(2*n-3)*(n-2)*a(n-1)=0. - R. J. Mathar, Apr 01 2015

A167559 The row sums of triangle A167557.

Original entry on oeis.org

1, 5, 46, 598, 10056, 207960, 5111856, 145721520, 4728528000, 172134253440, 6948393580800, 308058678385920, 14882061698380800, 778136819786726400, 43783908789625804800, 2638014817853286144000, 169455373956672285081600
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Comments

Triangle A167557 is the lower left triangle of the ED1 array A167546.

Crossrefs

A167546 is the ED1 array.
Equals row sums of triangle A167557.

Programs

  • Mathematica
    Table[Sum[4^(m - 1) *(m - 1)! * (n + m - 2)!/(2*m - 2)!, {m, 1, n}], {n, 1, 10}] (* G. C. Greubel, Jun 15 2016 *)

Formula

a(n) = Sum_{m=1..n} (4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)!).

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

Original entry on oeis.org

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

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Comments

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

Examples

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

Crossrefs

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

Programs

  • Magma
    [4^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Jul 20 2011
  • Maple
    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
  • Mathematica
    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!;
    

Formula

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)

Extensions

Edited by Karol A. Penson, Jan 22 2002

A167546 The ED1 array read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 2, 4, 1, 6, 12, 7, 1, 24, 48, 32, 10, 1, 120, 240, 160, 62, 13, 1, 720, 1440, 960, 384, 102, 16, 1, 5040, 10080, 6720, 2688, 762, 152, 19, 1, 40320, 80640, 53760, 21504, 6144, 1336, 212, 22, 1
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Comments

The coefficients in the upper right triangle of the ED1 array (m > n) were found with the a(n,m) formula while the coefficients in the lower left triangle of the ED1 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).
Our procedure for finding the coefficients in the lower left triangle can be compared with the procedure that De Smit and Lenstra used to fill in the hole in the middle of 'The Print Gallery' by M. C. Escher, see the links. In this lithograph Escher made use of the so-called Droste effect, hence we propose to call this square array of numbers the ED1 array.
For the ED2, ED3 and ED4 arrays see A167560, A167572 and A167584.

Examples

			The ED1 array begins with:
1, 1, 1, 1, 1, 1, 1, 1, 1, 1
1, 4, 7, 10, 13, 16, 19, 22, 25, 28
2, 12, 32, 62, 102, 152, 212, 282, 362, 452
6, 48, 160, 384, 762, 1336, 2148, 3240, 4654, 6432
24, 240, 960, 2688, 6144, 12264, 22200, 37320, 59208, 89664
120, 1440, 6720, 21504, 55296, 122880, 245640, 452880, 783144, 1285536
		

Crossrefs

A000012, A016777, 2*A005891, A167547, A167548 and A167549 equal the first sixth rows of the array.
A000142 equals the first column of the array.
A167550 equals the a(n, n+1) diagonal of the array.
A047053 equals the a(n, n) diagonal of the array.
A167558 equals the a(n+1, n) diagonal of the array.
A167551 equals the row sums of the ED1 array read by antidiagonals.
A167552 is a triangle related to the a(n) formulas of rows of the ED1 array.
A167556 is a triangle related to the GF(z) formulas of the rows of the ED1 array.
A167557 is the lower left triangle of the ED1 array.
Cf. A068424 (the (m-1)!/(m-n-1)! factor), A007680 (the (2*n-1)*(n-1)! factor).
Cf. A167560 (ED2 array), A167572 (ED3 array), A167584 (ED4 array).

Programs

  • Maple
    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-1+m-1)!/(2*m-2)! od; for m from n+1 to mmax do a(n,m):= (2*n-1)*(n-1)! + 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);
  • Mathematica
    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 - 1 + m - 1)!/(2*m - 2)!)]; For[m = n + 1, m <= mmax, m++, a[n, m] = (2*n - 1)*(n - 1)! + 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 *)

Formula

a(n,m) = (2*(m-1)!/(m-n-1)!)*Integral_{y>=0} sinh(y*(2*n-1))/cosh(y)^(2*m-1) for m > n.
The (n-1)-differences of the n-th array row lead to the recurrence relation
Sum_{k=0..n-1} (-1)^k*binomial(n-1,k)*a(n,m-k) = (2*n-1)*(n-1)!
which in its turn leads to, see also A167557,
a(n,m) = 4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)! for m <= n.
Showing 1-4 of 4 results.