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

A049311 Number of (0,1) matrices with n ones and no zero rows or columns, up to row and column permutations.

Original entry on oeis.org

1, 3, 6, 16, 34, 90, 211, 558, 1430, 3908, 10725, 30825, 90156, 273234, 848355, 2714399, 8909057, 30042866, 103859678, 368075596, 1335537312, 4958599228, 18820993913, 72980867400, 288885080660, 1166541823566, 4802259167367, 20141650236664
Offset: 1

Views

Author

Keywords

Comments

Also the number of bipartite graphs with n edges, no isolated vertices and a distinguished bipartite block, up to isomorphism.
The EULERi transform (A056156) is also interesting.
a(n) is also the number of non-isomorphic set multipartitions (multisets of sets) of weight n. - Gus Wiseman, Mar 17 2017

Examples

			E.g. a(2) = 3: two ones in same row, two ones in same column, or neither.
a(3) = 6 is coefficient of x^3 in (1/36)*((1 + x)^9 + 6*(1 + x)^3*(1 + x^2)^3 + 8*(1 + x^3)^3 + 9*(1 + x)*(1 + x^2)^4 + 12*(1 + x^3)*(1 + x^6))=1 + x + 3*x^2 + 6*x^3 + 7*x^4 + 7*x^5 + 6*x^6 + 3*x^7 + x^8 + x^9.
There are a(3) = 6 binary matrices with 3 ones, with no zero rows or columns, up to row and column permutation:
  [1 0 0] [1 1 0] [1 0] [1 1] [1 1 1] [1]
  [0 1 0] [0 0 1] [1 0] [1 0] ....... [1].
  [0 0 1] ....... [0 1] ............. [1]
Non-isomorphic representatives of the a(3)=6 set multipartitions are: ((123)), ((1)(23)), ((2)(12)), ((1)(1)(1)), ((1)(2)(2)), ((1)(2)(3)). - _Gus Wiseman_, Mar 17 2017
		

Crossrefs

