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 15 results. Next

A000332 Binomial coefficient binomial(n,4) = n*(n-1)*(n-2)*(n-3)/24.

Original entry on oeis.org

0, 0, 0, 0, 1, 5, 15, 35, 70, 126, 210, 330, 495, 715, 1001, 1365, 1820, 2380, 3060, 3876, 4845, 5985, 7315, 8855, 10626, 12650, 14950, 17550, 20475, 23751, 27405, 31465, 35960, 40920, 46376, 52360, 58905, 66045, 73815, 82251, 91390, 101270, 111930, 123410
Offset: 0

Views

Author

Keywords

Comments

Number of intersection points of diagonals of convex n-gon where no more than two diagonals intersect at any point in the interior.
Also the number of equilateral triangles with vertices in an equilateral triangular array of points with n rows (offset 1), with any orientation. - Ignacio Larrosa Cañestro, Apr 09 2002. [See Les Reid link for proof. - N. J. A. Sloane, Apr 02 2016] [See Peter Kagey link for alternate proof. - Sameer Gauria, Jul 29 2025]
Start from cubane and attach amino acids according to the reaction scheme that describes the reaction between the active sites. See the hyperlink on chemistry. - Robert G. Wilson v, Aug 02 2002
For n>0, a(n) = (-1/8)*(coefficient of x in Zagier's polynomial P_(2n,n)). (Zagier's polynomials are used by PARI/GP for acceleration of alternating or positive series.)
Figurate numbers based on the 4-dimensional regular convex polytope called the regular 4-simplex, pentachoron, 5-cell, pentatope or 4-hypertetrahedron with Schlaefli symbol {3,3,3}. a(n)=((n*(n-1)*(n-2)*(n-3))/4!). - Michael J. Welch (mjw1(AT)ntlworld.com), Apr 01 2004, R. J. Mathar, Jul 07 2009
Maximal number of crossings that can be created by connecting n vertices with straight lines. - Cameron Redsell-Montgomerie (credsell(AT)uoguelph.ca), Jan 30 2007
If X is an n-set and Y a fixed (n-1)-subset of X then a(n) is equal to the number of 4-subsets of X intersecting Y. - Milan Janjic, Aug 15 2007
Product of four consecutive numbers divided by 24. - Artur Jasinski, Dec 02 2007
The only prime in this sequence is 5. - Artur Jasinski, Dec 02 2007
For strings consisting entirely of 0's and 1's, the number of distinct arrangements of four 1's such that 1's are not adjacent. The shortest possible string is 7 characters, of which there is only one solution: 1010101, corresponding to a(5). An eight-character string has 5 solutions, nine has 15, ten has 35 and so on, congruent to A000332. - Gil Broussard, Mar 19 2008
For a(n)>0, a(n) is pentagonal if and only if 3 does not divide n. All terms belong to the generalized pentagonal sequence (A001318). Cf. A000326, A145919, A145920. - Matthew Vandermast, Oct 28 2008
Nonzero terms = row sums of triangle A158824. - Gary W. Adamson, Mar 28 2009
Except for the 4 initial 0's, is equivalent to the partial sums of the tetrahedral numbers A000292. - Jeremy Cahill (jcahill(AT)inbox.com), Apr 15 2009
If the first 3 zeros are disregarded, that is, if one looks at binomial(n+3, 4) with n>=0, then it becomes a 'Matryoshka doll' sequence with alpha=0: seq(add(add(add(i,i=alpha..k),k=alpha..n),n=alpha..m),m=alpha..50). - Peter Luschny, Jul 14 2009
For n>=1, a(n) is the number of n-digit numbers the binary expansion of which contains two runs of 0's. - Vladimir Shevelev, Jul 30 2010
For n>0, a(n) is the number of crossing set partitions of {1,2,..,n} into n-2 blocks. - Peter Luschny, Apr 29 2011
The Kn3, Ca3 and Gi3 triangle sums of A139600 are related to the sequence given above, e.g., Gi3(n) = 2*A000332(n+3) - A000332(n+2) + 7*A000332(n+1). For the definitions of these triangle sums, see A180662. - Johannes W. Meijer, Apr 29 2011
For n > 3, a(n) is the hyper-Wiener index of the path graph on n-2 vertices. - Emeric Deutsch, Feb 15 2012
Except for the four initial zeros, number of all possible tetrahedra of any size, having the same orientation as the original regular tetrahedron, formed when intersecting the latter by planes parallel to its sides and dividing its edges into n equal parts. - V.J. Pohjola, Aug 31 2012
a(n+3) is the number of different ways to color the faces (or the vertices) of a regular tetrahedron with n colors if we count mirror images as the same.
a(n) = fallfac(n,4)/4! is also the number of independent components of an antisymmetric tensor of rank 4 and dimension n >= 1. Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Does not satisfy Benford's law [Ross, 2012] - N. J. A. Sloane, Feb 12 2017
Number of chiral pairs of colorings of the vertices (or faces) of a regular tetrahedron with n available colors. Chiral colorings come in pairs, each the reflection of the other. - Robert A. Russell, Jan 22 2020
From Mircea Dan Rus, Aug 26 2020: (Start)
a(n+3) is the number of lattice rectangles (squares included) in a staircase of order n; this is obtained by stacking n rows of consecutive unit lattice squares, aligned either to the left or to the right, which consist of 1, 2, 3, ..., n squares and which are stacked either in the increasing or in the decreasing order of their lengths. Below, there is a staircase or order 4 which contains a(7) = 35 rectangles. [See the Teofil Bogdan and Mircea Dan Rus link, problem 3, under A004320]
_
||
|||_
|||_|_
|||_|_|
(End)
a(n+4) is the number of strings of length n on an ordered alphabet of 5 letters where the characters in the word are in nondecreasing order. E.g., number of length-2 words is 15: aa,ab,ac,ad,ae,bb,bc,bd,be,cc,cd,ce,dd,de,ee. - Jim Nastos, Jan 18 2021
From Tom Copeland, Jun 07 2021: (Start)
Aside from the zeros, this is the fifth diagonal of the Pascal matrix A007318, the only nonvanishing diagonal (fifth) of the matrix representation IM = (A132440)^4/4! of the differential operator D^4/4!, when acting on the row vector of coefficients of an o.g.f., or power series.
M = e^{IM} is the matrix of coefficients of the Appell sequence p_n(x) = e^{D^4/4!} x^n = e^{b. D} x^n = (b. + x)^n = Sum_{k=0..n} binomial(n,k) b_n x^{n-k}, where the (b.)^n = b_n have the e.g.f. e^{b.t} = e^{t^4/4!}, which is that for A025036 aerated with triple zeros, the first column of M.
See A099174 and A000292 for analogous relationships for the third and fourth diagonals of the Pascal matrix. (End)
For integer m and positive integer r >= 3, the polynomial a(n) + a(n + m) + a(n + 2*m) + ... + a(n + r*m) in n has its zeros on the vertical line Re(n) = (3 - r*m)/2 in the complex plane. - Peter Bala, Jun 02 2024

