cp's OEIS Frontend

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

Showing 1-4 of 4 results.

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

A334058 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 5n.

Original entry on oeis.org

1, 0, 1, 121, 4, 1, 124760, 1347, 18, 1, 486854621, 2001548, 8154, 52, 1, 5184423824705, 10231953233, 17045774, 35542, 121, 1, 123243726413573515, 134835947255262, 112619668659, 102416812, 124881, 246, 1, 5717986519188343198259, 3821094862609800013, 1820735766620673, 863827126967, 486979381, 375627, 455, 1
Offset: 0

Views

Author

Donovan Young, Apr 15 2020

Keywords

Comments

In this generalized game of memory n indistinguishable quintuples of matched cards are placed on the vertices of the path of length 5n. A polyomino is a quintuple on five adjacent vertices.
T(n,k) is the number of set partitions of {1..5n} into n sets of 5 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;
        121,       4,    1;
     124760,    1347,   18,  1;
  486854621, 2001548, 8154, 52, 1;
  ...
For n=2 and k=1 the polyomino must start either on the second, third, fourth, or fifth vertex of the path, otherwise the remaining quintuple will also form a polyomino; thus T(2,1) = 4.
		

Crossrefs

Row sums are A025037.

Programs

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

Formula

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

A361948 Array read by ascending antidiagonals. A(n, k) = Product_{j=0..k-1} binomial((j + 1)*n - 1, n - 1) if n >= 1, and A(0, k) = 1 for all k.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 10, 15, 1, 1, 1, 1, 35, 280, 105, 1, 1, 1, 1, 126, 5775, 15400, 945, 1, 1, 1, 1, 462, 126126, 2627625, 1401400, 10395, 1, 1, 1, 1, 1716, 2858856, 488864376, 2546168625, 190590400, 135135, 1, 1
Offset: 0

Views

Author

Peter Luschny, Apr 13 2023

Keywords

Comments

Row n gives the leading coefficients of the set partition polynomials of type n. The sequence of these polynomial sequences starts: A097805, A048993, A156289, A291451, A291452, ...

Examples

			Array A(n, k) starts:
  [0] 1, 1,   1,       1,           1,                 1, ...
  [1] 1, 1,   1,       1,           1,                 1, ...
  [2] 1, 1,   3,      15,         105,               945, ...  A001147
  [3] 1, 1,  10,     280,       15400,           1401400, ...  A025035
  [4] 1, 1,  35,    5775,     2627625,        2546168625, ...  A025036
  [5] 1, 1, 126,  126126,   488864376,     5194672859376, ...  A025037
  [6] 1, 1, 462, 2858856, 96197645544, 11423951396577720, ...  A025038
.
Triangle A(n-k, k) starts:
  [0] 1;
  [1] 1, 1;
  [2] 1, 1,  1;
  [3] 1, 1,  1,   1;
  [4] 1, 1,  3,   1,   1;
  [5] 1, 1, 10,  15,   1, 1;
  [6] 1, 1, 35, 280, 105, 1, 1;
		

Crossrefs

Cf. A060540 (subarray), A370407 (antidiagonal sums, row sums).
Cf. A001147 (row 2), A025035 (row 3), A025036 (row 4), A025037 (row 5), A025038 (row 6), A025039 (row 7), A025040 (row 8), A025041 (row 9).
Cf. A088218 (column 2), A060542 (column 3), A082368 (column 4), A322252 (column 5), A057599 (main diagonal).

Programs

  • Maple
    A := (n, k) -> mul(binomial((j + 1)*n - 1, n - 1), j = 0..k-1):
    seq(seq(A(n-k, k), k = 0..n), n = 0..9);
    # Alternative, using recursion:
    A := proc(n, k) local P; P := proc(n, k) option remember;
    if n = 0 then return x^k*k! fi; if k = 0 then 1 else add(binomial(n*k, n*j)*
    P(n,k-j)*x, j=1..k) fi end: coeff(P(n, k), x, k) / k! end:
    seq(print(seq(A(n, k), k = 0..5)), n = 0..6);
    # Alternative, using exponential generating function:
    egf := n -> ifelse(n=0, 1, exp(x^n/n!)): ser := n -> series(egf(n), x, 8*n):
    row := n -> local k; seq((n*k)!*coeff(ser(n), x, n*k), k = 0..6):
    for n from 0 to 6 do [n], row(n) od;  # Peter Luschny, Aug 15 2024
  • Mathematica
    A[n_, k_] := Product[Binomial[n (j + 1) - 1, n - 1], {j, 0, k - 1}]; Table[A[n - k, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Michael De Vlieger, Apr 13 2023 *)
  • SageMath
    def Arow(n, size):
        if n == 0: return [1] * size
        return [prod(binomial((j + 1)*n - 1, n - 1) for j in range(k)) for k in range(size)]
    for n in range(7): print(Arow(n, 7))
    # Alternative, using exponential generating function:
    def SetPolyLeadCoeff(m, n):
        x, z = var("x, z")
        if m == 0: return 1
        w = exp(2 * pi * I / m)
        o = sum(exp(z * w ** k) for k in range(m)) / m
        t = exp(x * (o - 1)).taylor(z, 0, m*n)
        p = factorial(m*n) * t.coefficient(z, m*n)
        return p.leading_coefficient(x)
    for m in range(7):
        print([SetPolyLeadCoeff(m, k) for k in range(6)])

Formula

A(n, k) = (1/k!) * [x^k] P(n, k), where P(n, k) = k!*x^k if n = 0 and otherwise 1 if k = 0 and otherwise Sum_{j=1..k} binomial(n*k, n*j)*P(n, k-j)*x.
A(n, k) = (n*k)!*[x^(n*k)] exp(x^n/n!) for n >= 1. - Peter Luschny, Aug 15 2024

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

Original entry on oeis.org

1, 0, 1, 0, 127, 126, 0, 255256, 381381, 126126, 0, 2979852651, 5447453786, 2956465512, 488864376, 0, 127156445503275, 264284637872750, 184292523727620, 52359004217520, 5194672859376
Offset: 0

Views

Author

Peter Luschny, Jul 21 2019

Keywords

Examples

			Triangle starts:
[0] [1]
[1] [0, 1]
[2] [0, 127, 126]
[3] [0, 255256, 381381, 126126]
[4] [0, 2979852651, 5447453786, 2956465512, 488864376]
[5] [0, 127156445503275, 264284637872750, 184292523727620, 52359004217520, 5194672859376]
[6] [0, 15160169962750251082, 34544220081315967665, 28276496764200664980, 10634436034307385300, 1865368063755476280, 123378675083039376]
		

Crossrefs

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

Programs

Formula

T(n, k) = T_{5}(n, k) where T_{m}(n, k) is defined in A326477.
Showing 1-4 of 4 results.