Programs

  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v,n,(-1)^(n-1)/n))))-1,-#v)}
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    K(q, t, k)={WeighT(Vec(sum(j=1, #q, gcd(t, q[j])*x^lcm(t, q[j])) + O(x*x^k), -k))}
    a(n)={my(s=0); forpart(q=n, s+=permcount(q)*polcoef(exp(x*Ser(sum(t=1, n, K(q, t, n)/t))), n)); s/n!} \\ Andrew Howroyd, Jan 16 2023

Formula

Calculate number of connected bipartite graphs + number of connected bipartite graphs with no duality automorphism, then apply EULER transform.
a(n) is the coefficient of x^n in the cycle index Z(S_n X S_n; 1+x, 1+x^2, ...), where S_n X S_n is Cartesian product of symmetric groups S_n of degree n.

Extensions

More terms and formula from Vladeta Jovovic, Jul 29 2000
a(19)-a(28) from Max Alekseyev, Jul 22 2009
a(29)-a(102) from Aliaksandr Siarhei, Dec 13 2013
Name edited by Gus Wiseman, Dec 18 2018

A120732 Number of square matrices with nonnegative integer entries and without zero rows or columns such that sum of all entries is equal to n.

Original entry on oeis.org

1, 1, 3, 15, 107, 991, 11267, 151721, 2360375, 41650861, 821881709, 17932031225, 428630422697, 11138928977049, 312680873171465, 9428701154866535, 303957777464447449, 10431949496859168189, 379755239311735494421
Offset: 0

Views

Author

Vladeta Jovovic, Aug 18 2006

Keywords

Examples

			From _Gus Wiseman_, Nov 14 2018: (Start)
The a(3) = 15 matrices:
  [3]
.
  [2 0] [1 1] [1 1] [1 0] [1 0] [0 2] [0 1] [0 1]
  [0 1] [1 0] [0 1] [1 1] [0 2] [1 0] [2 0] [1 1]
.
  [1 0 0] [1 0 0] [0 1 0] [0 1 0] [0 0 1] [0 0 1]
  [0 1 0] [0 0 1] [1 0 0] [0 0 1] [1 0 0] [0 1 0]
  [0 0 1] [0 1 0] [0 0 1] [1 0 0] [0 1 0] [1 0 0]
(End)
		

Crossrefs

Programs

  • Mathematica
    Table[1/n!*Sum[(-1)^(n-k)*StirlingS1[n,k]*Sum[(m!)^2*StirlingS2[k,m]^2,{m,0,k}],{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, May 07 2014 *)
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]]; Table[Length[Select[multsubs[Tuples[Range[n],2],n],Union[First/@#]==Union[Last/@#]==Range[Max@@First/@#]&]],{n,5}] (* Gus Wiseman, Nov 14 2018 *)

Formula

a(n) = (1/n!)*Sum_{k=0..n} (-1)^(n-k)*Stirling1(n,k)*A048144(k).
G.f.: Sum_{n>=0} Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*((1-x)^(-j)-1)^n.
a(n) ~ c * n! / (sqrt(n) * (log(2))^(2*n)), where c = 0.4670932578797312973586879293426... . - Vaclav Kotesovec, May 07 2014
In closed form, c = 2^(log(2)/2-2) / (log(2) * sqrt(Pi*(1-log(2)))). - Vaclav Kotesovec, May 03 2015
G.f.: Sum_{n>=0} (1-x)^n * (1 - (1-x)^n)^n. - Paul D. Hanna, Mar 26 2018

A104602 Number of square (0,1)-matrices with exactly n entries equal to 1 and no zero row or columns.

Original entry on oeis.org

1, 1, 2, 10, 70, 642, 7246, 97052, 1503700, 26448872, 520556146, 11333475922, 270422904986, 7016943483450, 196717253145470, 5925211960335162, 190825629733950454, 6543503207678564364, 238019066600097607402, 9153956822981328930170, 371126108428565106918404
Offset: 0

Views

Author

Ralf Stephan, Mar 27 2005

Keywords

Comments

Number of square (0,1)-matrices with exactly n entries equal to 1 and no zero row or columns, up to row and column permutation, is A057151(n). - Vladeta Jovovic, Mar 25 2006

Examples

			From _Gus Wiseman_, Nov 14 2018: (Start)
The a(3) = 10 matrices:
  [1 1] [1 1] [1 0] [0 1]
  [1 0] [0 1] [1 1] [1 1]
.
  [1 0 0] [1 0 0] [0 1 0] [0 1 0] [0 0 1] [0 0 1]
  [0 1 0] [0 0 1] [1 0 0] [0 0 1] [1 0 0] [0 1 0]
  [0 0 1] [0 1 0] [0 0 1] [1 0 0] [0 1 0] [1 0 0]
(End)
		

Crossrefs

Programs

  • Mathematica
    Table[1/n!*Sum[StirlingS1[n,k]*Sum[(m!)^2*StirlingS2[k, m]^2, {m, 0, k}],{k,0,n}],{n,1,20}] (* Vaclav Kotesovec, May 07 2014 *)
    Table[Length[Select[Subsets[Tuples[Range[n],2],{n}],Union[First/@#]==Union[Last/@#]==Range[Max@@First/@#]&]],{n,5}] (* Gus Wiseman, Nov 14 2018 *)

Formula

a(n) = (1/n!)*Sum_{k=0..n} Stirling1(n,k)*A048144(k). - Vladeta Jovovic, Mar 25 2006
G.f.: Sum_{n>=0} Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*((1+x)^j-1)^n. - Vladeta Jovovic, Mar 25 2006
a(n) ~ c * n! / (sqrt(n) * (log(2))^(2*n)), where c = 0.28889864564457451375789435201798... . - Vaclav Kotesovec, May 07 2014
In closed form, c = 1 / (log(2) * 2^(log(2)/2+2) * sqrt(Pi*(1-log(2)))). - Vaclav Kotesovec, May 03 2015
G.f.: Sum_{n>=0} ((1+x)^n - 1)^n / (1+x)^(n*(n+1)). - Paul D. Hanna, Mar 26 2018

Extensions

More terms from Vladeta Jovovic, Mar 25 2006
a(0)=1 prepended by Alois P. Heinz, Jan 14 2015

A054976 Number of binary n X n matrices with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

1, 3, 17, 179, 3835, 200082, 29610804, 13702979132, 20677458750966, 103609939177198046, 1745061194503344181714, 99860890306900024150675406, 19611238933283757244479826044874, 13340750149227624084760722122669739026, 31706433098827528779057124372265863803044450
Offset: 1

Views

Author

Vladeta Jovovic, May 27 2000

Keywords

Comments

Also the number of non-isomorphic set multipartitions (multisets of sets) with n parts and n vertices. - Gus Wiseman, Nov 18 2018

Examples

			From _Gus Wiseman_, Nov 18 2018: (Start)
Inequivalent representatives of the a(3) = 17 matrices:
  100 100 100 100 100 010 010 001 001 001 001 110 101 101 011 011 111
  100 010 001 011 011 001 101 001 101 011 111 101 011 011 011 111 111
  011 001 011 011 111 111 011 111 011 111 111 011 011 111 111 111 111
Non-isomorphic representatives of the a(1) = 1 through a(3) = 17 set multipartitions:
  {{1}}  {{1},{2}}      {{1},{2},{3}}
         {{2},{1,2}}    {{1},{1},{2,3}}
         {{1,2},{1,2}}  {{1},{3},{2,3}}
                        {{1},{2,3},{2,3}}
                        {{2},{1,3},{2,3}}
                        {{2},{3},{1,2,3}}
                        {{3},{1,3},{2,3}}
                        {{3},{3},{1,2,3}}
                        {{1,2},{1,3},{2,3}}
                        {{1},{2,3},{1,2,3}}
                        {{1,3},{2,3},{2,3}}
                        {{3},{2,3},{1,2,3}}
                        {{1,3},{2,3},{1,2,3}}
                        {{2,3},{2,3},{1,2,3}}
                        {{3},{1,2,3},{1,2,3}}
                        {{2,3},{1,2,3},{1,2,3}}
                        {{1,2,3},{1,2,3},{1,2,3}}
(End)
		

Crossrefs

Column sums of A057150.

Programs

Formula

a(n) = A002724(n) - 2*A002725(n-1) + A002724(n-1).

Extensions

More terms from David Wasserman, Mar 06 2002
Terms a(14) and beyond from Andrew Howroyd, Apr 11 2020

A057151 Number of square binary matrices with n ones, with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

1, 1, 2, 4, 8, 18, 41, 102, 252, 666, 1789, 5031, 14486, 43280, 132777, 420267, 1366307, 4566966, 15661086, 55081118, 198425478, 731661754, 2758808581, 10629386376, 41814350148, 167830018952, 686822393793, 2864024856054, 12162059027416, 52564545391789
Offset: 1

Views

Author

Vladeta Jovovic, Aug 14 2000

Keywords

Comments

Number of square binary matrices with n ones and with no zero rows or columns is A104602(n). - Vladeta Jovovic, Mar 25 2006
Also the number of non-isomorphic square set multipartitions (multisets of sets) of weight n. A multiset partition or hypergraph is square if its length (number of blocks or edges) is equal to its number of vertices. The weight of a multiset partition is the sum of sizes of its parts. - Gus Wiseman, Nov 16 2018

Examples

			There are 666 square binary matrices with 10 ones, with no zero rows or columns, up to row and column permutation: 33 of size 4 X 4, 248 of size 5 X 5, 288 of size 6 X 6, 79 of size 7 X 7, 15 of size 8 X 8, 2 of size 9 X 9 and 1 of size 10 X 10. Cf. A057150.
From _Gus Wiseman_, Nov 16 2018: (Start)
Non-isomorphic representatives of the a(1) = 1 through a(6) = 18 square set multipartitions:
  {1}  {1}{2}  {2}{12}    {12}{12}      {1}{23}{23}      {12}{13}{23}
               {1}{2}{3}  {1}{1}{23}    {2}{13}{23}      {1}{23}{123}
                          {1}{3}{23}    {2}{3}{123}      {13}{23}{23}
                          {1}{2}{3}{4}  {3}{13}{23}      {3}{23}{123}
                                        {3}{3}{123}      {1}{1}{1}{234}
                                        {1}{2}{2}{34}    {1}{1}{24}{34}
                                        {1}{2}{4}{34}    {1}{1}{4}{234}
                                        {1}{2}{3}{4}{5}  {1}{2}{34}{34}
                                                         {1}{3}{24}{34}
                                                         {1}{3}{4}{234}
                                                         {1}{4}{24}{34}
                                                         {1}{4}{4}{234}
                                                         {2}{4}{12}{34}
                                                         {3}{4}{12}{34}
                                                         {4}{4}{12}{34}
                                                         {1}{2}{3}{3}{45}
                                                         {1}{2}{3}{5}{45}
                                                         {1}{2}{3}{4}{5}{6}
(End)
		

Crossrefs

Extensions

More terms from Max Alekseyev, May 31 2007

A321721 Number of non-isomorphic non-normal semi-magic square multiset partitions of weight n.

Original entry on oeis.org

1, 1, 2, 2, 4, 2, 7, 2, 10, 7, 12, 2, 38, 2, 21, 46, 72, 2, 162, 2, 420, 415, 64, 2, 4987, 1858, 110, 9336, 45456, 2, 136018, 2, 1014658, 406578, 308, 3996977, 34937078, 2, 502, 28010167, 1530292965, 2, 508164038, 2, 54902992348, 51712929897, 1269, 2, 3217847072904, 8597641914, 9168720349613
Offset: 0

Views

Author

Gus Wiseman, Nov 18 2018

Keywords

Comments

A non-normal semi-magic square multiset partition of weight n is a multiset partition of weight n whose part sizes and vertex degrees are all equal to d, for some d|n.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
Also the number of nonnegative integer square matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, with row sums and column sums all equal to d, for some d|n.

Examples

			Non-isomorphic representatives of the a(2) = 2 through a(6) = 7 multiset partitions:
  {{11}}   {{111}}     {{1111}}       {{11111}}         {{111111}}
  {{1}{2}} {{1}{2}{3}} {{11}{22}}     {{1}{2}{3}{4}{5}} {{111}{222}}
                       {{12}{12}}                       {{112}{122}}
                       {{1}{2}{3}{4}}                   {{11}{22}{33}}
                                                        {{11}{23}{23}}
                                                        {{12}{13}{23}}
                                                        {{1}{2}{3}{4}{5}{6}}
Inequivalent representatives of the a(6) = 7 matrices:
  [6]
.
  [3 0] [2 1]
  [0 3] [1 2]
.
  [2 0 0] [2 0 0] [1 1 0]
  [0 2 0] [0 1 1] [1 0 1]
  [0 0 2] [0 1 1] [0 1 1]
.
  [1 0 0 0 0 0]
  [0 1 0 0 0 0]
  [0 0 1 0 0 0]
  [0 0 0 1 0 0]
  [0 0 0 0 1 0]
  [0 0 0 0 0 1]
Inequivalent representatives of the a(9) = 7 matrices:
  [9]
.
  [3 0 0] [3 0 0] [2 1 0] [2 1 0] [1 1 1]
  [0 3 0] [0 2 1] [1 1 1] [1 0 2] [1 1 1]
  [0 0 3] [0 1 2] [0 1 2] [0 2 1] [1 1 1]
.
  [1 0 0 0 0 0 0 0 0]
  [0 1 0 0 0 0 0 0 0]
  [0 0 1 0 0 0 0 0 0]
  [0 0 0 1 0 0 0 0 0]
  [0 0 0 0 1 0 0 0 0]
  [0 0 0 0 0 1 0 0 0]
  [0 0 0 0 0 0 1 0 0]
  [0 0 0 0 0 0 0 1 0]
  [0 0 0 0 0 0 0 0 1]
		

Crossrefs

Formula

a(p) = 2 for p prime corresponding to the 1 X 1 square [p] and the permutation matrices of size p X p with partition (1...10...0). - Chai Wah Wu, Jan 16 2019
a(n) = Sum_{d|n} A333733(d,n/d) for n > 0. - Andrew Howroyd, Apr 11 2020

Extensions

a(11)-a(13) from Chai Wah Wu, Jan 16 2019
a(14)-a(15) from Chai Wah Wu, Jan 20 2019
Terms a(16) and beyond from Andrew Howroyd, Apr 11 2020

A104601 Triangle T(r,n) read by rows: number of n X n (0,1)-matrices with exactly r entries equal to 1 and no zero row or columns.

Original entry on oeis.org

1, 0, 2, 0, 4, 6, 0, 1, 45, 24, 0, 0, 90, 432, 120, 0, 0, 78, 2248, 4200, 720, 0, 0, 36, 5776, 43000, 43200, 5040, 0, 0, 9, 9066, 222925, 755100, 476280, 40320, 0, 0, 1, 9696, 727375, 6700500, 13003620, 5644800, 362880, 0, 0, 0, 7480, 1674840
Offset: 1

Views

Author

Ralf Stephan, Mar 27 2005

Keywords

Examples

			1
0,2
0,4,6
0,1,45,24
0,0,90,432,120
0,0,78,2248,4200,720
0,0,36,5776,43000,43200,5040
0,0,9,9066,222925,755100,476280,40320
0,0,1,9696,727375,6700500,13003620,5644800,362880
0,0,0,7480,1674840,37638036,179494350,226262400,71850240,3628800
		

Crossrefs

Right-edge diagonals include A000142, A055602, A055603. Row sums are in A104602.
Column sums are in A048291. The triangle read by columns = A055599.

Programs

  • Mathematica
    t[r_, n_] := Sum[ Sum[ (-1)^(2n - d - k/d)*Binomial[n, d]*Binomial[n, k/d]*Binomial[k, r], {d, Divisors[k]}], {k, r, n^2}]; Flatten[ Table[t[r, n], {r, 1, 10}, {n, 1, r}]] (* Jean-François Alcover, Jun 27 2012, from formula *)
    Table[Length[Select[Subsets[Tuples[Range[n],2],{n}],Union[First/@#]==Union[Last/@#]==Range[k]&]],{n,6},{k,n}] (* Gus Wiseman, Nov 14 2018 *)

Formula

T(r, n) = Sum{l>=r, Sum{d|l, (-1)^(2n-d-l/d)*C(n, d)*C(n, l/d)*C(l, r) }}.
E.g.f.: Sum(((1+x)^n-1)^n*exp((1-(1+x)^n)*y)*y^n/n!,n=0..infinity). - Vladeta Jovovic, Feb 24 2008

A057149 Triangle T(n,k) of n X n binary matrices with k ones, with no zero rows or columns, up to row and column permutation.

Original entry on oeis.org

0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 2, 5, 4, 3, 1, 1, 0, 0, 0, 0, 1, 2, 11, 21, 34, 33, 33, 19, 14, 6, 3, 1, 1, 0, 0, 0, 0, 0, 1, 2, 14, 49, 131, 248, 410, 531, 601, 566, 474, 336, 222, 124, 67, 32, 16, 6, 3, 1, 1, 0, 0, 0, 0, 0, 0, 1, 2, 15, 69, 288, 840, 2144, 4488, 8317, 13160
Offset: 1

Views

Author

Vladeta Jovovic, Aug 14 2000

Keywords

Comments

Row sums give A054976.

Examples

			[0,1], [0,0,1,1,1], [0,0,0,1,2,5,4,3,1,1],...;
T(4,6)=11, i.e. there are 11 4 X 4 binary matrices with 6 ones, with no zero rows or columns, up to row and column permutation:
[0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1]
[0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 1 0]
[0 0 0 1] [0 0 1 0] [0 0 1 0] [0 0 1 1] [0 1 1 0] [0 0 1 1]
[1 1 1 0] [1 1 0 1] [1 1 1 0] [1 1 0 0] [1 0 1 0] [1 1 0 0]
and
[0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1] [0 0 0 1]
[0 0 1 0] [0 0 1 0] [0 0 1 0] [0 0 1 0] [0 0 1 0]
[0 1 0 0] [0 1 0 1] [0 1 0 1] [0 1 0 1] [1 1 0 0]
[1 0 1 1] [1 0 0 1] [1 0 1 0] [1 1 0 0] [1 1 0 0].
		

Crossrefs

A321724 Irregular triangle read by rows where T(n,k) is the number of non-isomorphic non-normal semi-magic square multiset partitions of weight n and length d = A027750(n, k).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 3, 5, 1, 1, 5, 1, 1, 3, 7, 1, 1, 1, 1, 4, 9, 12, 11, 1, 1, 1, 1, 4, 15, 1, 1, 13, 31, 1, 1, 5, 43, 22, 1, 1, 1, 1, 5, 22, 103, 30, 1, 1, 1, 1, 6, 106, 264, 42, 1, 1, 30, 383, 1, 1, 6, 56, 1, 1, 1, 1, 7, 45, 321, 2804, 1731, 77, 1
Offset: 1

Views

Author

Gus Wiseman, Nov 18 2018

Keywords

Comments

Also the number of nonnegative integer square matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, with row sums and column sums all equal to d.
A non-normal semi-magic square multiset partition of weight n is a multiset partition of weight n whose part sizes and vertex degrees are all equal to d, for some d|n.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Triangle begins:
  1
  1 1
  1 1
  1 2 1
  1 1
  1 2 3 1
  1 1
  1 3 5 1
  1 5 1
  1 3 7 1
Inequivalent representatives of the T(10,3) = 7 semi-magic squares (zeros not shown):
  [2    ] [2    ] [2    ] [2    ] [2    ] [11   ] [11   ]
  [ 2   ] [ 2   ] [ 2   ] [ 11  ] [ 11  ] [11   ] [1 1  ]
  [  2  ] [  2  ] [  11 ] [ 11  ] [ 1 1 ] [  11 ] [ 1 1 ]
  [   2 ] [   11] [  1 1] [   11] [  1 1] [  1 1] [  1 1]
  [    2] [   11] [   11] [   11] [   11] [   11] [   11]
		

Crossrefs

Formula

T(n,k) = A333733(d, n/d), where d = A027750(n, k). - Andrew Howroyd, Apr 11 2020

Extensions

a(28)-a(39) from Chai Wah Wu, Jan 16 2019
Terms a(40) and beyond from Andrew Howroyd, Apr 11 2020
Edited by Peter Munn, Mar 05 2025

A321609 Array read by antidiagonals: T(n,k) is the number of inequivalent binary n X n matrices with k ones, under row and column permutations.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 3, 1, 1, 0, 0, 1, 3, 1, 1, 0, 0, 1, 6, 3, 1, 1, 0, 0, 0, 7, 6, 3, 1, 1, 0, 0, 0, 7, 16, 6, 3, 1, 1, 0, 0, 0, 6, 21, 16, 6, 3, 1, 1, 0, 0, 0, 3, 39, 34, 16, 6, 3, 1, 1, 0, 0, 0, 1, 44, 69, 34, 16, 6, 3, 1, 1, 0, 0, 0, 1, 55, 130, 90, 34, 16, 6, 3, 1, 1
Offset: 0

Views

Author

Andrew Howroyd, Nov 14 2018

Keywords

Examples

			Array begins:
==========================================================
n\k| 0  1  2  3  4  5  6   7   8    9   10    11    12
---+------------------------------------------------------
0  | 1  0  0  0  0  0  0   0   0    0    0     0     0 ...
1  | 1  1  0  0  0  0  0   0   0    0    0     0     0 ...
2  | 1  1  3  1  1  0  0   0   0    0    0     0     0 ...
3  | 1  1  3  6  7  7  6   3   1    1    0     0     0 ...
4  | 1  1  3  6 16 21 39  44  55   44   39    21    16 ...
5  | 1  1  3  6 16 34 69 130 234  367  527   669   755 ...
6  | 1  1  3  6 16 34 90 182 425  870 1799  3323  5973 ...
7  | 1  1  3  6 16 34 90 211 515 1229 2960  6893 15753 ...
8  | 1  1  3  6 16 34 90 211 558 1371 3601  9209 24110 ...
9  | 1  1  3  6 16 34 90 211 558 1430 3825 10278 28427 ...
...
		

Crossrefs

Rows n=6..8 are A052370, A053304, A053305.
Main diagonal is A049311.
Row sums are A002724.
Cf. A052371 (as triangle), A057150, A246106, A318795.

Programs

  • Mathematica
    permcount[v_List] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
    c[p_List, q_List, k_] := SeriesCoefficient[Product[Product[(1 + O[x]^(k + 1) + x^LCM[p[[i]], q[[j]]])^GCD[p[[i]], q[[j]]], {j, 1, Length[q]}], {i, 1, Length[p]}], {x, 0, k}];
    M[m_, n_, k_] := Module[{s = 0}, Do[Do[s += permcount[p]*permcount[q]*c[p, q, k], {q, IntegerPartitions[n]}], {p, IntegerPartitions[m]}]; s/(m!*n!)]
    Table[M[n - k, n - k, k], {n, 0, 12}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Sep 10 2019, after Andrew Howroyd *)
  • PARI
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    c(p, q, k)={polcoef(prod(i=1, #p, prod(j=1, #q, (1 + x^lcm(p[i], q[j]) + O(x*x^k))^gcd(p[i], q[j]))), k)}
    M(m, n, k)={my(s=0); forpart(p=m, forpart(q=n, s+=permcount(p) * permcount(q) * c(p, q, k))); s/(m!*n!)}
    for(n=0, 10, for(k=0, 12, print1(M(n, n, k), ", ")); print); \\ Andrew Howroyd, Nov 14 2018

Formula

T(n,k) = T(k,k) for n > k.
T(n,k) = 0 for k > n^2.
Showing 1-10 of 12 results. Next