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

A005651 Sum of multinomial coefficients (n_1+n_2+...)!/(n_1!*n_2!*...) where (n_1, n_2, ...) runs over all integer partitions of n.

Original entry on oeis.org

1, 1, 3, 10, 47, 246, 1602, 11481, 95503, 871030, 8879558, 98329551, 1191578522, 15543026747, 218668538441, 3285749117475, 52700813279423, 896697825211142, 16160442591627990, 307183340680888755, 6147451460222703502, 129125045333789172825, 2841626597871149750951
Offset: 0

Views

Author

Keywords

Comments

This is the total number of hierarchies of n labeled elements arranged on 1 to n levels. A distribution of elements onto levels is "hierarchical" if a level l+1 contains <= elements than level l. Thus for n=4 the arrangement {1,2}:{3}{4} is not allowed. See also A140585. Examples: Let the colon ":" separate two consecutive levels l and l+1. Then n=2 --> 3: {1}{2}, {1}:{2}, {2}:{1}, n=3 --> 10: {1}{2}{3}, {1}{2}:{3}, {3}{1}:{2}, {2}{3}:{1}, {1}:{2}:{3}, {3}:{1}:{2}, {2}:{3}:{1}, {1}:{3}:{2}, {2}:{1}:{3}, {3}:{2}:{1}. - Thomas Wieder, May 17 2008
n identical objects are painted by dipping them into a long row of cans of paint of distinct colors. Begining with the first can and not skipping any cans k, 1<=k<=n, objects are dipped (painted) and not more objects are dipped into any subsequent can than were dipped into the previous can. The painted objects are then linearly ordered. - Geoffrey Critzer, Jun 08 2009
a(n) is the number of partitions of n where each part i is marked with a word of length i over an n-ary alphabet whose letters appear in alphabetical order and all n letters occur exactly once in the partition. a(3) = 10: 3abc, 2ab1c, 2ac1b, 2bc1a, 1a1b1c, 1a1c1b, 1b1a1c, 1b1c1a, 1c1a1b, 1c1b1a. - Alois P. Heinz, Aug 30 2015
Also the number of ordered set partitions of {1,...,n} with weakly decreasing block sizes. - Gus Wiseman, Sep 03 2018
The parity of a(n) is that of A000110(A000120(n)), so a(n) is even if and only if A000120(n) == 2 (mod 3). - Álvar Ibeas, Aug 11 2020

Examples

			For n=3, say the first three cans in the row contain red, white, and blue paint respectively. The objects can be painted r,r,r or r,r,w or r,w,b and then linearly ordered in 1 + 3 + 6 = 10 ways. - _Geoffrey Critzer_, Jun 08 2009
From _Gus Wiseman_, Sep 03 2018: (Start)
The a(3) = 10 ordered set partitions with weakly decreasing block sizes:
  {{1},{2},{3}}
  {{1},{3},{2}}
  {{2},{1},{3}}
  {{2},{3},{1}}
  {{3},{1},{2}}
  {{3},{2},{1}}
  {{2,3},{1}}
  {{1,2},{3}}
  {{1,3},{2}}
  {{1,2,3}}
(End)
		

References

  • Abramowitz and Stegun, Handbook, p. 831, column labeled "M_1".
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 126.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Main diagonal of: A226873, A261719, A309973.
Row sums of: A226874, A262071, A327803.
Column k=1 of A309951.
Column k=0 of A327801.

