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-10 of 19 results. Next

A325305 Irregular triangular array, read by rows: T(n,k) is the sum of the products of distinct multinomial coefficients (n_1 + n_2 + n_3 + ...)!/(n_1! * n_2! * n_3! * ...) taken k at a time, where (n_1, n_2, n_3, ...) runs over all integer partitions of n (n >= 0, 0 <= k <= A070289(n)).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 2, 1, 10, 27, 18, 1, 47, 718, 4416, 10656, 6912, 1, 246, 20545, 751800, 12911500, 100380000, 304200000, 216000000, 1, 1602, 929171, 260888070, 39883405500, 3492052425000, 177328940580000, 5153150631600000, 82577533320000000, 669410956800000000, 2224399449600000000, 1632586752000000000, 1, 11271
Offset: 0

Views

Author

Petros Hadjicostas, Sep 05 2019

Keywords

Comments

This array was inspired by R. H. Hardin's recurrences for the columns of array A212855. Row n has length A070289(n) + 1.
This array differs from array A309951 starting at row n = 7. Array A309951 calculates a similar sum of products of multinomial coefficients, but the multinomial coefficients do not have to be distinct. Row n of array A309951 has length A000041(n) + 1, i.e., one more than the number of partitions of n.
Let P_n be the set of all lists a = (a_1, a_2,..., a_n) of integers a_i >= 0, i = 1,..., n such that 1*a_1 + 2*a_2 + ... + n*a_n = n; i.e., P_n is the set all integer partitions of n. (We use a different notation for partitions than the one in the name of T(n,k).) Then |P_n| = A000041(n) for n >= 0.
For n = 1..6, all the multinomial coefficients n!/((1!)^a_1 * (2!)^a_2 * ... * (n!)^a^n) corresponding to lists (a_1,...,a_n) in P_n are distinct; that is, A000041(n) = A070289(n) for n = 1..6.
For n = 7, the partitions (a_1, a_2, a_3, a_4, a_5, a_6, a_7) = (0, 2, 1, 0, 0, 0, 0) (i.e., 2 + 2 + 3) and (a_1, a_2, a_3, a_4, a_5, a_6, a_7) = (3, 0, 0, 1, 0, 0, 0) (i.e., 1 + 1 + 1 + 4) give the same multinomial coefficient: 210 = 7!/(2!2!3!) = 7!/(1!1!1!4!). Hence, A000041(7) > A070289(7).
Looking at the multinomial coefficients of the integer partitions of n = 8, 9, 10 on pp. 831-832 of Abramowitz and Stegun (1964), we see that, even in these cases, we have A000041(n) > A070289(n).

Examples

			Triangle begins as follows:
  [n=0]: 1,   1;
  [n=1]: 1,   1;
  [n=2]: 1,   3,     2;
  [n=3]: 1,  10,    27,     18;
  [n=4]: 1,  47,   718,   4416,    10656,      6912;
  [n=5]: 1, 246, 20545, 751800, 12911500, 100380000, 304200000, 216000000;
  ...
For example, when n = 3, the integer partitions of 3 are 3, 1+2, 1+1+1, and the corresponding multinomial coefficients are 3!/3! = 1, 3!/(1!2!) = 3, and 3!/(1!1!1!) = 6. They are all distinct. Then T(n=3, k=0) = 1, T(n=3, k=1) = 1 + 3 + 6 = 10, T(n=3, k=2) = 1*3 + 1*6 + 3*6 = 27, and T(n=3, k=3) = 1*3*6 = 18.
Consider the list [1, 7, 21, 35, 42, 105, 140, 210, 420, 630, 840, 1260, 2520, 5040] of the A070289(7) = 15 - 1 = 14 distinct multinomial coefficients corresponding to the 15 integer partitions of 7. Then  T(7,0) = 1, T(7,1) = 11271 (sum of the coefficients), T(7,2) = 46169368 (sum of products of every two different coefficients), T(7,3) = 92088653622 (sum of products of every three different coefficients), and so on. Finally, T(7,14) = 2372695722072874920960000000000 = product of these coefficients.
		

Crossrefs

Programs

  • Maple
    g:= proc(n, i) option remember; `if`(n=0 or i=1, [n!], [{map(x->
          binomial(n, i)*x, g(n-i, min(n-i, i)))[], g(n, i-1)[]}[]])
        end:
    b:= proc(n, m) option remember; `if`(n=0, 1,
          expand(b(n-1, m)*(g(m$2)[n]*x+1)))
        end:
    T:= n->(p->seq(coeff(p, x, i), i=0..degree(p)))(b(nops(g(n$2)), n)):
    seq(T(n), n=0..7);  # Alois P. Heinz, Sep 05 2019
  • Mathematica
    g[n_, i_] := g[n, i] = If[n == 0 || i == 1, {n!}, Union[Map[Function[x, Binomial[n, i] x], g[n - i, Min[n - i, i]]], g[n, i - 1]]];
    b[n_, m_] := b[n, m] = If[n == 0, 1, b[n - 1, m] (g[m, m][[n]] x + 1)];
    T[n_] := CoefficientList[b[Length[g[n, n]], n], x];
    T /@ Range[0, 7] // Flatten (* Jean-François Alcover, May 06 2020, after Maple *)

Formula

Sum_{k=0..A070289(n)} (-1)^k * T(n,k) = 0.