Examples

			a(5) = 5 from the five independent components of an antisymmetric tensor A of rank 4 and dimension 5, namely A(1,2,3,4), A(1,2,3,5), A(1,2,4,5), A(1,3,4,5) and A(2,3,4,5). See the Dec 10 2015 comment. - _Wolfdieter Lang_, Dec 10 2015
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 74, Problem 8.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 70.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 7.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 294.
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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).
  • Charles W. Trigg, Mathematical Quickies, New York: Dover Publications, Inc., 1985, p. 53, #191.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 127.

Crossrefs

binomial(n, k): A161680 (k = 2), A000389 (k = 5), A000579 (k = 6), A000580 (k = 7), A000581 (k = 8), A000582 (k = 9).
Cf. A000217, A000292, A007318 (column k = 4).
Cf. A158824.
Cf. A006008 (Number of ways to color the faces (or vertices) of a regular tetrahedron with n colors when mirror images are counted as two).
Cf. A104712 (third column, k=4).
See A269747 for a 3-D analog.
Cf. A006008 (oriented), A006003 (achiral) tetrahedron colorings.
Row 3 of A325000, col. 4 of A007318.

Programs

  • GAP
    A000332 := List([1..10^2], n -> Binomial(n, 4)); # Muniru A Asiru, Oct 16 2017
    
  • Magma
    [Binomial(n,4): n in [0..50]]; // Vincenzo Librandi, Nov 23 2014
    
  • Maple
    A000332 := n->binomial(n,4); [seq(binomial(n,4), n=0..100)];
  • Mathematica
    Table[ Binomial[n, 4], {n, 0, 45} ] (* corrected by Harvey P. Dale, Aug 22 2011 *)
    Table[(n-4)(n-3)(n-2)(n-1)/24, {n, 100}] (* Artur Jasinski, Dec 02 2007 *)
    LinearRecurrence[{5,-10,10,-5,1}, {0,0,0,0,1}, 45] (* Harvey P. Dale, Aug 22 2011 *)
    CoefficientList[Series[x^4 / (1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Nov 23 2014 *)
  • PARI
    a(n)=binomial(n,4);
    
  • Python
    # Starts at a(3), i.e. computes n*(n+1)*(n+2)*(n+3)/24
    # which is more in line with A000217 and A000292.
    def A000332():
        x, y, z, u = 1, 1, 1, 1
        yield 0
        while True:
            yield x
            x, y, z, u = x + y + z + u + 1, y + z + u + 1, z + u + 1, u + 1
    a = A000332(); print([next(a) for i in range(41)]) # Peter Luschny, Aug 03 2019
    
  • Python
    print([n*(n-1)*(n-2)*(n-3)//24 for n in range(50)])
    # Gennady Eremin, Feb 06 2022

Formula

a(n) = n*(n-1)*(n-2)*(n-3)/24.
G.f.: x^4/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = n*a(n-1)/(n-4). - Benoit Cloitre, Apr 26 2003, R. J. Mathar, Jul 07 2009
a(n) = Sum_{k=1..n-3} Sum_{i=1..k} i*(i+1)/2. - Benoit Cloitre, Jun 15 2003
Convolution of natural numbers {1, 2, 3, 4, ...} and A000217, the triangular numbers {1, 3, 6, 10, ...}. - Jon Perry, Jun 25 2003
a(n) = A110555(n+1,4). - Reinhard Zumkeller, Jul 27 2005
a(n+1) = ((n^5-(n-1)^5) - (n^3-(n-1)^3))/24 - (n^5-(n-1)^5-1)/30; a(n) = A006322(n-2)-A006325(n-1). - Xavier Acloque, Oct 20 2003; R. J. Mathar, Jul 07 2009
a(4*n+2) = Pyr(n+4, 4*n+2) where the polygonal pyramidal numbers are defined for integers A>2 and B>=0 by Pyr(A, B) = B-th A-gonal pyramid number = ((A-2)*B^3 + 3*B^2 - (A-5)*B)/6; For all positive integers i and the pentagonal number function P(x) = x*(3*x-1)/2: a(3*i-2) = P(P(i)) and a(3*i-1) = P(P(i) + i); 1 + 24*a(n) = (n^2 + 3*n + 1)^2. - Jonathan Vos Post, Nov 15 2004
First differences of A000389(n). - Alexander Adamchuk, Dec 19 2004
For n > 3, the sum of the first n-2 tetrahedral numbers (A000292). - Martin Steven McCormick (mathseq(AT)wazer.net), Apr 06 2005 [Corrected by Doug Bell, Jun 25 2017]
Starting (1, 5, 15, 35, ...), = binomial transform of [1, 4, 6, 4, 1, 0, 0, 0, ...]. - Gary W. Adamson, Dec 28 2007
Sum_{n>=4} 1/a(n) = 4/3, from the Taylor expansion of (1-x)^3*log(1-x) in the limit x->1. - R. J. Mathar, Jan 27 2009
A034263(n) = (n+1)*a(n+4) - Sum_{i=0..n+3} a(i). Also A132458(n) = a(n)^2 - a(n-1)^2 for n>0. - Bruno Berselli, Dec 29 2010
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5); a(0)=0, a(1)=0, a(2)=0, a(3)=0, a(4)=1. - Harvey P. Dale, Aug 22 2011
a(n) = (binomial(n-1,2)^2 - binomial(n-1,2))/6. - Gary Detlefs, Nov 20 2011
a(n) = Sum_{k=1..n-2} Sum_{i=1..k} i*(n-k-2). - Wesley Ivan Hurt, Sep 25 2013
a(n) = (A000217(A000217(n-2) - 1))/3 = ((((n-2)^2 + (n-2))/2)^2 - (((n-2)^2 + (n-2))/2))/(2*3). - Raphie Frank, Jan 16 2014
Sum_{n>=0} a(n)/n! = e/24. Sum_{n>=3} a(n)/(n-3)! = 73*e/24. See A067764 regarding the second ratio. - Richard R. Forberg, Dec 26 2013
Sum_{n>=4} (-1)^(n+1)/a(n) = 32*log(2) - 64/3 = A242023 = 0.847376444589... . - Richard R. Forberg, Aug 11 2014
4/(Sum_{n>=m} 1/a(n)) = A027480(m-3), for m>=4. - Richard R. Forberg, Aug 12 2014
E.g.f.: x^4*exp(x)/24. - Robert Israel, Nov 23 2014
a(n+3) = C(n,1) + 3*C(n,2) + 3*C(n,3) + C(n,4). Each term indicates the number of ways to use n colors to color a tetrahedron with exactly 1, 2, 3, or 4 colors.
a(n) = A080852(1,n-4). - R. J. Mathar, Jul 28 2016
From Gary W. Adamson, Feb 06 2017: (Start)
G.f.: Starting (1, 5, 14, ...), x/(1-x)^5 can be written
as (x * r(x) * r(x^2) * r(x^4) * r(x^8) * ...) where r(x) = (1+x)^5;
as (x * r(x) * r(x^3) * r(x^9) * r(x^27) * ...) where r(x) = (1+x+x^2)^5;
as (x * r(x) * r(x^4) * r(x^16) * r(x^64) * ...) where r(x) = (1+x+x^2+x^3)^5;
... (as a conjectured infinite set). (End)
From Robert A. Russell, Jan 22 2020: (Start)
a(n) = A006008(n) - a(n+3) = (A006008(n) - A006003(n)) / 2 = a(n+3) - A006003(n).
a(n+3) = A006008(n) - a(n) = (A006008(n) + A006003(n)) / 2 = a(n) + A006003(n).
a(n) = A007318(n,4).
a(n+3) = A325000(3,n). (End)
Product_{n>=5} (1 - 1/a(n)) = cosh(sqrt(15)*Pi/2)/(100*Pi). - Amiram Eldar, Jan 21 2021

Extensions

Some formulas that referred to another offset corrected by R. J. Mathar, Jul 07 2009

A025035 Number of partitions of { 1, 2, ..., 3n } into sets of size 3.

Original entry on oeis.org

1, 1, 10, 280, 15400, 1401400, 190590400, 36212176000, 9161680528000, 2977546171600000, 1208883745669600000, 599606337852121600000, 356765771022012352000000, 250806337028474683456000000, 205661196363349240433920000000, 194555491759728381450488320000000
Offset: 0

Views

Author

Keywords

Comments

Row sums of A157703. - Johannes W. Meijer, Mar 07 2009
Number of bottom-row-increasing column-strict arrays of size 3 X n. - Ran Pan, Apr 10 2015
a(n) is the number of rooted semi-labeled or phylogenetic trees with n interior vertices and each interior vertex having out-degree 3. - Albert Alejandro Artiles Calix, Aug 12 2016

Examples

			G.f. = 1 + x + 10*x^2 + 280*x^3 + 15400*x^4 + 1401400*x^5 + ...
		

References

  • Erdos, Peter L., and L A. Szekely. "Applications of antilexicographic order. I. An enumerative theory of trees." Academic Press Inc. (1989): 488-96. Web. 4 July 2016.

Crossrefs

Column k=3 of A060540.

Programs

  • Magma
    [Factorial(3*n)/(Factorial(n)*6^n): n in [0..20]]; // Vincenzo Librandi, Apr 10 2015
  • Maple
    a := pochhammer(n+1, 2*n)/6^n: seq(a(n), n=0..15); # Peter Luschny, Nov 18 2019
  • Mathematica
    Select[Range[0, 39]! CoefficientList[Series[Exp[x^3/3!], {x, 0, 39}],  x], # > 0 &]  (* Geoffrey Critzer, Sep 24 2011 *)
    Table[(3 n)!/(n! (3!)^n), {n, 0, 15}] (* Michael De Vlieger, Aug 14 2016 *)
    a[ n_] := With[{m = 3 n}, If[ m < 0, 0, m! SeriesCoefficient[Exp[x^3/3!], {x, 0, m}]]]; (* Michael Somos, Nov 25 2016 *)
  • PARI
    {a(n) = if( n<0, 0, (3*n)! / n! / 6^n)}; /* Michael Somos, Mar 26 2003 */
    
  • PARI
    {a(n) = if( n<0, 0, prod( i=0, n-1, binomial( 3*n - 3*i, 3)) / n!)}; /* Michael Somos, Feb 15 2011 */
    
  • Sage
    [rising_factorial(n+1,2*n)/6^n for n in (0..15)] # Peter Luschny, Jun 26 2012
    

Formula

a(n) = (3*n)!/(n!*(3!)^n). - Christian G. Bower, Sep 01 1998
Integral representation as n-th moment of a positive function on the positive axis: a(n) = Integral_{x>=0} x^n*sqrt(2/(3*x))*BesselK(1/3, 2*sqrt(2*x)/3)/Pi dx, for n>=0. - Karol A. Penson, Oct 05 2005
E.g.f.: exp(x^3/3!) (with interpolated zeros). - Paul Barry, May 26 2003
a(n) = Product_{i=0..n-1} binomial(3*n-3*i,3) / n! (equivalent to Christian Bower formula). - Olivier Gérard, Feb 14 2011
2*a(n) - (3*n-1)*(3*n-2)*a(n-1) = 0. - R. J. Mathar, Dec 03 2012
a(n) ~ sqrt(3)*9^n*n^(2*n)/(2^n*exp(2*n)). - Ilya Gutkovskiy, Aug 12 2016
a(n) = Pochhammer(n + 1, 2*n)/6^n. - Peter Luschny, Nov 18 2019

A060540 Square array read by antidiagonals downwards: T(n,k) = (n*k)!/(k!^n*n!), (n>=1, k>=1), the number of ways of dividing nk labeled items into n unlabeled boxes with k items in each box.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 10, 15, 1, 1, 35, 280, 105, 1, 1, 126, 5775, 15400, 945, 1, 1, 462, 126126, 2627625, 1401400, 10395, 1, 1, 1716, 2858856, 488864376, 2546168625, 190590400, 135135, 1, 1, 6435, 66512160, 96197645544, 5194672859376, 4509264634875, 36212176000, 2027025, 1
Offset: 1

Views

Author

Henry Bottomley, Apr 02 2001

Keywords

Comments

The Copeland link gives the associations of this entry with the operator calculus of Appell Sheffer polynomials, the combinatorics of simple set partitions encoded in the Faa di Bruno formula for composition of analytic functions (formal Taylor series), the Pascal matrix, and the geometry of the n-dimensional simplices (hypertriangles, or hypertetrahedra). These, in turn, are related to simple instances of the application of the exponential formula / principle / schema giving the number of not-necessarily-connected objects composed from an ensemble of connected objects. - Tom Copeland, Jun 09 2021

Examples

			Array begins:
  1,   1,       1,          1,             1,                 1, ...
  1,   3,      10,         35,           126,               462, ...
  1,  15,     280,       5775,        126126,           2858856, ...
  1, 105,   15400,    2627625,     488864376,       96197645544, ...
  1, 945, 1401400, 2546168625, 5194672859376, 11423951396577720, ...
  ...
		

Crossrefs

Main diagonal is A057599.
Related to A057599, see also A096126 and A246048.
Cf. A060358, A361948 (includes row/col 0).
Cf. A000217, A000292, A000332, A000389, A000579, A000580, A007318, A036040, A099174, A133314, A132440, A135278 (associations in Copeland link).

Programs

  • Mathematica
    T[n_, k_] := (n*k)!/(k!^n*n!);
    Table[T[n-k+1, k], {n, 1, 10}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jun 29 2018 *)
  • PARI
    { i=0; for (m=1, 20, for (n=1, m, k=m - n + 1; write("b060540.txt", i++, " ", (n*k)!/(k!^n*n!))); ) } \\ Harry J. Smith, Jul 06 2009

Formula

T(n,k) = (n*k)!/(k!^n*n!) = T(n-1,k)*A060543(n,k) = A060538(n,k)/k!.
T(n,k) = Product_{j=2..n} binomial(j*k-1,k-1). - M. F. Hasler, Aug 22 2014

Extensions

Definition reworded by M. F. Hasler, Aug 23 2014

A260876 Number of m-shape set partitions, square array read by ascending antidiagonals, A(m,n) for m, n >= 0.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 4, 5, 5, 1, 1, 11, 31, 15, 7, 1, 1, 36, 365, 379, 52, 11, 1, 1, 127, 6271, 25323, 6556, 203, 15, 1, 1, 463, 129130, 3086331, 3068521, 150349, 877, 22, 1, 1, 1717, 2877421, 512251515, 3309362716, 583027547, 4373461, 4140, 30
Offset: 0

Views

Author

Peter Luschny, Aug 02 2015

Keywords

Comments

A set partition of m-shape is a partition of a set with cardinality m*n for some n >= 0 such that the sizes of the blocks are m times the parts of the integer partitions of n.
If m = 0, all possible sizes are zero. Thus the number of set partitions of 0-shape is the number of integer partitions of n (partition numbers A000041).
If m = 1, the set is {1, 2, ..., n} and the set of all possible sizes are the integer partitions of n. Thus the number of set partitions of 1-shape is the number of set partitions (Bell numbers A000110).
If m = 2, the set is {1, 2, ..., 2n} and the number of set partitions of 2-shape is the number of set partitions into even blocks A005046.
From Petros Hadjicostas, Aug 06 2019: (Start)
Irwin (1916) proved the following combinatorial result: Assume r_1, r_2, ..., r_n are positive integers and we have r_1*r_2*...*r_n objects. We divide them into r_1 classes of r_2*r_3*...*r_n objects each, then each class into r_2 subclasses of r_3*...*r_n objects each, and so on. We call each such classification, without reference to order, a "classification" par excellence. He proved that the total number of classifications is (r_1*r_2*...*r_n)!/( r1! * (r_2!)^(r_1) * (r_3!)^(r_1*r_2) * ... (r_n!)^(r_1*r_2*...*r_{n-1}) ).
Apparently, this problem appeared in Carmichael's "Theory of Numbers".
This result can definitely be used to prove some special cases of my conjecture below. (End)

Examples

			[ n ] [0  1   2       3        4           5              6]
[ m ] ------------------------------------------------------
[ 0 ] [1, 1,  2,      3,       5,          7,            11]  A000041
[ 1 ] [1, 1,  2,      5,      15,         52,           203]  A000110
[ 2 ] [1, 1,  4,     31,     379,       6556,        150349]  A005046
[ 3 ] [1, 1, 11,    365,   25323,    3068521,     583027547]  A291973
[ 4 ] [1, 1, 36,   6271, 3086331, 3309362716, 6626013560301]  A291975
        A260878,A309725, ...
For example the number of set partitions of {1,2,...,9} with sizes in [9], [6,3] and [3,3,3] is 1, 84 and 280 respectively. Thus A(3,3) = 365.
Formatted as a triangle:
[1]
[1, 1]
[1, 1,   2]
[1, 1,   2,    3]
[1, 1,   4,    5,     5]
[1, 1,  11,   31,    15,    7]
[1, 1,  36,  365,   379,   52,  11]
[1, 1, 127, 6271, 25323, 6556, 203, 15]
.
From _Peter Luschny_, Aug 14 2019: (Start)
For example consider the case n = 4. There are five integer partitions of 4:
  P = [[4], [3, 1], [2, 2], [2, 1, 1], [1, 1, 1, 1]]. The shapes are m times the parts of the integer partitions: S(m) = [[4m], [3m, m], [2m, 2m], [2m, m, m], [m, m, m, m]].
* In the case m = 1 we look at set partitions of {1, 2, 3, 4} with sizes in  [[4], [3, 1], [2, 2], [2, 1, 1], [1, 1, 1, 1]] which gives rise to [1, 4, 3, 6, 1] with sum 15.
* In the case m = 2 we look at set partitions of {1, 2, .., 8} with sizes in [[8], [6, 2], [4, 4], [4, 2, 2], [2, 2, 2, 2]] which gives rise to [1, 28, 35, 210, 105] with sum 379.
* In the case m = 0 we look at set partitions of {} with sizes in [[0], [0, 0], [0, 0], [0, 0, 0], [0, 0, 0, 0]] which gives rise to [1, 1, 1, 1, 1] with sum 5 (because the only partition of the empty set is the set that contains the empty set, thus from the definition T(0,4) = Sum_{S(0)} card({0}) = A000041(4) = 5).
If n runs through 0, 1, 2,... then the result is an irregular triangle in which the n-th row lists multinomials for partitions of [m*n] which have only parts which are multiples of m. These are the triangles A080575 (m = 1), A257490 (m = 2), A327003 (m = 3), A327004 (m = 4). In the case m = 0 the triangle is A000012 subdivided into rows of length A000041. See the cross references how this case integrates into the full picture.
(End)
		

Crossrefs

-----------------------------------------------------------------
[m] | multi- | sum of | main | by | comple- |
| nomials | rows | diagonal | size | mentary |
-----------------------------------------------------------------
Cf. A326996 (main diagonal), A260883 (ordered), A260875 (complementary).
Columns include A000012, A260878, A309725.

Programs

  • Maple
    A:= proc(m, n) option remember; `if`(m=0, combinat[numbpart](n),
          `if`(n=0, 1, add(binomial(m*n-1, m*k-1)*A(m, n-k), k=1..n)))
        end:
    seq(seq(A(d-n, n), n=0..d), d=0..10);  # Alois P. Heinz, Aug 14 2019
  • Mathematica
    A[m_, n_] := A[m, n] = If[m==0, PartitionsP[n], If[n==0, 1, Sum[Binomial[m n - 1, m k - 1] A[m, n - k], {k, 1, n}]]];
    Table[Table[A[d - n, n], {n, 0, d}], {d, 0, 10}] // Flatten (* Jean-François Alcover, Dec 06 2019, after Alois P. Heinz *)
  • SageMath
    def A260876(m, n):
        shapes = ([x*m for x in p] for p in Partitions(n))
        return sum(SetPartitions(sum(s), s).cardinality() for s in shapes)
    for m in (0..4): print([A260876(m,n) for n in (0..6)])

Formula

From Petros Hadjicostas, Aug 02 2019: (Start)
A(m, 2) = 1 + (1/2) * binomial(2*m, m) for m >= 1.
A(m, 3) = 1 + binomial(3*m, m) + (3*m)!/(6 * (m!)^3) for m >= 1.
A(m, 4) = (1/4!) * multinomial(4*m, [m, m, m, m]) + (1/2) * multinomial(4*m, [2*m, m, m]) + multinomial(4*m, [m, 3*m]) + (1/2) * multinomial(4*m, [2*m, 2*m]) + 1 for m >= 1.
Conjecture: For n >= 0, let P be the set of all possible lists (a_1,...,a_n) of nonnegative integers such that a_1*1 + a_2*2 + ... + a_n*n = n. Consider terms of the form multinomial(n*m, m*[1,..., 1,2,..., 2,..., n,..., n])/(a_1! * a_2! * ... * a_n!), where in the list [1,...,1,2,...,2,...,n,...,n] the number 1 occurs a_1 times, 2 occurs a_2 times, ..., and n occurs a_n times. (Here a_n = 0 or 1.) Summing these terms over P we get A(m, n) provided m >= 1. (End)
Conjecture for a recurrence: A(m, n) = Sum_{k = 0..n-1} binomial(m*n - 1, m*k) * A(m, k) with A(m, 0) = 1 for m >= 1 and n >= 0. (Unfortunately, the recurrence does not hold for m = 0.) - Petros Hadjicostas, Aug 12 2019

A291451 Triangle read by rows, expansion of e.g.f. exp(x*(exp(z)/3 + 2*exp(-z/2)* cos(z*sqrt(3)/2)/3 - 1)), nonzero coefficients of z.

Original entry on oeis.org

1, 0, 1, 0, 1, 10, 0, 1, 84, 280, 0, 1, 682, 9240, 15400, 0, 1, 5460, 260260, 1401400, 1401400, 0, 1, 43690, 7128576, 99379280, 285885600, 190590400, 0, 1, 349524, 193360720, 6600492080, 42549306800, 76045569600, 36212176000
Offset: 0

Views

Author

Peter Luschny, Sep 07 2017

Keywords

Examples

			Triangle starts:
[1]
[0, 1]
[0, 1,    10]
[0, 1,    84,     280]
[0, 1,   682,    9240,    15400]
[0, 1,  5460,  260260,  1401400,   1401400]
[0, 1, 43690, 7128576, 99379280, 285885600, 190590400]
		

Crossrefs

Cf. A048993 (m=1), A156289 (m=2), this seq. (m=3), A291452 (m=4).
Diagonal: A000012 (m=1), A001147 (m=2), A025035 (m=3), A025036 (m=4).
Row sums: A000110 (m=1), A005046 (m=2), A291973 (m=3), A291975 (m=4).
Alternating row sums: A000587 (m=1), A260884 (m=2), A291974 (m=3), A291976 (m=4).

Programs

  • Maple
    CL := (f, x) -> PolynomialTools:-CoefficientList(f, x):
    A291451_row := proc(n) exp(x*(exp(z)/3+2*exp(-z/2)*cos(z*sqrt(3)/2)/3-1)):
    series(%, z, 66): CL((3*n)!*coeff(series(%,z,3*(n+1)),z,3*n),x) end:
    for n from 0 to 7 do A291451_row(n) od;
    # Alternative:
    A291451row := proc(n) local P; P := proc(m, n) option remember;
    if n = 0 then 1 else add(binomial(m*n, m*k)*P(m, n-k)*x, k=1..n) fi end:
    CL(P(3, n), x); seq(%[k+1]/k!, k=0..n) end: # Peter Luschny, Sep 03 2018
  • Mathematica
    P[m_, n_] := P[m, n] = If[n == 0, 1, Sum[Binomial[m*n, m*k]*P[m, n - k]*x, {k, 1, n}]];
    row[n_] := Module[{cl = CoefficientList[P[3, n], x]}, Table[cl[[k + 1]]/k!, {k, 0, n}]];
    Table[row[n], {n, 0, 7}] // Flatten (* Jean-François Alcover, Jul 23 2019, after Peter Luschny *)

A291452 Triangle read by rows, expansion of e.g.f. exp(x*(cos(z) + cosh(z) - 2)/2), nonzero coefficients of z.

Original entry on oeis.org

1, 0, 1, 0, 1, 35, 0, 1, 495, 5775, 0, 1, 8255, 450450, 2627625, 0, 1, 130815, 35586525, 727476750, 2546168625, 0, 1, 2098175, 2941884000, 181262956875, 1932541986375, 4509264634875
Offset: 0

Views

Author

Peter Luschny, Sep 07 2017

Keywords

Examples

			Triangle starts:
[1]
[0, 1]
[0, 1,      35]
[0, 1,     495,       5775]
[0, 1,    8255,     450450,      2627625]
[0, 1,  130815,   35586525,    727476750,    2546168625]
[0, 1, 2098175, 2941884000, 181262956875, 1932541986375, 4509264634875]
		

Crossrefs

Cf. A048993 (m=1), A156289 (m=2), A291451 (m=3), this seq. (m=4).
Diagonal: A000012 (m=1), A001147 (m=2), A025035 (m=3), A025036 (m=4).
Row sums: A000110 (m=1), A005046 (m=2), A291973 (m=3), A291975 (m=4).
Alternating row sums: A000587 (m=1), A260884 (m=2), A291974 (m=3), A291976 (m=4).

Programs

  • Maple
    CL := (f,x) -> PolynomialTools:-CoefficientList(f,x):
    A291452_row := proc(n) exp(x*(cos(z)+cosh(z)-2)/2):
    series(%, z, 88): CL((4*n)!*coeff(series(%,z,4*(n+1)),z,4*n),x) end:
    for n from 0 to 7 do A291452_row(n) od;
    # Alternative:
    A291452row := proc(n) local P; P := proc(m, n) option remember;
    if n = 0 then 1 else add(binomial(m*n, m*k)*P(m, n-k)*x, k=1..n) fi end:
    CL(P(4, n), x); seq(%[k+1]/k!, k=0..n) end: # Peter Luschny, Sep 03 2018
  • Mathematica
    P[m_, n_] := P[m, n] = If[n == 0, 1, Sum[Binomial[m*n, m*k]*P[m, n - k]*x, {k, 1, n}]];
    row[n_] := Module[{cl = CoefficientList[P[4, n], x]}, Table[cl[[k + 1]]/k!, {k, 0, n}]];
    Table[row[n], {n, 0, 6}] // Flatten (* Jean-François Alcover, Jul 23 2019, after Peter Luschny *)

A334057 Triangle read by rows: T(n,k) is the number of configurations with exactly k polyomino matchings in a generalized game of memory played on the path of length 4n.

Original entry on oeis.org

1, 0, 1, 31, 3, 1, 5474, 288, 12, 1, 2554091, 72026, 1476, 31, 1, 2502018819, 43635625, 508610, 5505, 65, 1, 4456194509950, 52673302074, 394246455, 2559565, 16710, 120, 1, 13077453070386914, 111562882654972, 580589062179, 2504572910, 10288390, 43806, 203, 1
Offset: 0

Views

Author

Donovan Young, Apr 15 2020

Keywords

Comments

In this generalized game of memory n indistinguishable quadruples of matched cards are placed on the vertices of the path of length 4n. A polyomino is a quadruple on four adjacent vertices.
T(n,k) is the number of set partitions of {1..4n} into n sets of 4 with k of the sets being a contiguous set of elements. - Andrew Howroyd, Apr 16 2020

Examples

			The first few rows of T(n,k) are:
        1;
        0,     1;
       31,     3,    1;
     5474,   288,   12,  1;
  2554091, 72026, 1476, 31, 1;
  ...
For n=2 and k=1 the polyomino must start either on the second vertex of the path, the third, or the fourth, otherwise the remaining quadruple will also form a polyomino; thus T(2,1) = 3.
		

Crossrefs

Row sums are A025036.

Programs

  • Mathematica
    CoefficientList[Normal[Series[Sum[y^j*(4*j)!/24^j/j!/(1+y*(1-z))^(4*j+1),{j,0,20}],{y,0,20}]],{y,z}]
  • PARI
    T(n,k)={sum(j=0, n-k, (-1)^(n-j-k)*(n+3*j)!/(24^j*j!*(n-j-k)!*k!))} \\ Andrew Howroyd, Apr 16 2020

Formula

G.f.: Sum_{j>=0} (4*j)! * y^j / (j! * 24^j * (1+(1-z)*y)^(4*j+1)).
T(n,k) = Sum_{j=0..n-k} (-1)^(n-j-k)*(n+3*j)!/(24^j*j!*(n-j-k)!*k!). - Andrew Howroyd, Apr 16 2020

A200472 Triangle T(n,k) is the number of ways to assign n people to k unlabeled groups of equal size.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 1, 3, 0, 1, 1, 0, 0, 0, 1, 1, 10, 15, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 35, 0, 105, 0, 0, 0, 1, 1, 0, 280, 0, 0, 0, 0, 0, 1, 1, 126, 0, 0, 945, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 462, 5775, 15400, 0, 10395, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Dennis P. Walsh, Nov 18 2011

Keywords

Comments

If k is not a factor of n, T(n,k) = 0. If k is a factor of n, T(n,k) = (n!/k!)/(n/k)!^k. If n is a multiple of k, we may obtain T(n,k) by arranging all n people in an ordered line, which can be done in n! ways. Peel off the first n/k people for "group 1", the next n/k people for "group 2", ..., and the last n/k people for "group k". Since the k groups are actually unlabeled, we must divide n! by k! Also, since the ordering of the n/k people within each of the k groups is not of importance, we must now divide by (n/k)!^k. Therefore, T(n,k) = (n!/k!)/(n/k)!^k.
Also, T(2n,n) provide the sequence consisting of the products of consecutive odd integers.

Examples

			Triangle T(n,k) begins
1;
1,   1;
1,   0,    1;
1,   3,    0,     1;
1,   0,    0,     0,   1;
1,  10,   15,     0,   0,     1;
1,   0,    0,     0,   0,     0,  1;
1,  35,    0,   105,   0,     0,  0,  1;
1,   0,  280,     0,   0,     0,  0,  0,  1;
1, 126,    0,     0, 945,     0,  0,  0,  0,  1;
1,   0,    0,     0,   0,     0,  0,  0,  0,  0,  1;
1, 462, 5775, 15400,   0, 10395,  0,  0,  0,  0,  0,  1;
...
T(6,2) = 10 since there are 10 ways to assign 6 people (A,B,C,D,E,F) into 2 groups of size 3. The assignments are {A,B,C}|{D,E,F}, {A,B,D}|{C,E,F}, {A,B,E}|{C,D,F}, {A,B,F}|{C,D,E}, {A,C,D}|{B,E,F}, {A,C,E}|{B,D,F}, {A,C,F}|{B,D,E}, {B,C,D}|{A,E,F}, {B,C,E}|{A,D,F}, and {B,C,F}|{A,D,E}.
		

Crossrefs

T(2n,n) is A001147(n).
T(3n,n) is A025035(n).
T(4n,n) is A025036(n).
Row sums of T(n,k) provide A038041(n).
A200473 is A200472 with zeros removed.

Programs

  • Maple
    T:= (n, k)-> `if`(modp(n, k)=0, n!/(k!*((n/k)!)^k), 0):
    seq(seq(T(n, k), k=1..n), n=1..20);
  • Mathematica
    nn=11;s=Sum[Exp[y x^i/i!]-1,{i,1,nn}];Range[0,nn]!CoefficientList[Series[s,{x,0,nn}],{x,y}]//Grid  (* Geoffrey Critzer, Sep 15 2012 *)
  • PARI
    T(n,k) = if(n%k!=0, 0, (n!/k!)/((n/k)!)^k );
    for (n=1,15, for (k=1,n, print1(T(n,k),", "));print());
    /* Joerg Arndt, Sep 16 2012 */

Formula

For k that divide n, T(n,k) = (n!/k!)/((n/k)!)^k; otherwise, T(n,k) = 0.
E.g.f. when k is fixed: (1/k!) sum(j>=1, (x^j/j!)^k ).
E.g.f. for T(n*r,n): exp(x^r/r!).
T(2n,n) = (2n-1)!! = (2n-1)(2n-3)...(3)(1).

A326585 Coefficients of polynomials related to ordered set partitions. Triangle read by rows, T_{m}(n, k) for m = 4 and 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 36, 35, 0, 12046, 17820, 5775, 0, 16674906, 30263480, 16216200, 2627625, 0, 65544211366, 135417565890, 93516348900, 26189163000, 2546168625, 0, 588586227465426, 1334168329550300, 1083314031995250, 402794176785000, 69571511509500, 4509264634875
Offset: 0

Views

Author

Peter Luschny, Jul 21 2019

Keywords

Examples

			Triangle starts:
[0] [1]
[1] [0, 1]
[2] [0, 36, 35]
[3] [0, 12046, 17820, 5775]
[4] [0, 16674906, 30263480, 16216200, 2627625]
[5] [0, 65544211366, 135417565890, 93516348900, 26189163000, 2546168625]
[6] [0, 588586227465426, 1334168329550300, 1083314031995250, 402794176785000, 69571511509500, 4509264634875]
		

Crossrefs

Row sums A243665. Main diagonal A025036.
A129062 (m=1, associated with A131689), A326477 (m=2, associated with A241171), A326587 (m=3, associated with A278073), this sequence (m=4, associated with A278074).

Programs

Formula

T(n, k) = T_{4}(n, k) where T_{m}(n, k) is defined in A326477.

A360098 Square array read by antidiagonals upwards: T(n,k) is the number of ways of choosing nonnegative numbers for k n-sided dice, k >= 0, n >= 1, so that summing the faces can give any integer from 0 to n^k - 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 15, 1, 1, 1, 1, 1, 7, 1, 105, 1, 1, 1, 1, 1, 1, 71, 1, 945, 1, 1, 1, 1, 1, 10, 1, 1001, 1, 10395, 1, 1, 1, 1, 1, 3, 280, 1, 18089, 1, 135135, 1, 1, 1, 1, 1, 7, 15, 15400, 1, 398959, 1
Offset: 1

Views

Author

William P. Orrick, Jan 25 2023

Keywords

Comments

T(n,k) depends on n only through its prime signature. (See A118914.) For example, for fixed k, T(p*q^2,k) will be the same for any pair of distinct primes p and q. Hence we may define T(n,k) = R(s(n),k), where s(n) is the prime signature of n.
Also the number of Krasner factorizations of (x^(n^k)-1) / (x-1) into k polynomials each having n nonzero terms all with coefficient +1. (Krasner and Ranulac, 1937)

Examples

			There are 3 ways to assign numbers to two 4-sided dice:
 {{0, 1, 2, 3}, {0, 4, 8, 12}},
 {{0, 1, 8, 9}, {0, 2, 4,  6}},
 {{0, 1, 4, 5}, {0, 2, 8, 10}}.
The northwest corner of T(n,k) begins:
  1   1    1    1       1          1             1 ...     (s(1)  = {})
  1   1    1    1       1          1             1 ...     (s(2)  = {1})
  1   1    1    1       1          1             1 ...     (s(3)  = {1})
  1   1    3   15     105        945         10395 ...     (s(4)  = {2})
  1   1    1    1       1          1             1 ...     (s(5)  = {1})
  1   1    7   71    1001      18089        398959 ...     (s(6)  = {1,1})
  1   1    1    1       1          1             1 ...     (s(7)  = {1})
  1   1   10  280   15400    1401400     190590400 ...     (s(8)  = {3})
  1   1    3   15     105        945         10395 ...     (s(9)  = {2})
  1   1    7   71    1001      18089        398959 ...     (s(10) = {1,1})
  1   1    1    1       1          1             1 ...     (s(11) = {1})
  1   1   42  3660 614040  169200360   69444920160 ...     (s(12) = {1,2})
  1   1    1    1       1          1             1 ...     (s(13) = {1})
  1   1    7   71    1001      18089        398959 ...     (s(14) = {1,1})
  1   1    7   71    1001      18089        398959 ...     (s(15) = {1,1})
  1   1   35 5775 2627625 2546168625 4509264634875 ...     (s(16) = {4})
  ...
		

Crossrefs

For rows of index n = p^j, p prime, or equivalently, for rows of signature {j} we have T(p^2,k) = R({2},k) = A001147(k), T(p^3,k) = R({3},k) = A025035(k), T(p^4,k) = R({4},k) = A025036(k), and, generally, T(p^j,k) = R({j},k) = the k-th element of the j-th column of the square array A060540.
For n = p * q, p and q distinct primes, we have T(p*q,k) = R({1,1},k) = |A002119(k)|.
Column correspondences are T(n,2) = A273013(n) and T(n,3) = A131514(n).
Cf. A118914.

Programs

  • SageMath
    @cached_function
    def r(i,M):
        kminus1 = len(M)
        u = tuple([1 for j in range(kminus1)])
        if i > 1 and M == u:
            return(1)
        elif M != u:
            divList = divisors(i)[:-1]
            return(sum(r(M[j],tuple(sorted(M[:j]+tuple([d])+M[j+1:]))) for d in divList for j in range(kminus1)))
    def T(n,k):
        if n == 1 or k == 0:
            return(1)
        else:
            return(r(n,tuple([n for j in range(k-1)]))) / factorial(k-1)

Formula

Use M to denote a (k-1)-element multiset of positive integers. Let U denote the (k-1)-element multiset whose elements all equal 1 and let N denote the (k-1)-element multiset whose elements all equal n. For i in M, let M_{i,j} denote the result of replacing i with j in M. Then T(1,k) = T(n,0) = 1, while for n > 1 and k > 0 we have T(n,k) = r(n,N) / (k-1)! where r(i,M) is given by the recurrence
r(i,U) = 1 for i > 1,
r(i,M) = Sum_{m in M} Sum_{d|i,d
Showing 1-10 of 15 results. Next