Programs

  • Maple
    A005651b := proc(k) add( d/(d!)^(k/d),d=numtheory[divisors](k)) ; end proc:
    A005651 := proc(n) option remember; local k ; if n <= 1 then 1; else (n-1)!*add(A005651b(k)*procname(n-k)/(n-k)!, k=1..n) ; end if; end proc:
    seq(A005651(k), k=0..10) ; # R. J. Mathar, Jan 03 2011
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0 or i=1, n!,
          b(n, i-1) +binomial(n, i)*b(n-i, min(n-i, i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..25);  # Alois P. Heinz, Aug 29 2015, Dec 12 2016
  • Mathematica
    Table[Total[n!/Map[Function[n, Apply[Times, n! ]], IntegerPartitions[n]]], {n, 0, 20}] (* Geoffrey Critzer, Jun 08 2009 *)
    Table[Total[Apply[Multinomial, IntegerPartitions[n], {1}]], {n, 0, 20}] (* Jean-François Alcover and Olivier Gérard, Sep 11 2014 *)
    b[n_, i_, t_] := b[n, i, t] = If[t==1, 1/n!, Sum[b[n-j, j, t-1]/j!, {j, i, n/t}]]; a[n_] := If[n==0, 1, n!*b[n, 0, n]]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Nov 20 2015, after Alois P. Heinz *)
  • Maxima
    a(m,n):=if n=m then 1 else sum(binomial(n,k)*a(k,n-k),k,m,(n/2))+1;
    makelist(a(1,n),n,0,17); /* Vladimir Kruchinin, Sep 06 2014 */
    
  • PARI
    a(n)=my(N=n!,s);forpart(x=n,s+=N/prod(i=1,#x,x[i]!));s \\ Charles R Greathouse IV, May 01 2015
    
  • PARI
    { my(n=25); Vec(serlaplace(prod(k=1, n, 1/(1-x^k/k!) + O(x*x^n)))) } \\ Andrew Howroyd, Dec 20 2017

Formula

E.g.f.: 1 / Product (1 - x^k/k!).
a(n) = Sum_{k=1..n} (n-1)!/(n-k)!*b(k)*a(n-k), where b(k) = Sum_{d divides k} d*d!^(-k/d). - Vladeta Jovovic, Oct 14 2002
a(n) ~ c * n!, where c = Product_{k>=2} 1/(1-1/k!) = A247551 = 2.52947747207915264... . - Vaclav Kotesovec, May 09 2014
a(n) = S(n,1), where S(n,m) = sum(k=m..n/2 , binomial(n,k)*S(n-k,k))+1, S(n,n)=1, S(n,m)=0 for nVladimir Kruchinin, Sep 06 2014
E.g.f.: exp(Sum_{k>=1} Sum_{j>=1} x^(j*k)/(k*(j!)^k)). - Ilya Gutkovskiy, Jun 18 2018

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Mar 29 2003

A000275 Coefficients of a Bessel function (reciprocal of J_0(z)); also pairs of permutations with rise/rise forbidden.

Original entry on oeis.org

1, 1, 3, 19, 211, 3651, 90921, 3081513, 136407699, 7642177651, 528579161353, 44237263696473, 4405990782649369, 515018848029036937, 69818743428262376523, 10865441556038181291819, 1923889742567310611949459, 384565973956329859109177427, 86180438505835750284241676121
Offset: 0

Views

Author

Keywords

Comments

a(n) has the Lucas property, namely a(n) is congruent to a(n_0)a(n_1)...a(n_k) modulo p for any prime p where n_0,n_1,... are the base p digits of n. (Carlitz via McIntosh)

Examples

			From _Peter Bala_, Aug 08 2011: (Start)
a(3) = 19: The 19 pairs of permutations in the group S_3 x S_3 with no common rises correspond to the zero entries in the table below.
  ======================================
   Number of common rises in S_3 x S_3
  ======================================
     | 123   132   213   231   312   321
  ======================================
  123|  2     1     1     1     1     0
  132|  1     1     0     1     0     0
  213|  1     0     1     0     1     0
  231|  1     1     0     1     0     0
  312|  1     0     1     0     1     0
  321|  0     0     0     0     0     0
(End)
G.f. = 1 + x + 3*x^2 + 19*x^3 + 211*x^4 + 3651*x^5 + 90921*x^6 + ...
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Row 2 of A212855.
Cf. A055133 (absolute value of column 0 of triangle), A192721 (column 1), A115368.
Column k=1 of A340986.

Programs

  • Maple
    A000275 := proc(n) sum(z^k/k!^2, k = 0..infinity);
    series(%^x, z=0, n+1): n!^2*coeff(%,z,n); add(abs(coeff(%,x,k)), k=0..n) end:
    seq(A000275(n), n=0..17); # Peter Luschny, May 27 2017
  • Mathematica
    a[0] = 1; a[n_] := a[n] = Sum[(-1)^(r+n+1)*Binomial[n, r]^2 a[r], {r, 0, n-1}]; Table[a[n], {n, 0, 17}] (* Jean-François Alcover, Aug 05 2013 *)
    CoefficientList[Series[1/BesselJ[0,Sqrt[4*x]], {x, 0, 20}], x]* Range[0, 20]!^2 (* Vaclav Kotesovec, Mar 02 2014 *)
    a[ n_] := If[ n < 0, 0, (n! 2^n)^2 SeriesCoefficient[ 1 / BesselJ[ 0, x], {x, 0, 2 n}]]; (* Michael Somos, Aug 20 2015 *)
  • PARI
    {a(n) = if( n<0, 0, n!^2 * 4^n * polcoeff( 1 / besselj(0, x + x * O(x^(2*n))), 2*n))}; /* Michael Somos, May 17 2004 */

Formula

a(n) = Sum_{r=0..n-1} (-1)^(r+n+1) binomial(n, r)^2 a(r), if n > 0.
Sum_{n>=0} a(n) * x^n / n!^2 = 1 / J_0(sqrt(4*x)). - _Michael Somos, May 17 2004
From Peter Bala, Aug 08 2011: (Start)
Conjectural formula: 1 = Sum_{n>=0} a(n)*x^n*Sum_{k>=0} binomial(n+k,k)^2*(-x)^k.
Apart from the initial term, first column of A192721. (End)
E.g.f.: 1/J_0(sqrt(4*x)) = 1 + x/Q(0), where Q(k) = (k+1)^2 - x + (k+1)^2*x/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Dec 06 2013
a(n) ~ c * (n!)^2 / r^n, where r = (1/4)*BesselJZero[0,1]^2 = 1.4457964907366961302939989396139517587678604516... and c = 1/(sqrt(r) * BesselJ(1, 2*sqrt(r))) = 1.60197469692804662664846689139151227422675123376219... - Vaclav Kotesovec, Mar 02 2014, updated Apr 01 2018

Extensions

More terms from Christian G. Bower, Apr 25 2000

A070289 Number of distinct values of multinomial coefficients ( n / (p1, p2, p3, ...) ) where (p1, p2, p3, ...) runs over all partitions of n.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 14, 20, 27, 36, 47, 64, 79, 102, 125, 157, 193, 243, 296, 366, 441, 538, 639, 773, 911, 1092, 1294, 1532, 1799, 2131, 2475, 2901, 3369, 3935, 4554, 5292, 6084, 7033, 8087, 9292, 10617, 12198, 13880, 15874, 18039, 20541, 23263, 26414, 29838
Offset: 0

Views

Author

Naohiro Nomoto, May 12 2002

Keywords

Crossrefs

Programs

  • Maple
    b:= proc(n,i) option remember;
          if n=0 then {1} elif i<1 then {} else {b(n, i-1)[],
             seq(map(x-> x*i!^j, b(n-i*j, i-1))[], j=1..n/i)} fi
        end:
    a:= n-> nops(b(n, n)):
    seq(a(n), n=0..50);  # Alois P. Heinz, Aug 14 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, {1}, If[i<1, {}, Union[Join[b[n, i-1], Flatten[ Table[Function[{x}, x*i!^j] /@ b[n-i*j, i-1], {j, 1, n/i}]]]]]]; a[n_] := Length[b[n, n]]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Mar 23 2015, after Alois P. Heinz *)
  • Sage
    def A070289(n):
        P = Partitions(n)
        M = set(multinomial(list(x)) for x in P)
        return len(M)
    [A070289(n) for n in range(20)]
    # Joerg Arndt, Aug 14 2012

Formula

a(n) = A215520(n,n) = A215521(2*n,n). - Alois P. Heinz, Nov 08 2012

Extensions

Terms a(n) for n >= 45 corrected by Joerg Arndt and Alois P. Heinz, Aug 14 2012

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

Original entry on oeis.org

1, 1, 7, 163, 8983, 966751, 179781181, 53090086057, 23402291822743, 14687940716402023, 12645496977257273257, 14490686095184389113277, 21557960797148733086439949, 40776761007750226749220637461, 96332276574683758035941025907591
Offset: 0

Views

Author

R. H. Hardin, May 28 2012

Keywords

Examples

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

Crossrefs

Programs

  • Maple
    A212856 := proc(n) sum(z^k/k!^3, k = 0..infinity);
    series(%^x, z=0, n+1): n!^3*coeff(%,z,n); add(abs(coeff(%,x,k)), k=0..n) end:
    seq(A212856(n), n=0..14); # Peter Luschny, May 27 2017
    # second Maple program:
    a:= proc(n) option remember; `if`(n=0, 1, -add(
          binomial(n, j)^3*(-1)^j*a(n-j), j=1..n))
        end:
    seq(a(n), n=0..15);  # Alois P. Heinz, Apr 26 2020
  • Mathematica
    f[0] = 1; f[n_] := f[n] = Sum[(-1)^(n+k+1)*f[k]*Binomial[n, k]^2/(n-k)!, {k, 0, n-1}]; a[n_] := f[n]*n!; Array[a, 14] (* Jean-François Alcover, Feb 27 2018, after Daniel Suteu *)

Formula

a(n) = f(n) * n!, where f(0) = 1, f(n) = Sum_{k=0..n-1} (-1)^(n+k+1) * f(k) * binomial(n, k)^2 / (n-k)!. - Daniel Suteu, Feb 23 2018
a(n) = (n!)^3 * [x^n] 1 / (1 + Sum_{k>=1} (-x)^k / (k!)^3). - Seiichi Manyama, Jul 18 2020
a(n) ~ c * n!^3 / r^n, where r = 1.16151549806386358435938834554462085598002... is the root of the equation HypergeometricPFQ[{}, {1, 1}, -r] = 0 and c = 1.182760720067731330743886867947078139186402925891650811631774628... - Vaclav Kotesovec, Sep 16 2020

Extensions

a(0)=1 prepended by Alois P. Heinz, Apr 26 2020

A212850 Number of n X 3 arrays with rows being permutations of 0..2 and no column j greater than column j-1 in all rows.

Original entry on oeis.org

1, 19, 163, 1135, 7291, 45199, 275563, 1666495, 10038331, 60348079, 362442763, 2175719455, 13057505371, 78354598159, 470156286763, 2821023814015, 16926401164411, 101559181827439, 609357415487563, 3656151466494175
Offset: 1

Views

Author

R. H. Hardin, May 28 2012

Keywords

Comments

From Petros Hadjicostas, Aug 25 2019: (Start)
Both formulas below follow from the theory in the documentation of array A309951. We have Sum_{s = 0..A000041(3)} (-1)^s * A309951(3,s) * a(n-s) = 0, i.e., a(n) - 10*a(n-1) - 27*a(n-2) + 18*a(n-3) = 0 for n >= 4. This is a consequence of Eq. (6) on p. 248 of Abramson and Promislow (1978), where we let t=0 in the equation.
In the explicit formula by Vaclav Kotesovec below, a(n) = 6^n - 2*3^n + 1^n, the numbers 1, 3, 6 (that are raised to the n-th power) are the multinomial coefficients of the A000041(3) = 3 integer partitions of 3: 1 = 3!/3!, 3 = 3!/(1!2!), 6 = 3!/(1!1!1!).
(End)

Examples

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

Crossrefs

Formula

Empirical: a(n) = 10*a(n-1) - 27*a(n-2) + 18*a(n-3).
Explicit formula: a(n) = 6^n - 2*3^n + 1. - Vaclav Kotesovec, May 31 2012

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

Original entry on oeis.org

1, 211, 8983, 271375, 7225951, 182199871, 4479288703, 108787179775, 2626338801151, 63217691436031, 1519452489242623, 36493601345048575, 876167372044132351, 21031868446675976191, 504811062363654815743, 12116020140998121291775, 290791139166323355287551
Offset: 1

Views

Author

R. H. Hardin, May 28 2012

Keywords

Comments

Column 4 of A212855.
From Petros Hadjicostas, Aug 25 2019: (Start)
All formulas below follow from the theory in the documentation of array A309951.
We have Sum_{s = 0..A000041(4)} (-1)^s * A309951(4,s) * a(n-s) = 0, i.e., a(n) - 47*a(n-1) + 718*a(n-2) - 4416*a(n-3) + 10656*a(n-4) - 6912*a(n-5) = 0 for n >= 6. This is a consequence of Eq. (6) on p. 248 of Abramson and Promislow (1978).
Note that in R. J. Mathar's formula a(n) = 24^n + 6^n - 3*12^n + 2*4^n - 1^n, the numbers 1, 4, 12, 6, and 24 (that are raised to the n-th power) are the multinomial coefficients of the A000041(4) = 5 integer partitions of 4: 4!/4! = 1, 4!/(1!3!) = 4, 12 = 4!/(1!1!2!), 6 = 4!/(2!2!), 24 = 4!/(1!1!1!1!).
Note also that these numbers appear also in the denominator of the Colin Barker's g.f.: (1 - x)*(1 - 4*x)*(1 - 6*x)*(1 - 12*x)*(1 - 24*x) = 1 - 47*x + 718*x^2 - 4416*x^3 + 10656*x^4 - 6912*x^5. (End)

Examples

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

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, 4];
    Table[a[n], {n, 1, 15}] (* Jean-François Alcover, Apr 01 2024, after Alois P. Heinz in A212855 *)

Formula

Empirical: a(n) = 47*a(n-1) - 718*a(n-2) + 4416*a(n-3) - 10656*a(n-4) + 6912*a(n-5).
Empirical: a(n) = 24^n + 6^n - 3*12^n + 2*4^n - 1. R. J. Mathar, Jun 25 2012
Empirical g.f.: x*(1 + 164*x - 216*x^2 - 3744*x^3) / ((1 - x)*(1 - 4*x)*(1 - 6*x)*(1 - 12*x)*(1 - 24*x)). - Colin Barker, Jul 21 2018

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