A212855 T(n,k) = number of n X k arrays with rows being permutations of 0..k-1 and no column j greater than column j-1 in all rows (n, k >= 1).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 19, 7, 1, 1, 211, 163, 15, 1, 1, 3651, 8983, 1135, 31, 1, 1, 90921, 966751, 271375, 7291, 63, 1, 1, 3081513, 179781181, 158408751, 7225951, 45199, 127, 1, 1, 136407699, 53090086057, 191740223841, 21855093751, 182199871, 275563, 255, 1
Offset: 1

Views

Author

R. H. Hardin, May 28 2012

Keywords

Comments

In other words, there are no "column rises", where a "column rise" means a pair of adjacent columns where each entry in the left column is strictly less than the adjacent entry in the right column.
This is R(n,k,0) in [Abramson-Promislow].
From Petros Hadjicostas, Sep 09 2019: (Start)
As stated above, in the notation of Abramson and Promislow (1978), we have T(n,k) = R(n, k, t=0).
Let P_k be the set of all lists a = (a_1, a_2, ..., a_k) of integers a_i >= 0, i = 1, ..., k, such that 1*a_1 + 2*a_2 + ... + k*a_k = k; i.e., P_k is the set all integer partitions of k. Then |P_k| = A000041(k).
From Eq. (6), p. 248, in Abramson and Promislow (1978), with t=0, we get T(n,k) = Sum_{a in P_k} (-1)^(k - Sum_{j=1..k} a_j) * (a_1 + a_2 + ... + a_k)!/(a_1! * a_2! * ... * a_k!) * (k! / ((1!)^a_1 * (2!)^a_2 * ... * (k!)^a_k))^n.
The integer partitions of k = 1..10 are listed on pp. 831-832 of Abramowitz and Stegun (1964). We see that, for k = 1..6, the corresponding multinomial coefficients k! / ((1!)^a_1 * (2!)^a_2 * ... * (k!)^a_k) are all distinct; that is, A070289(k) = A000041(k) and A309951(k,s) = A325305(k,s) for s = 0..A000041(k). For 7 <= k <= 10, this is not true anymore; i.e., A070289(k) < A000041(k) for 7 <= k <= 10 (and we conjecture that this is the case for all k >= 7).
From the theory of difference equations, we see that Abramson and Promislow's Eq. (6) on p. 248 (with t=0) implies that Sum_{s = 0..A070289(k)} (-1)^s * A325305(k,s) * T(n-s,k) = 0 for n >= A070289(k) + 1. For k = 1..5, these recurrences give R. H. Hardin's empirical recurrences shown in the Formula section below.
We also have Sum_{s = 0..A000041(k)} (-1)^s * A309951(k,s) * T(n-s,k) = 0 for n >= A000041(k) + 1, but for k >= 7, the recurrence we get (for column k) may not necessarily be minimal.
To derive the recurrence for row n, let y=0 in Eq. (8), p. 249, of Abramson and Promislow (1978). We get 1 + Sum_{k >= 1} T(n,k)*x^k/(k!)^n = 1/f_n(-x), where f_n(x) = Sum_{i >= 0} (x^i/(i!)^n). Matching coefficients, we get Sum_{s = 1..k} T(n,s) * (-1)^(s-1) * binomial(k,s)^n = 1, from which the recurrence in the Formula section follows.
(End)

Examples

			Some solutions for n=3 and k=4:
  2 1 3 0    1 3 0 2    3 0 2 1    1 3 0 2    1 3 2 0
  2 0 1 3    1 3 0 2    3 1 2 0    1 0 3 2    1 3 0 2
  2 3 0 1    3 0 2 1    2 3 1 0    2 0 3 1    3 1 0 2
Table starts:
  1  1     1         1             1                  1                       1
  1  3    19       211          3651              90921                 3081513
  1  7   163      8983        966751          179781181             53090086057
  1 15  1135    271375     158408751       191740223841         429966316953825
  1 31  7291   7225951   21855093751    164481310134301     2675558106868421881
  1 63 45199 182199871 2801736968751 128645361626874561 14895038886845467640193
		

Crossrefs

Cf. A000012 (row 1), A000275 (row 2), A212856 (row 3), A212857 (row 4), A212858 (row 5), A212859 (row 6), A212860 (row 7).
Cf. A000012 (column 1), A000225 (column 2), A212850 (column 3), A212851 (column 4), A212852 (column 5), A212853 (column 6), A212854 (column 7).
Cf. A000041, A070289 (order of minimal recurrence for column k), A192721, A212806 (main diagonal), A309951, A325305.

