A071549 -id:A071549 - cp's OEIS frontend

A187783 De Bruijn's triangle, T(m,n) = (m*n)!/(n!^m) read by downward antidiagonals.

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 6, 6, 1, 1, 1, 20, 90, 24, 1, 1, 1, 70, 1680, 2520, 120, 1, 1, 1, 252, 34650, 369600, 113400, 720, 1, 1, 1, 924, 756756, 63063000, 168168000, 7484400, 5040, 1

Offset: 0

Robert G. Wilson v, Jan 05 2013

From Tilman Piesk, Oct 28 2014: (Start)
Number of permutations of a multiset that contains m different elements n times. These multisets have the signatures A249543(m,n-1) for m>=1 and n>=2.
In an m-dimensional Pascal tensor (the generalization of a symmetric Pascal matrix) P(x1,...,xn) = (x1+...+xn)!/(x1!*...*xn!), so the main diagonal of an m-dimensional Pascal tensor is D(n) = (m*n)!/(n!^m). These diagonals are the rows of this array (with m>0), which begins like this:
  m\n:0   1       2            3                4                    5
  0:  1   1       1            1                1                    1 ... A000012;
  1:  1   1       1            1                1                    1 ... A000012;
  2:  1   2       6           20               70                  252 ... A000984;
  3:  1   6      90         1680            34650               756756 ... A006480;
  4:  1  24    2520       369600         63063000          11732745024 ... A008977;
  5:  1 120  113400    168168000     305540235000      623360743125120 ... A008978;
  6:  1 720 7484400 137225088000 3246670537110000 88832646059788350720 ... A008979;
with columns: A000142 (n=1), A000680 (n=2), A014606 (n=3), A014608 (n=4), A014609 (n=5).
A089759 is the transpose of this matrix. A034841 is its diagonal. A141906 is its lower triangle. A120666 is the upper triangle of this matrix with indices starting from 1. A248827 are the diagonal sums (or the row sums of the triangle).
(End)

			T(3,5) = (3*5)!/(5!^3) = 756756 = A014609(3) = A006480(5) is the number of permutations of a multiset that contains 3 different elements 5 times, e.g., {1,1,1,1,1,2,2,2,2,2,3,3,3,3,3}.

Tilman Piesk, First 54 rows of the triangle, flattened
T. Chappell, A. Lascoux, S. Ole Warnaar, and W. Zudilin, Logarithmic and complex constant term identities, arXiv:1112.3130 [math.CO], 2012.
Tilman Piesk, Array for indices 0..16
Tilman Piesk, PHP code used to create the b-file
Tilman Piesk, Illustration of the multisets for m,n=0..4
Wikipedia, Permutations of multisets, Pascal matrix and simplex

Cf. A089759 (transposed), A141906 (subtriangle), A120666 (subtriangle transposed), A060538 (1st row/column removed).
Rows: A000012, A000984, A006480, A006480, A008978, A008979, A071549, A071550, A071551, A071552.
Columns: A000012, A000142, A000680, A014606, A014608, A014609, A248814, A172603, A172609, A172613.
Main diagonal gives: A034841.
Row sums of the triangle: A248827.

[Factorial(k*(n-k))/(Factorial(n-k))^k: k in [0..n], n in [0..10]]; // G. C. Greubel, Dec 26 2022

T[n_, k_]:= (k*n)!/(n!)^k; Table[T[n, k-n], {k, 9}, {n, 0, k-1}]//Flatten