Programs

  • Maple
    A212855_row := proc(m,len) proc(n,m) sum(z^k/k!^m, k = 0..infinity);
    series(%^x, z=0, n+1): n!^m*coeff(%,z,n); [seq(coeff(%,x,k),k=0..n)] end;
    seq(add(abs(k), k=%(j,m)), j=1..len) end:
    for n from 1 to 6 do A212855_row(n,7) od; # Peter Luschny, May 26 2017
    # second Maple program:
    T:= proc(n, k) option remember; `if`(k=0, 1, -add(
          binomial(k, j)^n*(-1)^j*T(n, k-j), j=1..k))
        end:
    seq(seq(T(n, 1+d-n), n=1..d), d=1..10);  # Alois P. Heinz, Apr 26 2020
  • Mathematica
    rows = 9;
    row[m_, len_] := Module[{p, s0, s1, s2}, p = Function[{n, m0}, s0 = Sum[ z^k/k!^m0, {k, 0, n}]; s1 = Series[s0^x, {z, 0, n+1}] // Normal; s2 = n!^m0*Coefficient[s1, z, n]; Table[Coefficient[s2, x, k], {k, 0, n}]]; Table[Sum[Abs[k], {k, p[j, m]}], {j, 1, len}]];
    T = Table[row[n, rows+1], {n, 1, rows}];
    Table[T[[n-k+1, k]], {n, 1, rows}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Feb 27 2018, after Peter Luschny *)

Formula

Empirical recurrence for column k:
k=1: a(n) = 1*a(n-1).
k=2: a(n) = 3*a(n-1) - 2*a(n-2).
k=3: a(n) = 10*a(n-1) - 27*a(n-2) + 18*a(n-3).
k=4: a(n) = 47*a(n-1) - 718*a(n-2) + 4416*a(n-3) - 10656*a(n-4) + 6912*a(n-5).
k=5: a(n) = 246*a(n-1) - 20545*a(n-2) + 751800*a(n-3) - 12911500*a(n-4) + 100380000*a(n-5) - 304200000*a(n-6) + 216000000*a(n-7).
[All the "empirical" recurrences above are correct. See the comments above.]
From Benoit Jubin, May 29 2012: (Start)
T(n,1) = T(1,n) = 1.
T(n,2) = 2^n - 1 since the only n X 2 matrix with rows permutations of {0,1} which has a column rise is the one where all rows are [0,1].
(k!)^n*(1 - (k-1)/2^n) <= T(n,k) <= (k!)^n (the first inequality is (11) in the Abramson-Promislow reference, the second is trivial). (End)
For r >= 1, A(n, r) = Sum_{k=0..n} |[x^k] n!^r [z^n] S(r, z)^x| where S(r, z) = Sum_{k>=0} z^k/k!^r. - Peter Luschny, Feb 27 2018
From Petros Hadjicostas, Sep 09 2019: (Start)
Recurrence for column k: Sum_{s = 0..A070289(k)} (-1)^s * A325305(k,s) * T(n-s,k) = 0 for n >= A070289(k) + 1.
Recurrence for row n: T(n,k) = (-1)^(k-1) + Sum_{s = 1..k-1} T(n,s) * (-1)^(k-s-1) * binomial(k,s)^n for k >= 1.
(End)
Sum_{k>=1} T(n,k)*z^k/(k!)^n = 1/E_n(-z) -1 where E_n(z) = Sum_{k>=0} z^k/(k!)^n. - Geoffrey Critzer, Apr 28 2023

A212852 Number of n X 5 arrays with rows being permutations of 0..4 and no column j greater than column j-1 in all rows.

Original entry on oeis.org

1, 3651, 966751, 158408751, 21855093751, 2801736968751, 347190069843751, 42328368099218751, 5119530150996093751, 616756797369980468751, 74155772004699902343751, 8907394925520999511718751
Offset: 1

Views

Author

R. H. Hardin, May 28 2012

Keywords

Comments

Column 5 of A212855.
From Petros Hadjicostas, Sep 06 2019: (Start)
Let P_5 be the set of all lists b = (b_1, b_2, b_3, b_4, b_5) of integers b_i >= 0, i = 1, ..., 5, such that 1*b_1 + 2*b_2 + 3*b_3 + 4*b_4 + 5*b_5 = 5; i.e., P_5 is the set all integer partitions of 5. Then |P_5| = A000041(5) = 7.
From Eq. (6), p. 248, in Abramson and Promislow (1978), we get a(n) = A212855(n,5) = Sum_{b in P_5} (-1)^(5 - Sum_{j=1..5} b_j) * (b_1 + b_2 + b_3 + b_4 + b_5)!/(b_1! * b_2! * b_3! * b_4! * b_5!) * (5! / ((1!)^b_1 * (2!)^b_2 * (3!)^b_3 * (4!)^b_4 * (5!)^b_5))^n.
The integer partitions of 5 are listed on p. 831 of Abramowitz and Stegun (1964). We see that the corresponding multinomial coefficients 5! / ((1!)^b_1 * (2!)^b_2 * (3!)^b_3 * (4!)^b_4 * (5!)^b_5) are all distinct; that is, A070289(5) = A000041(5) = 7.
Using the integer partitions of 5 and the above formula for a(n), we may derive R. J. Mathar's formula below.
(End)

Examples

			Some solutions for n=3
..0..3..1..2..4....0..2..4..1..3....0..1..4..3..2....0..2..3..4..1
..1..0..4..3..2....1..0..3..2..4....1..3..0..4..2....0..4..3..1..2
..2..4..1..3..0....1..2..0..4..3....3..1..4..0..2....4..0..1..3..2
		

Crossrefs

Programs

  • Mathematica
    T[n_, k_] := T[n, k] = If[k == 0, 1, -Sum[Binomial[k, j]^n*(-1)^j*T[n, k - j], {j, 1, k}]];
    a[n_] := T[n, 5];
    Table[a[n], {n, 1, 12}] (* Jean-François Alcover, Apr 01 2024, after Alois P. Heinz in A212855 *)

Formula

Empirical: a(n) = 246*a(n-1) -20545*a(n-2) +751800*a(n-3) -12911500*a(n-4) +100380000*a(n-5) -304200000*a(n-6) +216000000*a(n-7).
Empirical: a(n) = -2*5^n + 3*20^n - 4*60^n + 120^n + 3*30^n - 2*10^n + 1. R. J. Mathar, Jun 25 2012
Sum_{s = 0..7} (-1)^s * A325305(5, s) * a(n-s) = 0 for n >= 8. (This is the same as R. H. Hardin's recurrence above, and it follows from Eq. (6) (with t=0), p. 248, in Abramson and Promislow (1978).) - Petros Hadjicostas, Sep 06 2019

A212853 Number of n X 6 arrays with rows being permutations of 0..5 and no column j greater than column j-1 in all rows.

Original entry on oeis.org

1, 90921, 179781181, 191740223841, 164481310134301, 128645361626874561, 96426023622482278621, 70816637331790329140481, 51492108377805402906874141, 37256471170472317193421713601, 26890352949868734582700237312861
Offset: 1

Views

Author

R. H. Hardin, May 28 2012

Keywords

Comments

Column 6 of A212855.
From Petros Hadjicostas, Sep 08 2019: (Start)
Let P_6 be the set of all lists b = (b_1, b_2, b_3, b_4, b_5, b_6) of integers b_i >= 0, i = 1, ..., 6, such that 1*b_1 + 2*b_2 + 3*b_3 + 4*b_4 + 5*b_5 + 6*b_6 = 6; i.e., P_6 is the set all integer partitions of 6. Then |P_6| = A000041(6) = 11.
From Eq. (6), p. 248, in Abramson and Promislow (1978), with t=0, we get a(n) = A212855(n,6) = Sum_{b in P_6} (-1)^(6-Sum_{j=1..6} b_j) * (b_1 + b_2 + b_3 + b_4 + b_5 + b_6)!/(b_1! * b_2! * b_3! * b_4! * b_5! * b_6!) * (6! / ((1!)^b_1 * (2!)^b_2 * (3!)^b_3 * (4!)^b_4 * (5!)^b_5 * (6!)^b_6))^n.
The integer partitions of 6 are listed on p. 831 of Abramowitz and Stegun (1964). We see that the corresponding multinomial coefficients 6! / ((1!)^b_1 * (2!)^b_2 * (3!)^b_3 * (4!)^b_4 * (5!)^b_5 * (6!)^b_6) are all distinct; that is, A070289(6) = A000041(6) = 11 and A309951(6,s) = A325305(6,s) for s = 0..11. (Compare with the comments for A212854.)
Using the information about partitions of 6 in Eq. (6) (with t=0), p. 248, of Abramson and Promislow (1978), we may derive the explicit equation for a(n) shown below.
Using standard results from the theory of difference equations (since the solution is known explicitly), we may derive R. H. Hardin's empirical recurrence. The recurrence is equivalent to Sum_{s = 0..11} (-1)^s * A325305(6,s) * a(n-s) = 0 for n >= 12.
(End)

Examples

			Some solutions for n=3:
  0 3 1 4 2 5   0 3 1 4 2 5   0 3 1 4 2 5   0 3 1 4 2 5
  3 0 2 4 5 1   1 3 0 4 5 2   4 0 3 1 2 5   0 1 5 2 3 4
  1 2 4 0 3 5   5 0 4 2 3 1   2 1 5 4 3 0   3 1 5 0 4 2
		

Crossrefs

Programs

  • Mathematica
    T[n_, k_] := T[n, k] = If[k == 0, 1, -Sum[Binomial[k, j]^n*(-1)^j*T[n, k - j], {j, 1, k}]];
    a[n_] := T[n, 6];
    Table[a[n], {n, 1, 12}] (* Jean-François Alcover, Apr 01 2024, after Alois P. Heinz in A212855 *)

Formula

Empirical: a(n) = 1602*a(n-1) - 929171*a(n-2) + 260888070*a(n-3) - 39883405500*a(n-4) + 3492052425000*a(n-5) - 177328940580000*a(n-6) + 5153150631600000*a(n-7) - 82577533320000000*a(n-8) + 669410956800000000*a(n-9) - 2224399449600000000*a(n-10) + 1632586752000000000*a(n-11) for n >= 12. [It is correct; see the comments above.]
a(n) = -1 + 2*6^n + 2*15^n + 20^n - 3*30^n - 6*60^n - 90^n + 4*120^n + 6*180^n - 5*360^n + 720^n for n >= 1. - Petros Hadjicostas, Sep 08 2019

A212854 Number of n X 7 arrays with rows being permutations of 0..6 and no column j greater than column j-1 in all rows.

Original entry on oeis.org

1, 3081513, 53090086057, 429966316953825, 2675558106868421881, 14895038886845467640193, 78785944892341703819175577, 406643086764765052892275303425, 2073826171428339544452057104498041
Offset: 1

Views

Author

R. H. Hardin, May 28 2012

Keywords

Comments

From Petros Hadjicostas, Aug 25 2019: (Start)
We have a(m) = R(m,n=7,t=0) = A212855(m,7) for m >= 1, where R(m,n,t) = LHS of Eq. (6) of Abramson and Promislow (1978, p. 248).
Let P_7 be the set of all lists b = (b_1, b_2,..., b_7) of integers b_i >= 0, i = 1, ..., 7 such that 1*b_1 + 2*b_2 + ... + 7*b_7 = 7; i.e., P_7 is the set all integer partitions of 7. Then |P_7| = A000041(7) = 15.
We have a(m) = A212855(m,7) = Sum_{b in P_7} (-1)^(7 - Sum_{j=1..7} b_j) * (b_1 + b_2 + ... + b_7)!/(b_1! * b_2! * ... * b_7!) * (7! / ((1!)^b_1 * (2!)^b_2 * ... * (7!)^b_7))^m.
The integer partitions of 7 are listed on p. 831 of Abramowitz and Stegun (1964). We see that, when (b_1, b_2, ..., b_7) = (0, 2, 1, 0, 0, 0, 0) or (3, 0, 0, 1, 0, 0, 0) (i.e., we have the partitions 2+2+3 and 1+1+1+4), the corresponding multinomial coefficients are 210 = 7!/(2!2!3!) = 7!/(1!1!1!4!), so the number of terms in the expression for a(m) is |P_7| - 1 = 15 - 1 = 14 (see below in the Formula section).
Let M_7 := [1, 7, 21, 35, 42, 105, 140, 210, 420, 630, 840, 1260, 2520, 5040] be the A070289(7) = 15 - 1 = 14 distinct multinomial coefficients corresponding to the 15 integer partitions of 7 in P_7. The characteristic equation of the recurrence for a(m) is f(x) := Product_{r in M_7} (x-r) = Sum_{i = 0..14} (-1)^{14-i} * c_i * x^i. It turns out that c_{14} = 1, c_{13} = 11271, c_{12} = 46169368, c_{11} = 92088653622, and so on (see R. H. Hardin's recurrence below), and c_0 = 2372695722072874920960000000000 = product of elements in M_7.
It follows that a(m) satisfies the recurrence Sum_{i = 0..14} (-1)^{14-i} * c_i * a(m-i) = 0, which is equivalent to R. H. Hardin's empirical recurrence below.
If we count the multinomial coefficient 210 twice in the characteristic equation (since it corresponds to two different integer partitions of 7) then we get (x-210)*f(x) = Sum_{i = 0..15} (-1)^{15-i} * d_i * x^i, where (d_0, d_1, ..., d_15) is row k = 7 in irregular triangular array A309951. We have d_{15} = 1, d_{14} = 11481, ..., d_0 = 498266101635303733401600000000000 (see Alois P. Heinz's b-file for A309951 with entries 37 to 52). Note that d_0 = 210 * c_0.
We then have Sum_{s = 0..15} (-1)^s * A309951(7, s) * a(m-s) = 0 for m >= 16. The latter recurrence is of order 15, and it is not minimal (as opposed to the one below by R. H. Hardin, which is of order 14 and minimal).
(End)

Examples

			Some solutions for n=3
..0..3..4..1..5..2..6....0..3..4..1..5..2..6....0..3..4..1..5..2..6
..1..0..3..5..2..6..4....1..0..3..2..4..5..6....1..0..4..2..5..6..3
..5..2..1..0..6..3..4....4..6..5..1..0..3..2....2..4..0..6..3..5..1
		

Crossrefs

Programs

  • Mathematica
    T[n_, k_] := T[n, k] = If[k == 0, 1, -Sum[Binomial[k, j]^n*(-1)^j*T[n, k - j], {j, 1, k}]];
    a[n_] := T[n, 7];
    Table[a[n], {n, 1, 12}] (* Jean-François Alcover, Apr 01 2024, after Alois P. Heinz in A212855 *)

Formula

Empirical: a(n) = 11271*a(n-1) -46169368*a(n-2) +92088653622*a(n-3) -100896701243149*a(n-4) +64220064122517975*a(n-5) -24283767237355832850*a(n-6) +5479502670227877007500*a(n-7) -734487423806273666445000*a(n-8) +57519812656973505919500000*a(n-9) -2547756421856270328438000000*a(n-10) +60760702040873540340600000000*a(n-11) -700874827794270417254400000000*a(n-12) +3015300813467611878720000000000*a(n-13) -2372695722072874920960000000000*a(n-14). [It is correct; see the comments above and one of the formulas below.]
a(n) = 1 - 2*7^n - 2*21^n - 2*35^n + 3*42^n + 6*105^n + 3*140^n - 210^n - 12*420^n - 4*630^n + 5*840^n + 10*1260^n - 6*2520^n + 5040^n. - Petros Hadjicostas, Aug 25 2019
Sum_{s = 0..14} (-1)^s * A325305(7, s) * a(n-s) = 0 for n >= 15. (This is the same as R. H. Hardin's recurrence above, and it follows from Eq. (6), p. 248, in Abramson and Promislow (1978) with t=0.) - Petros Hadjicostas, Sep 06 2019

A309951 Irregular triangular array, read by rows: T(n,k) is the sum of the products of multinomial coefficients (n_1 + n_2 + n_3 + ...)!/(n_1! * n_2! * n_3! * ...) taken k at a time, where (n_1, n_2, n_3, ...) runs over all integer partitions of n (n >= 0, 0 <= k <= A000041(n)).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 2, 1, 10, 27, 18, 1, 47, 718, 4416, 10656, 6912, 1, 246, 20545, 751800, 12911500, 100380000, 304200000, 216000000, 1, 1602, 929171, 260888070, 39883405500, 3492052425000, 177328940580000, 5153150631600000, 82577533320000000, 669410956800000000, 2224399449600000000, 1632586752000000000, 1, 11481
Offset: 0

Views

Author

Petros Hadjicostas, Aug 25 2019

Keywords

Comments

This array was inspired by R. H. Hardin's recurrences for the columns of array A212855. Rows k=1 to k=5 are due to him, while the remaining rows were computed by Alois P. Heinz.
Row n has length A000041(n) + 1, i.e., one more than the number of partitions of n.
Let R(m,n) := R(m,n,t=0) = A212855(m,n) for m,n >= 1, where R(m,n,t) = LHS of Eq. (6) of Abramson and Promislow (1978, p. 248).
Let P_n be the set of all lists a = (a_1, a_2,..., a_n) of integers a_i >= 0, i = 1,..., n such that 1*a_1 + 2*a_2 + ... + n*a_n = n; i.e., P_n is the set all integer partitions of n. (We use a different notation for partitions than the one in the name of T(n,k).) Then |P_n| = A000041(n) for n >= 0.
We have R(m,n) = A212855(m,n) = Sum_{a in P_n} (-1)^(n - Sum_{j=1..n} a_j) * (a_1 + a_2 + ... + a_n)!/(a_1! * a_2! * ... * a_n!) * (n! / ((1!)^a_1 * (2!)^a_2 * ... * (n!)^a_n))^m.
The recurrence of R. H. Hardin for column n of array A212855 is Sum_{s = 0..|P_n|} (-1)^s * T(n,s) * R(m-s,n) = 0 for n >= 1 and m >= |P_n| + 1.
The above recurrence is correct for all n >= 1, but it is not always a minimal one. For example, it seems to be the minimal one for n = 1,...,6, but not for n = 7 (see A212854). It seems to be minimal whenever every two different partitions of n give different multinomial coefficients.
For n = 7, the partitions (a_1, a_2, a_3, a_4, a_5, a_6, a_7) = (0, 2, 1, 0, 0, 0, 0) (i.e., 2 + 2 + 3) and (a_1, a_2, a_3, a_4, a_5, a_6, a_7) = (3, 0, 0, 1, 0, 0, 0) (i.e., 1 + 1 + 1 + 4) give the same multinomial coefficient: 210 = 7!/(2!2!3!) = 7!/(1!1!1!4!). Hence, to find the minimal recurrence for n = 7, we count 210 only once in the set of multinomial coefficients: 1, 7, 21, 35, 42, 105, 140, 210, 420, 630, 840, 1260, 2520, 5040. Then the absolute value of the coefficient of a(n-1) in the minimal recurrence is the sum of these multinomial coefficients (i.e., 11271); the absolute value of the coefficient of a(n-2) in the minimal recurrence is the sum of products of every two of them (i.e., 46169368), and so on.
Looking at the multinomial coefficients of the integer partitions of n = 8, 9, 10 on pp. 831-832 of Abramowitz and Stegun (1964), we see that, even in these cases, the above recurrence is not the minimal one. The number of distinct multinomial coefficients among the integer partitions of n is given by A070289.

Examples

			Triangle begins as follows:
  [n=0]: 1,   1;
  [n=1]: 1,   1;
  [n=2]: 1,   3,     2;
  [n=3]: 1,  10,    27,     18;
  [n=4]: 1,  47,   718,   4416,    10656,      6912;
  [n=5]: 1, 246, 20545, 751800, 12911500, 100380000, 304200000, 216000000;
  ...
For example, when n = 3, the integer partitions of 3 are 3, 1+2, 1+1+1, and the corresponding multinomial coefficients are 3!/3! = 1, 3!/(1!2!) = 3, and 3!/(1!1!1!) = 6. Then T(n=3, k=0) = 1, T(n=3, k=1) = 1 + 3 + 6 = 10, T(n=3, k=2) = 1*3 + 1*6 + 3*6 = 27, and T(n=3, k=3) = 1*3*6 = 18.
Since |P_3| = A000041(3) = 3, the recurrence of _R. H. Hardin_ for column n = 3 of array A212855 is T(3,0)*R(m,3) - T(3,1)*R(m-1,3) + T(3,2)*R(m-2,3) - T(3,3)*R(m-3,3) = 0; i.e., R(m,3) - 10*R(m-1,3) + 27*R(m-2,3) - 18*R(m-3,3) = 0 for m >= 4. We have the initial conditions R(m=1,3) = 1, R(m=2,3) = 19, and R(m=3,3) = 163. Thus, R(m,3) = 6^m - 2*3^m + 1 = A212850(m) for m >= 1. See the documentation of array A212855.
		

Crossrefs

Rightmost terms in rows give A309972.

Programs

  • Maple
    g:= proc(n, i) option remember; `if`(n=0 or i=1, [n!], [map(x->
          binomial(n, i)*x, g(n-i, min(n-i, i)))[], g(n, i-1)[]])
        end:
    b:= proc(n, m) option remember; `if`(n=0, 1,
          expand(b(n-1, m)*(g(m$2)[n]*x+1)))
        end:
    T:= n->(p->seq(coeff(p, x, i), i=0..degree(p)))(b(nops(g(n$2)), n)):
    seq(T(n), n=0..7);  # Alois P. Heinz, Aug 25 2019
  • Mathematica
    g[n_, i_] := g[n, i] = If[n==0 || i==1, {n!}, Join[Binomial[n, i]*#& /@ g[n - i, Min[n - i, i]], g[n, i - 1]]];
    b[n_, m_] := b[n, m] = If[n==0, 1, Expand[b[n-1, m]*(g[m, m][[n]]*x+1)]];
    T[n_] := CoefficientList[b[Length[g[n, n]], n], x];
    T /@ Range[0, 7] // Flatten (* Jean-François Alcover, Feb 18 2021, after Alois P. Heinz *)

Formula

Sum_{k=0..A000041(n)} (-1)^k * T(n,k) = 0.

A212806 Number of n X n matrices in which each row is a permutation of [1..n] and which contain no column rises.

Original entry on oeis.org

1, 3, 163, 271375, 21855093751, 128645361626874561, 78785944892341703819175577, 6795588328283070704898044776213094655, 107414633522643325764587104395687638119674465944431, 392471529081605251407320880492124164530148025908765037878553312273, 407934916447631403509359040563002566177814886353044858592046202746464825839911293037
Offset: 1

Views

Author

N. J. A. Sloane, May 27 2012

Keywords

Comments

A column rise in a matrix M = (m_{i,j}) is a value of j such that m_{i,j} < m_{i,j+1} for all i = 1..n.
From Petros Hadjicostas, Aug 26 2019: (Start)
Let R(m,n) := R(m,n,t=0) = A212855(m,n) for m,n >= 1, where R(m,n,t) = LHS of Eq. (6) of Abramson and Promislow (1978, p. 248).
Let P_n be the set of all lists b = (b_1, b_2,..., b_n) of integers b_i >= 0, i = 1,..., n, such that 1*b_1 + 2*b_2 + ... + n*b_n = n; i.e., P_n is the set all integer partitions of n. Then |P_n| = A000041(n) for n >= 0.
We have a(n) = R(n,n) = A212855(n,n) = Sum_{b in P_n} (-1)^(n - Sum_{j=1..n} b_j) * (b_1 + b_2 + ... + b_n)!/(b_1! * b_2! * ... * b_n!) * (n! / ((1!)^b_1 * (2!)^b_2 * ... * (n!)^b_n)^n.
(End)

Examples

			For n=2 the three matrices are [12/21], [21/12], [21/21] (but not [12/12]).
From _Petros Hadjicostas_, Aug 26 2019: (Start)
For example, when n = 3, the integer partitions of 3 are 3, 1+2, and 1+1+1, with corresponding (b_1, b_2, b_3) notation (0,0,1), (1,1,0), and (3,0,0). The corresponding multinomial coefficients are 3!/3! = 1, 3!/(1!*2!) = 3, and 3!/(1!*1!*1!) = 6, while the corresponding quantities (b_1 + b_2 + b_3)!/(b_1!*b_2!*b_3!) are 1, 2, and 1. The corresponding exponents of -1 (i.e., n - Sum_{j=1..n} b_j) are 3 - (0+0+1) = 2, 3 - (1+1+0) = 1, and 3 - (3+0+0) = 0.
It follows that a(n) = (-1)^2 * 1 * 1^3 + (-1)^1 * 2 * 3^3 + (-1)^0 * 1 * 6^3 = 163.
(End)
		

Crossrefs

Programs

  • Maple
    A212806 := proc(n) sum(z^k/k!^n, k=0..infinity);
    series(%^x, z=0, n+1): n!^n*coeff(%,z,n); add(abs(coeff(%,x,k)),k=0..n) end:
    seq(A212806(n), n=1..11); # Peter Luschny, May 27 2017
  • Mathematica
    a[n_] := Module[{s0, s1, s2}, s0 = Sum[z^k/k!^n, {k, 0, n}]; s1 =  Series[s0^x, {z, 0, n + 1}] // Normal; s2 = n!^n*Coefficient[s1, z, n]; Sum[Abs[Coefficient[s2, x, k]], {k, 0, n}]]; Array[a, 11] (* Jean-François Alcover, Feb 27 2018, after Peter Luschny *)
    T[n_, k_] := T[n, k] = If[k == 0, 1, -Sum[Binomial[k, j]^n*(-1)^j*T[n, k-j], {j, 1, k}]];
    a[n_] := T[n, n];
    Table[a[n], {n, 1, 12}] (* Jean-François Alcover, Apr 01 2024, after Alois P. Heinz in A212855 *)

Formula

Abramson and Promislow give a g.f. for R(m,n,t), the number of m X n matrices in which each row is a permutation of [1..n] and which contain exactly t column rises:
1 + Sum_{n>=1} Sum_{t=0..n-1} R(m,n,t) y^t x^n/(n!)^m = (y-1)/(y-f(x(y-1))) where f(x) = Sum_{i>=0} x^i/(i!)^m.

Extensions

Corrected by R. H. Hardin, May 28 2012

A376661 Frequency of the most common number among the multinomial coefficients n!/(x_1! * ... * x_k!) for all partitions (x_1, ..., x_k) of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4, 4, 4, 6, 6, 7, 8, 9, 11, 11, 13, 13, 14, 15, 16, 18, 19, 20, 23, 24, 26, 27, 30, 33, 37, 40, 43, 49, 52, 57, 64, 68, 76, 79, 87, 93, 99, 109, 116, 125, 135, 143, 157, 171, 191, 206, 223, 238, 254, 276, 291
Offset: 0

Views

Author

Pontus von Brömssen, Oct 02 2024

Keywords

Comments

Frequency of the most common number in row n of A036038 (for n >= 1) or A078760.
The sequence is nondecreasing, because a set of partitions of n-1 with a common multinomial coefficient can be extended to a set of partitions of n with a common multinomial coefficient by adding a unit part to each partition. It appears that a(n) > a(n-1) for n >= 28.
The sequence is unbounded. To see this, note that the sets of parts (1,1,1,4) and (2,2,3) of a partition can be exchanged without affecting the value of the multinomial coefficient, because 1+1+1+4 = 2+2+3 and 1!*1!*1!*4! = 2!*2!*3!. In particular, a((7*k)!/24^k) >= k+1 from the partitions 7*k = (3*j)*1 + j*4 + (2*(k-j))*2 + (k-j)*3 for 0 <= j <= k.

Examples

			For n = 7, the only number that appears more than once in row 7 of A036038 is 210, which appears twice: 210 = 7!/(2!*2!*3!) = 7!/(1!*1!*1!*4!). Hence, a(7) = 2.
		

Crossrefs

A210237 Triangle of distinct values M(n) of multinomial coefficients for partitions of n in increasing order of n and M(n).

Original entry on oeis.org

1, 1, 2, 1, 3, 6, 1, 4, 6, 12, 24, 1, 5, 10, 20, 30, 60, 120, 1, 6, 15, 20, 30, 60, 90, 120, 180, 360, 720, 1, 7, 21, 35, 42, 105, 140, 210, 420, 630, 840, 1260, 2520, 5040, 1, 8, 28, 56, 70, 168, 280, 336, 420, 560, 840, 1120, 1680, 2520, 3360, 5040, 6720
Offset: 1

Views

Author

Sergei Viznyuk, Mar 18 2012

Keywords

Comments

Differs from A036038 after a(37). To illustrate where the difference comes from, consider 4,1,1,1 and 3,2,2 are two different partitions of 7 having the same value of multinomial coefficient M(n)=n!/(m1!*m2!*...*mk!)=210.
There is no known formula for M(n) sequence, however the asymptotic behavior has been studied, see the paper by Andrews, Knopfmacher, and Zimmermann.
The number of terms per row (for each value of n starting with n=1) forms sequence A070289.

Examples

			Trianglebegins:
  1;
  1, 2;
  1, 3,  6;
  1, 4,  6, 12, 24;
  1, 5, 10, 20, 30,  60, 120;
  1, 6, 15, 20, 30,  60,  90, 120, 180, 360, 720;
  1, 7, 21, 35, 42, 105, 140, 210, 420, 630, 840, 1260, 2520, 5040;
  ...
Thus for n=4 (fourth row) the distinct values of multinomial coefficients are:
  4!/(4!) = 1
  4!/(3!1!) = 4
  4!/(2!2!) = 6
  4!/(2!1!1!) = 12
  4!/(1!1!1!1!) = 24
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i<2, {1},
          {seq(map(x-> x*i!^j, b(n-i*j, i-1))[], j=0..n/i)})
        end:
    T:= n-> sort([map(x-> n!/x, b(n, n))[]])[]:
    seq(T(n), n=1..10);  # Alois P. Heinz, Aug 13 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0 || i<2, {1}, Union[Flatten @ Table[(#*i!^j&) /@ b[n-i*j, i-1], {j, 0, n/i}]]]; T[n_] := Sort[Flatten[n!/#& /@ b[n, n]] ]; Table[T[n], {n, 1, 10}] // Flatten (* Jean-François Alcover, Feb 05 2017, after Alois P. Heinz *)

A215520 Number T(n,k) of distinct values of multinomial coefficients M(n;lambda), where lambda ranges over all partitions of n with largest part <= k; triangle T(n,k), n>=1, 1<=k<=n, read by rows.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 3, 4, 5, 1, 3, 5, 6, 7, 1, 4, 7, 9, 10, 11, 1, 4, 8, 10, 12, 13, 14, 1, 5, 10, 14, 17, 18, 19, 20, 1, 5, 12, 16, 21, 23, 25, 26, 27, 1, 6, 14, 20, 27, 29, 32, 34, 35, 36, 1, 6, 16, 22, 32, 35, 40, 43, 45, 46, 47, 1, 7, 19, 28, 40, 45, 52, 57, 60, 62, 63, 64
Offset: 1

Views

Author

Alois P. Heinz, Aug 14 2012

Keywords

Comments

Differs from A026820 first at T(7,4).
Differs from A213008 first at T(8,3).

Examples

			T(3,2) =  2 = |{3!/(2!*1!), 3!/(1!*1!*1!)}| = |{3, 6}|.
T(5,2) =  3 = |{30, 60, 120}|.
T(7,4) = 10 = |{35, 105, 140, 210, 420, 630, 840, 1260, 2520, 5040}|.
T(8,3) = 10 = |{560, 1120, 1680, 2520, 3360, 5040, 6720, 10080, 20160, 40320}|.
T(9,2) =  5 = |{22680, 45360, 90720, 181440, 362880}|.
Triangle T(n,k) begins:
  1;
  1,  2;
  1,  2,  3;
  1,  3,  4,  5;
  1,  3,  5,  6,  7;
  1,  4,  7,  9, 10, 11;
  1,  4,  8, 10, 12, 13, 14;
  1,  5, 10, 14, 17, 18, 19, 20;
  1,  5, 12, 16, 21, 23, 25, 26, 27;
  1,  6, 14, 20, 27, 29, 32, 34, 35, 36;
		

Crossrefs

Main diagonal gives: A070289.
Columns k=1-3 give: A000012, A008619 (for n>1), A001399 (for n>2).

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, {1}, `if`(i<1, {},
          {b(n, i-1)[], seq(map(x-> x*i!^j, b(n-i*j, i-1))[], j=1..n/i)}))
        end:
    T:= (n, k)-> nops(b(n, k)):
    seq(seq(T(n, k), k=1..n), n=1..14);
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, {1}, If[i < 1, {}, Join[b[n, i - 1], Table[ b[n - i*j, i - 1] *i!^j, {j, 1, n/i}] // Flatten]] // Union]; T[n_, k_] := Length[b[n, k]]; Table[Table[T[n, k], {k, 1, n}], {n, 1, 14}] // Flatten (* Jean-François Alcover, Jan 21 2015, after Alois P. Heinz *)
Showing 1-10 of 19 results. Next