def A187783(n,k): return gamma(k*(n-k)+1)/(factorial(n-k))^k
flatten([[A187783(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Dec 26 2022

T(m,n) = (m*n)!/(n!)^m.

A060540(m,n) = T(m,n)/m! . - R. J. Mathar, Jun 21 2023

Row m=0 prepended by Tilman Piesk, Oct 28 2014

A008979 a(n) = (6n)!/(n!)^6.

Original entry on oeis.org

1, 720, 7484400, 137225088000, 3246670537110000, 88832646059788350720, 2670177736637149247308800, 85722533226982363751829504000, 2889253496242619386328267523990000, 101097362223624462291180422369532000000, 3644153415887633116359073848179365185734400, 134567406165969006655507763343147223231094784000

Offset: 0

N. J. A. Sloane

Vincenzo Librandi, Table of n, a(n) for n = 0..100
R. Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.

Cf. A000897, A000984, A002894, A002897, A006480, A008977, A008978, A071549, A071550, A071551, A071552, A186420, A188662.

[Factorial(6*n)/Factorial(n)^6: n in [0..20]]; // Vincenzo Librandi, Aug 13 2014

seq( (6*n)!/(n!)^6, n=0..20); # G. C. Greubel, Feb 17 2020

Table[(6 n)!/(n)!^6, {n, 0, 20}] (* Vincenzo Librandi, Aug 13 2014 *)

vector(21, n, my(m=n-1); (6*m)!/(m!)^6 ) \\ G. C. Greubel, Feb 17 2020

[factorial(6*n)/factorial(n)^6 for n in (0..20)] # G. C. Greubel, Feb 17 2020

From Peter Bala, Jul 12 2016: (Start)
a(n) = binomial(2*n,n)*binomial(3*n,n)*binomial(4*n,n)*
  binomial(5*n,n)*binomial(6*n,n) = ( [x^n](1 + x)^(2*n) ) * ( [x^n](1 + x)^(3*n) ) * ( [x^n](1 + x)^(4*n) ) * ( [x^n](1 + x)^(5*n) ) * ( [x^n](1 + x)^(6*n) ) = [x^n](F(x)^(720*n)), where F(x) = 1 + x + 4478*x^2 + 53085611*x^3 + 926072057094*x^4 + 19977558181209910*x^5 + 493286693783478576177*x^6 + ... appears to have integer coefficients. For similar results see A000897, A002894, A002897, A006480, A008977, A008978, A186420 and A188662. (End)
a(n) ~ 3^(6*n+1/2)*4^(3*n-1)/(Pi*n)^(5/2). - Ilya Gutkovskiy, Jul 12 2016
From Peter Bala, Feb 14 2020: (Start)
a(m*p^k) == a(m*p^(k-1)) ( mod p^(3*k) ) for prime p >= 5 and positive integers m and k - apply Mestrovic, equation 39, p. 12.
a(n) = [(x*y*z*u*v)^n] (1 + x + y + z + u + v)^(6*n). (End)

A071550 a(n) = (8n)!/n!^8.

Original entry on oeis.org

1, 40320, 81729648000, 369398958888960000, 2390461829733887910000000, 18975581770994682860770223800320, 171889289584866507880743491472699801600

Offset: 0

Benoit Cloitre, May 30 2002

nonn

Vincenzo Librandi, Table of n, a(n) for n = 0..100

Sequences (k*n)!/n!^k: A000984 (k = 2), A006480 (k = 3), A008977 (k = 4), A008978 (k = 5), A008979 (k = 6), A071549 (k = 7), A071551 (k = 9), A071552 (k = 10).

[Factorial(8*n)/Factorial(n)^8: n in [0..20]]; // Vincenzo Librandi, Aug 13 2014

Table[(8 n)!/(n)!^8, {n, 0, 20}] (* Vincenzo Librandi, Aug 13 2014 *)

A071551 a(n) = (9n)!/n!^9.

Original entry on oeis.org

1, 362880, 12504636144000, 1080491954750208000000, 140810154080474667338550000000, 23183587808948692737291767860055162880, 4439413043841128802009762476941510771390464000

Offset: 0

Benoit Cloitre, May 30 2002

nonn

Vincenzo Librandi, Table of n, a(n) for n = 0..100

Sequences (k*n)!/n!^k: A000984 (k = 2), A006480 (k = 3), A008977 (k = 4), A008978 (k = 5), A008979 (k = 6), A071549 (k = 7), A071550 (k = 8), A071552 (k = 10).

[Factorial(9*n)/Factorial(n)^9: n in [0..20]]; // Vincenzo Librandi, Aug 13 2014

Table[(9n)!/(n)!^9, {n, 0, 20}] (* Vincenzo Librandi, Aug 13 2014 *)

A071552 a(n) = (10n)!/n!^10.

Original entry on oeis.org

1, 3628800, 2375880867360000, 4386797336285844480000000, 12868639981414579848070084500000000, 49120458506088132224064306071170476903628800

Offset: 0

Benoit Cloitre, May 30 2002

nonn

Vincenzo Librandi, Table of n, a(n) for n = 0..100

Sequences (k*n)!/n!^k: A000984 (k = 2), A006480 (k = 3), A008977 (k = 4), A008978 (k = 5), A008979 (k = 6), A071549 (k = 7), A071550 (k = 8), A071551 (k = 9).

[Factorial(10*n)/Factorial(n)^10: n in [0..20]]; // Vincenzo Librandi, Aug 13 2014

Table[(10n)!/(n)!^10, {n, 0, 20}] (* Vincenzo Librandi, Aug 13 2014 *)

A352652 a(n) = ( binomial(7n,2n)binomial(7n/2,2n)binomial(2n,n)^2 ) / binomial(7n/2,n)^2.

Original entry on oeis.org

1, 30, 2860, 343200, 45643500, 6435891280, 942422020540, 141696569678400, 21724714133822700, 3381208130986900500, 532553441617598475360, 84695057996350934903680, 13578009523892192555221500, 2191530567314796197691108600, 355765014009052303028935320000

Offset: 0

Peter Bala, Apr 03 2022

We write x! as shorthand for Gamma(x+1) and binomial(x,y) as shorthand for x!/(y!*(x-y)!) = Gamma(x+1)/(Gamma(y+1)*Gamma(x-y+1)).
Given two sequences of numbers c = (c_1, c_2, ..., c_K) and d = (d_1, d_2, ..., d_L) where c_1 + ... + c_K = d_1 + ... + d_L  we can define the factorial ratio sequence u_n(c, d) = (c_1*n)!*(c_2*n)!* ... *(c_K*n)!/ ( (d_1*n)!*(d_2*n)!* ... *(d_L*n)! ) and ask whether it is integral for all n >= 0. The integer L - K is called the height of the sequence. Bober completed the classification of integral factorial ratio sequences of height 1. Soundararajan gives many examples of two-parameter families of integral factorial ratio sequences of height 2. Usually, it is assumed that the c's and d's are integers but here we allow for some of the c's and d's to be rational numbers. See A276098 and the cross references for further examples of factorial ratios of this type.
Conjecture: the supercongruences a(n*p^k) == a(n*p^(k-1)) (mod p^(3*k)) hold for all primes p >= 5 and positive integers n and k. The case n = k = 1 is easily proved.

			Examples of supercongruences:
a(11) - a(1) = 84695057996350934903680 - 30 = 2*(5^2)*(11^3)*23*593* 3671*5693*4464799 == 0 (mod 11^3)
a(2*7) - a(2) = 355765014009052303028935320000 - 2860 = (2^2)*5*(7^3)*11* 269*3307*375101*14129010228023 == 0 (mod 7^3)

Peter Bala, Some integer ratios of factorials
J. W. Bober, Factorial ratios, hypergeometric series, and a family of step functions, 2007, arXiv:0709.1977v1 [math.NT], 2007; J. London Math. Soc., 79, Issue 2, (2009), 422-444.
K. Soundararajan, Integral factorial ratios: irreducible examples with height larger than 1, Phil. Trans. R. Soc. A378: 2018044, 2019.

Cf. A275654, A276098, A276100, A276101, A276102, A352651.

a := n -> if n = 0 then 1 elif n = 1 then 30 else
7*(5*n-2)*(5*n-4)*(5*n-6)*(5*n-8)*(7*n-1)*(7*n-3)*(7*n-5)*(7*n-9)*(7*n-11)*(7*n-13)/(3*n^2*(n-1)^2*(3*n-2)*(3*n-4)*(5*n-1)*(5*n- 3)*(5*n -7)*(5*n-9)) *a(n-2) end if:
seq(a(n), n = 0..20);

from math import factorial
from sympy import factorial2
def A352652(n): return int(factorial(7*n)*factorial2(5*n)**2//factorial(5*n)//factorial2(7*n)//factorial2(3*n)//factorial(n)**2) # Chai Wah Wu, Aug 08 2023

a(n) = (5/3)*Sum_{k = 0..n} (-1)^(n+k)*binomial(7*n,n-k)*binomial(5*n+k-1,k)^2 for n >= 1 (this formula shows 3*a(n) is integral; how to show a(n) is integral?).
a(n) = (5/3)*Sum_{k = 0..n} binomial(4*n-k-2,n-k)*binomial(5*n-1,k)^2 for n >= 1.
a(n) = (7*n)!*(5*n/2)!^2/((5*n)!*(7*n/2)!*(3*n/2)!*n!^2!).
a(n) = (5/3) * [x^n] ( (1 - x)^(2*n) * P(5*n-1,(1 + x)/(1 - x)) ) for n >= 1, where P(n,x) denotes the n-th Legendre polynomial.
a(n) = (5/3)*(-1)^n*binomial(7*n,n)*hypergeom([-n, 5*n, 5*n], [1, 6*n+1], 1) for n >= 1.
a(n) ~ sqrt(15)/Pi * 7^(7*n/2)/3^(3*n/2) * ( 1/(6*n) - 29/(945*n^2) + 841/(297675*n^3) + O(1/n^4) ).
a(n) = 7*(5*n-2)*(5*n-4)*(5*n-6)*(5*n-8)*(7*n-1)*(7*n-3)*(7*n-5)*(7*n-9)*(7*n-11)*(7*n-13)/(3*n^2*(n-1)^2*(3*n-2)*(3*n-4)*(5*n-1)*(5*n- 3)*(5*n -7)*(5*n-9)) * a(n-2) with a(0) = 1 and a(1) = 30.
a(n)*A275654(n) = (7*n)!/(n!^4*(3*n)!) = A071549(n)/A006480(n).
a(p) == 30 (mod p^3) for all primes p >= 5.

A367571 a(n) = Product_{k=0..n} (7*k)! / k!^7.

Original entry on oeis.org

1, 5040, 3432645216000, 626489905645044080640000000, 41646279370357699257014919153469440000000000000, 1200992054275801322636044235924808416678612164215512865177600000000000000

Offset: 0

Vaclav Kotesovec, Nov 23 2023

nonn

In general, for m > 1, Product_{k=0..n} (m*k)! / k!^m ~ A^(m - 1/m) * exp(m*n/2 - m/12 + 1/(12*m) - n/2) * m^(m*n^2/2 + m*n/2 - 1/(12*m) + n/2) * n^(-m*n/2 - m/3 + 1/(12*m) + n/2 + 1/4) * (2*Pi)^(-m*n/2 - m/4 + n/2 + 1/4) / Product_{j=1..m-1} Gamma(j/m)^(j/m), where A is the Glaisher-Kinkelin constant A074962.

Cf. A000178, A071549, A271947.

Cf. A007685, A367567, A367568, A367569, A367570.

Table[Product[(7*k)!/k!^7, {k, 0, n}], {n, 0, 10}]
Table[Product[Binomial[7*k,k] * Binomial[6*k,k] * Binomial[5*k,k] * Binomial[4*k,k] * Binomial[3*k,k] * Binomial[2*k,k], {k, 0, n}], {n, 0, 10}]

a(n) = Product_{k=0..n} binomial(7*k,k) * binomial(6*k,k) * binomial(5*k,k) * binomial(4*k,k) * binomial(3*k,k) * binomial(2*k,k).
a(n) = A271947(n) / A000178(n)^7.
a(n) ~ A^(48/7) * 7^(7*n^2/2 + 4*n - 1/84) * exp(3*n - 4/7) / (Gamma(1/7)^(1/7) * Gamma(2/7)^(2/7) * Gamma(3/7)^(3/7) * Gamma(4/7)^(4/7) * Gamma(5/7)^(5/7) * Gamma(6/7)^(6/7) * n^(3*n + 29/14) * (2*Pi)^(3*n + 3/2)), where A is the Glaisher-Kinkelin constant A074962.

A269125 Number of sequences with n copies each of 1,2,...,7 avoiding the pattern 12...7.

Original entry on oeis.org

1, 5039, 674641325, 171248808285596, 54391509124298646609, 18675035107760751336633904, 6541943110420293280017597411002, 2285085934263252199073238394141449534, 789007694236868030155290867470985247486401, 268568796699793322439632085075432144258760935329

Offset: 0

Alois P. Heinz, Feb 19 2016

nonn

Alois P. Heinz, Table of n, a(n) for n = 0..100

Row n=7 of A269129.

Cf. A071549, A268843.

a(n) = A071549(n) - A268843(n).

A198806 Number of closed paths of length n whose steps are 14th roots of unity, U_14(n).

Original entry on oeis.org

1, 0, 14, 0, 546, 0, 32900, 10080, 2570050, 2540160, 238935564, 465696000, 25142196156, 76886409600, 2900343069624, 12211317518400, 359067702643650, 1915829643087360, 47006105030584700, 300455419743198720, 6437718469449262996

Offset: 0

Simon Plouffe, Oct 30 2011

nonn

U_14(n), (comment in article): For each m >= 1, the sequence (U_m(N)), N >= 0 is P-recursive but is not algebraic when m > 2.

Andrew Howroyd, Table of n, a(n) for n = 0..200
Gilbert Labelle and Annie Lacasse, Closed paths whose steps are roots of unity, in FPSAC 2011, Reykjavik, Iceland DMTCS proc. AO, 2011, 599-610.

Cf. A002898, A070190, A071549.

seq(n)={Vec(serlaplace(sum(k=0, n, if(k,2,1)*(x^k*besseli(k, 2*x + O(x^(n-k+1)))/k!)^7)))} \\ Andrew Howroyd, Nov 01 2018

E.g.f.: BesselI(0,2*x)^7 + 2*Sum_{k>=1} BesselI(k,2*x)^7. - Andrew Howroyd, Nov 01 2018

cp's OEIS Frontend

A187783 De Bruijn's triangle, T(m,n) = (m*n)!/(n!^m) read by downward antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

SageMath

Formula

Extensions

A008979 a(n) = (6n)!/(n!)^6.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A071550 a(n) = (8n)!/n!^8.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

A071551 a(n) = (9n)!/n!^9.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

A071552 a(n) = (10n)!/n!^10.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

A352652 a(n) = ( binomial(7*n,2*n)*binomial(7*n/2,2*n)*binomial(2*n,n)^2 ) / binomial(7*n/2,n)^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Python

Formula

A367571 a(n) = Product_{k=0..n} (7*k)! / k!^7.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Formula

A269125 Number of sequences with n copies each of 1,2,...,7 avoiding the pattern 12...7.

Views

Author

Keywords

A352652 a(n) = ( binomial(7n,2n)binomial(7n/2,2n)binomial(2n,n)^2 ) / binomial(7n/2,n)^2.