A120733 -id:A120733 - cp's OEIS frontend

A007716 Number of polynomial symmetric functions of matrix of order n under separate row and column permutations.

1, 1, 4, 10, 33, 91, 298, 910, 3017, 9945, 34207, 119369, 429250, 1574224, 5916148, 22699830, 89003059, 356058540, 1453080087, 6044132794, 25612598436, 110503627621, 485161348047, 2166488899642, 9835209912767, 45370059225318, 212582817739535, 1011306624512711

Offset: 0

Colin Mallows

Also, the number of nonnegative integer n X n matrices with sum of elements equal to n, under row and column permutations (cf. A120733).
This is a two-dimensional generalization of the partition function (A000041), which equals the number of length n vectors of nonnegative integers with sum n, equivalent under permutations. - Franklin T. Adams-Watters, Sep 19 2011
Also number of non-isomorphic multiset partitions of weight n. - Gus Wiseman, Sep 19 2011

			The 10 non-isomorphic multiset partitions of weight 3 are {{1, 1, 1}}, {{1, 1, 2}}, {{1, 2, 3}}, {{1}, {1, 1}}, {{1}, {1, 2}}, {{1}, {2, 2}}, {{1}, {2, 3}}, {{1}, {1}, {1}}, {{1}, {1}, {2}}, {{1}, {2}, {3}}.

Andrew Howroyd, Table of n, a(n) for n = 0..50 (terms 0..30 from Seiichi Manyama)

Main diagonal of A318795.

Cf. A053307, A052365, A052366, A052367, A052372, A052373, A049311, A054688, A000041.

permcount[v_] := 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_, q_, k_] := SeriesCoefficient[1/Product[(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!)];
a[n_] := a[n] = M[n, n, n];
Table[Print[n, " ", a[n]]; a[n], {n, 0, 18}] (* Jean-François Alcover, May 03 2019, after Andrew Howroyd *)

\\ See A318795
a(n) = M(n,n,n); \\ Andrew Howroyd, Sep 03 2018

EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,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)={EulerT(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, Mar 29 2020

a(n) is the coefficient of x^n in the cycle index Z(S_n X S_n; x_1, x_2, ...) if we replace x_i with 1+x^i+x^(2*i)+x^(3*i)+x^(4*i)+..., where S_n X S_n is the Cartesian product of symmetric groups S_n of degree n. - Vladeta Jovovic, Mar 09 2000

More terms from Vladeta Jovovic, Jun 28 2000
a(19)-a(25) from Max Alekseyev, Jan 22 2010
a(0)=1 prepended by Alois P. Heinz, Feb 03 2019
a(26)-a(27) from Seiichi Manyama, Nov 23 2019

A138178 Number of symmetric 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, 9, 33, 125, 531, 2349, 11205, 55589, 291423, 1583485, 8985813, 52661609, 319898103, 2000390153, 12898434825, 85374842121, 580479540219, 4041838056561, 28824970996809, 210092964771637, 1564766851282299, 11890096357039749, 92151199272181629

Offset: 0

Vladeta Jovovic, Mar 03 2008

Number of normal semistandard Young tableaux of size n, where a tableau is normal if its entries span an initial interval of positive integers. - Gus Wiseman, Feb 23 2018

			a(4) = 33 because there are 1 such matrix of type 1 X 1, 7 matrices of type 2 X 2, 15 of type 3 X 3 and 10 of type 4 X 4, cf. A138177.
From _Gus Wiseman_, Feb 23 2018: (Start)
The a(3) = 9 normal semistandard Young tableaux:
1   1 2   1 3   1 2   1 1   1 2 3   1 2 2   1 1 2   1 1 1
2   3     2     2     2
3
(End)
From _Gus Wiseman_, Nov 14 2018: (Start)
The a(4) = 33 matrices:
[4]
.
[30][21][20][11][10][02][01]
[01][10][02][11][03][20][12]
.
[200][200][110][101][100][100][100][100][011][010][010][010][001][001][001]
[010][001][100][010][020][011][010][001][100][110][101][100][020][010][001]
[001][010][001][100][001][010][002][011][100][001][010][002][100][101][110]
.
[1000][1000][1000][1000][0100][0100][0010][0010][0001][0001]
[0100][0100][0010][0001][1000][1000][0100][0001][0100][0010]
[0010][0001][0100][0010][0010][0001][1000][1000][0010][0100]
[0001][0010][0001][0100][0001][0010][0001][0100][1000][1000]
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..500

Row sums of A138177.
Cf. A007716, A120733, A135588, A296188.
Cf. A057151, A104601, A104602, A120732, A316983, A320796, A321401, A321405, A321407.

gf:= proc(j) local k, n; add(add((-1)^(n-k) *binomial(n, k) *(1-x)^(-k) *(1-x^2)^(-binomial(k, 2)), k=0..n), n=0..j) end: a:= n-> coeftayl(gf(n+1), x=0, n): seq(a(n), n=0..25); # Alois P. Heinz, Sep 25 2008

Table[Sum[SeriesCoefficient[1/(2^(k+1)*(1-x)^k*(1-x^2)^(k*(k-1)/2)),{x,0,n}],{k,0,Infinity}],{n,0,20}]  (* Vaclav Kotesovec, Jul 03 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],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],Sort[Reverse/@#]==#]&]],{n,5}] (* Gus Wiseman, Nov 14 2018 *)

G.f.: Sum_{n>=0} Sum_{k=0..n} (-1)^(n-k)*C(n,k)*(1-x)^(-k)*(1-x^2)^(-C(k,2)).

G.f.: Sum_{n>=0} 2^(-n-1)*(1-x)^(-n)*(1-x^2)^(-C(n,2)). - Vladeta Jovovic, Dec 09 2009

More terms from Alois P. Heinz, Sep 25 2008

A101370 Number of zero-one matrices with n ones and no zero rows or columns.

Original entry on oeis.org

1, 4, 24, 196, 2016, 24976, 361792, 5997872, 111969552, 2324081728, 53089540992, 1323476327488, 35752797376128, 1040367629940352, 32441861122796672, 1079239231677587264, 38151510015777089280, 1428149538870997774080, 56435732691153773665280

Offset: 1

Peter J. Cameron, Jan 14 2005

a(n) = (1/(4*n!)) * Sum_{r, s>=0} (r*s)_n / 2^(r+s), where (m)_n is the falling factorial m * (m-1) * ... * (m-n+1). [Maia and Mendez]

			a(2)=4:
[1 1] [1] [1 0] [0 1]
..... [1] [0 1] [1 0]
From _Gus Wiseman_, Nov 14 2018: (Start)
The a(3) = 24 matrices:
  [111]
.
  [11][11][110][101][10][100][011][01][010][001]
  [10][01][001][010][11][011][100][11][101][110]
.
  [1][10][10][10][100][100][01][01][010][01][010][001][001]
  [1][10][01][01][010][001][10][10][100][01][001][100][010]
  [1][01][10][01][001][010][10][01][001][10][100][010][100]
(End)

Georg Cantor, Gesammelte Abhandlungen mathematischen und philosophischen Inhalts, p. 435 (IV, 4. Mitteilungen zur Lehre vom Transfiniten, VIII Nr. 13), Springer, Berlin. [Rainer Rosenthal, Apr 10 2007]

Alois P. Heinz, Table of n, a(n) for n = 1..400
P. J. Cameron, D. A. Gewurz and F. Merola, Product action, Discrete Math., 308 (2008), 386-394.
Giulio Cerbai and Anders Claesson, Enumerative aspects of Caylerian polynomials, arXiv:2411.08426 [math.CO], 2024. See pp. 3, 19.
Loïc Foissy, Claudia Malvenuto, and Frédéric Patras, Matrix symmetric and quasi-symmetric functions and noncommutative representation theory, arXiv:2503.14417 [math.CO], 2025. See p. 20.
M. Maia and M. Mendez, On the arithmetic product of combinatorial species, arXiv:math/0503436 [math.CO], 2005.

Cf. A000670 (the sequence P(n)), A049311 (row and column permutations allowed), A120733, A122725, A135589, A283877, A321446, A321587.

P:=function(n) return Sum([1..n],x->Stirling2(n,x)*Factorial(x)); end;

F:=function(n) return Sum([1..n],x->(-1)^(n-x)*Stirling1(n,x)*P(x)^2)/Factorial(n); end;

m = 17; a670[n_] = Sum[ StirlingS2[n, k]*k!, {k, 0, n}]; Rest[ CoefficientList[ Series[ Sum[ a670[n]^2*(Log[1 + x]^n/n!), {n, 0, m}], {x, 0, m}], x]] (* Jean-François Alcover, Sep 02 2011, after g.f.  *)
Table[Length[Select[Subsets[Tuples[Range[n],2],{n}],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#]]&]],{n,5}] (* Gus Wiseman, Nov 14 2018 *)

{A000670(n)=sum(k=0,n,stirling(n, k,2)*k!)}
{a(n)=polcoeff(sum(m=0,n,A000670(m)^2*log(1+x+x*O(x^n))^m/m!),n)}
/* Paul D. Hanna, Nov 07 2009 */

a(n) = (Sum s(n, k) * P(k)^2)/n!, where P(n) is the number of labeled total preorders on {1, ..., n} (A000670), s are signed Stirling numbers of the first kind.
G.f.: Sum_{m>=0,n>=0} Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*((1+x)^j-1)^m. - Vladeta Jovovic, Mar 25 2006
Inverse binomial transform of A007322. - Vladeta Jovovic, Aug 17 2006
G.f.: Sum_{n>=0} 1/(2-(1+x)^n)/2^(n+1). - Vladeta Jovovic, Sep 23 2006
G.f.: Sum_{n>=0} A000670(n)^2*log(1+x)^n/n! where 1/(1-x) = Sum_{n>=0} A000670(n)*log(1+x)^n/n!. - Paul D. Hanna, Nov 07 2009
a(n) ~ n! / (2^(2+log(2)/2) * (log(2))^(2*(n+1))). - Vaclav Kotesovec, Dec 31 2013

A114736 Number of planar partitions of n where parts strictly decrease along each row and column.

Original entry on oeis.org

1, 1, 1, 3, 4, 6, 10, 15, 22, 33, 49, 70, 102, 146, 205, 290, 405, 561, 779, 1071, 1463, 1999, 2714, 3667, 4946, 6641, 8880, 11848, 15753, 20870, 27586, 36354, 47766, 62621, 81878, 106785, 138975, 180449, 233778, 302270, 390027, 502256, 645603, 828330, 1060851

Offset: 0

Franklin T. Adams-Watters, Mar 16 2006

nonn

If these partitions are "flattened" into a simple partition, the resulting partitions are those for which any part size present with multiplicity k implies the presence of at least k(k-1)/2 larger parts. E.g., [3,1|1] flattens to [3,1^2], 1 has multiplicity 2, so there must be at least 2*1/2 = 1 part larger than 1 - which is the 3.

			For n = 5, we have the 6 partitions [5], [4,1], [4|1], [3,2], [3|2] and [3,1|1].
From _Gus Wiseman_, Nov 15 2018: (Start)
The a(6) = 10 plane partitions:
  6   5 1   4 2   3 2 1
.
  5   4 1   4   3 2   3 1
  1   1     2   1     2
.
  3
  2
  1
(End)

B. Gordon, Multirowed partitions with strict decrease along columns (Notes on plane partitions IV.), Symposia Amer. Math. Soc. 19 (1971) 91-100.

Alois P. Heinz, Table of n, a(n) for n = 0..85

Cf. A000009, A000219, A001970, A007716, A068313, A117433, A120733, A319646, A321645, A321652, A321653, A321655.

prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
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],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],And@@(OrderedQ[#,Greater]&/@prs2mat[#]),And@@(OrderedQ[#,Greater]&/@Transpose[prs2mat[#]])]&]],{n,5}] (* Gus Wiseman, Nov 15 2018 *)

Clarified definition, added 30 terms and reference. - Dennis K Moore, Jan 12 2011

a(40)-a(44) from Alois P. Heinz, Sep 26 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

Vladeta Jovovic, Aug 18 2006

nonn

			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)

G. C. Greubel, Table of n, a(n) for n = 0..250

Cf. A007716, A048291, A054976, A057149, A057150, A057151, A104601, A104602, A120733, A138178, A316983, A319616.

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 *)

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

A117433 Number of planar partitions of n with all part sizes distinct.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 11, 15, 21, 35, 41, 59, 75, 103, 149, 187, 243, 321, 413, 527, 735, 895, 1165, 1467, 1885, 2335, 2997, 3853, 4765, 5977, 7473, 9269, 11531, 14255, 17537, 22201, 26897, 33233, 40613, 50027, 60637, 74459, 89963, 109751, 134407, 162117, 195859

Offset: 0

Franklin T. Adams-Watters, Mar 16 2006, Apr 01 2008

nonn

Matches A072706 for n < 10, since a unimodal composition into distinct parts can be placed uniquely as a hook. Starting with n = 10, additional partitions are possible (starting with [4,3|2,1] and [4,2|3,1]).

			From _Gus Wiseman_, Nov 15 2018: (Start)
The a(10) = 35 strict plane partitions (A = 10):
  A  64  73  82  532  91  541  631  721  4321
.
  9  54  63  72  432  8  53  71  431  7  43  52  61  421  6  42  51
  1  1   1   1   1    2  2   2   2    3  21  3   3   3    4  31  4
.
  7  6  5  43  42  5  41
  2  3  4  2   3   3  3
  1  1  1  1   1   2  2
.
  4
  3
  2
  1
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..1000 (first 100 terms from Franklin T. Adams-Watters)
OEIS Wiki, Plane partitions

Cf. A000219, A072706, A117434, A000009.

Cf. A001970, A007716, A068313, A114736, A120733, A319646, A321645, A321652, A321653, A321655, A321659, A321660.

b:= proc(n, i) b(n, i):= `if`(n=0, [1], `if`(i<1, [], zip((x, y)
      -> x+y, b(n, i-1), `if`(i>n, [], [0, b(n-i, i-1)[]]), 0)))
    end:
g:= proc(n) g(n):= `if`(n<2, 1, (n-1)*g(n-2) +g(n-1)) end:
a:= proc(n) b(n, n); add(%[i]*g(i-1), i=1..nops(%)) end:
seq(a(n), n=0..60);  # Alois P. Heinz, Nov 18 2012

prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
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],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],UnsameQ@@DeleteCases[Join@@prs2mat[#],0],And@@(OrderedQ[#,Greater]&/@prs2mat[#]),And@@(OrderedQ[#,Greater]&/@Transpose[prs2mat[#]])]&]],{n,5}] (* Gus Wiseman, Nov 15 2018 *)
zip[f_, x_List, y_List, z_] := With[{m = Max[Length[x], Length[y]]}, f[PadRight[x, m, z], PadRight[y, m, z]]];
b[n_, i_] := b[n, i] = If[n == 0, {1}, If[i < 1, {}, zip[Plus, b[n, i - 1], If[i > n, {}, Join[{0}, b[n - i, i - 1]]], 0]]];
g[n_] := g[n] = If[n < 2, 1, (n - 1)*g[n - 2] + g[n - 1]];
a[n_] := With[{bn = b[n, n]}, Sum[bn[[i]]*g[i - 1], {i, 1, Length[bn]}]];
Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Dec 05 2023, after Alois P. Heinz *)

a(n) = Sum_{k=1..floor((sqrt(8*n+1)-1)/2)} A000085(k)*A008289(n,k).

A321719 Number of non-normal semi-magic squares with sum of entries equal to n.

Original entry on oeis.org

1, 1, 3, 7, 28, 121, 746, 5041, 40608, 362936, 3635017, 39916801, 479206146, 6227020801, 87187426839, 1307674521272, 20923334906117, 355687428096001, 6402415241245577, 121645100408832001, 2432905938909013343, 51090942176372298027, 1124001180562929946213

Offset: 0

Gus Wiseman, Nov 18 2018

nonn

A non-normal semi-magic square is a nonnegative integer matrix with row sums and column sums all equal to d, for some d|n.

Squares must be of size k X k where k is a divisor of n. This implies that a(p) = p! + 1 for p prime since the only allowable squares are of sizes 1 X 1 and p X p. The 1 X 1 square is [p], the p X p squares are necessarily permutation matrices and there are p! permutation matrices of size p X p. Also, a(n) >= n! + 1 for n > 1. - Chai Wah Wu, Jan 13 2019

			The a(3) = 7 semi-magic squares:
  [3]
.
  [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]

Andrew Howroyd, Table of n, a(n) for n = 0..100
Wikipedia, Magic square
Index entries for sequences related to magic squares

Cf. A006052, A120732, A120733, A257493, A271103, A319056, A319616.

Cf. A321717, A321718, A321720, A321721, A321722, A321723, A321724.

prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
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],And[Union[First/@#]==Range[Max@@First/@#]==Union[Last/@#],SameQ@@Total/@prs2mat[#],SameQ@@Total/@Transpose[prs2mat[#]]]&]],{n,5}]

a(p) = p! + 1 for p prime and a(n) >= n! + 1 for n > 1 (see comment above). - Chai Wah Wu, Jan 13 2019

a(n) = Sum_{d|n} A257493(d, n/d) for n > 0. - Andrew Howroyd, Apr 11 2020

a(7) from Chai Wah Wu, Jan 13 2019
a(6) corrected and a(8)-a(15) added by Chai Wah Wu, Jan 14 2019
a(16)-a(19) from Chai Wah Wu, Jan 16 2019
Terms a(20) and beyond from Andrew Howroyd, Apr 11 2020

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

Ralf Stephan, Mar 27 2005

nonn

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

			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)

Alois P. Heinz, Table of n, a(n) for n = 0..400
H. Cheballah, S. Giraudo, R. Maurice, Combinatorial Hopf algebra structure on packed square matrices, arXiv preprint arXiv:1306.6605 [math.CO], 2013-2015.
M. Maia and M. Mendez, On the arithmetic product of combinatorial species, arXiv:math/0503436 [math.CO], 2005.

Row sums of triangle A104601.

Cf. A048291, A049311, A054976, A057150, A057151, A101370, A120732, A120733, A138178, A316983, A319616.

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 *)

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

More terms from Vladeta Jovovic, Mar 25 2006

a(0)=1 prepended by Alois P. Heinz, Jan 14 2015

A261781 Number T(n,k) of compositions of n where each part i is marked with a word of length i over a k-ary alphabet whose letters appear in alphabetical order and all k letters occur at least once in the composition; triangle T(n,k), n >= 0, 0 <= k <= n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 2, 3, 0, 4, 16, 13, 0, 8, 66, 132, 75, 0, 16, 248, 924, 1232, 541, 0, 32, 892, 5546, 13064, 13060, 4683, 0, 64, 3136, 30720, 114032, 195020, 155928, 47293, 0, 128, 10888, 162396, 893490, 2327960, 3116220, 2075948, 545835

Offset: 0

Alois P. Heinz, Aug 31 2015

From Vaclav Kotesovec, Oct 14 2017: (Start)
Conjecture: For k > 0 the recurrence order for column k is equal to k*(k+1)/2.
Column k > 0 is asymptotic to c(k) * d(k)^n, where c(k) and d(k) are constants (dependent only on k).
k                           c(k)                            d(k)
1  A131577(n) ~ 0.50000000000000000000000000 * 2.00000000000000000000000000^n.
2  A293579(n) ~ 0.60355339059327376220042218 * 3.41421356237309504880168872^n.
3  A293580(n) ~ 0.64122035031051210658648604 * 4.84732210186307263951891624^n.
4  A293581(n) ~ 0.66065168848540565019767995 * 6.28521350788324520158143964^n.
5  A293582(n) ~ 0.67250239588725756267924287 * 7.72502395887257562679242875^n.
6  A293583(n) ~ 0.68048292906885160660288253 * 9.16579514882621927923459043^n.
7  A293584(n) ~ 0.68622254929933439577377124 * 10.6071156901906815408327973^n.
8  A293585(n) ~ 0.69054873168854973836384871 * 12.0487797070167958138215794^n.
9  A293586(n) ~ 0.69392626461456654033893782 * 13.4906727630621977261008808^n.
10 A293587(n) ~ 0.69663630864564830007443110 * 14.9327261729129660014886221^n.
---
Conjecture: d(k+1) - d(k) tends to 1/log(2).
d(2) - d(1) = 1.414213562373095048801688724209698...
d(3) - d(2) = 1.433108539489977590717227522340838...
d(4) - d(3) = 1.437891406020172562062523400686067...
d(5) - d(4) = 1.439810450989330425210989107036901...
d(6) - d(5) = 1.440771189953643652442161677346934...
d(7) - d(6) = 1.441320541364462261598206961226199...
d(8) - d(7) = 1.441664016826114272988782079622148...
d(9) - d(8) = 1.441893056045401912279301345910755...
d(10)- d(9) = 1.442053409850768275387741352145193...
1 / log(2)  = 1.442695040888963407359924681001892...
(End)

			A(3,2) = 16: 3aab, 3abb, 2aa1b, 2ab1a, 2ab1b, 2bb1a, 1a2ab, 1a2bb, 1b2aa, 1b2ab, 1a1a1b, 1a1b1a, 1a1b1b, 1b1a1a, 1b1a1b, 1b1b1a.
Triangle T(n,k) begins:
  1;
  0,  1;
  0,  2,    3;
  0,  4,   16,    13;
  0,  8,   66,   132,     75;
  0, 16,  248,   924,   1232,    541;
  0, 32,  892,  5546,  13064,  13060,   4683;
  0, 64, 3136, 30720, 114032, 195020, 155928, 47293;
  ...

Alois P. Heinz, Rows n = 0..140, flattened
E. Munarini, M. Poneti, and S. Rinaldi, Matrix compositions, JIS 12 (2009) 09.4.8, Table 2.

Columns k=0..10 give A000007, A131577, A293579, A293580, A293581, A293582, A293583, A293584, A293585, A293586, A293587.
Row sums give A120733.
Main diagonal gives A000670.
T(2n,n) gives A261784.
T(n+1,n)/2 gives A083385.
Cf. A261719 (same for partitions), A261780.

A:= proc(n, k) option remember; `if`(n=0, 1,
      add(A(n-j, k)*binomial(j+k-1, k-1), j=1..n))
    end:
T:= (n, k)-> add(A(n, k-i)*(-1)^i*binomial(k, i), i=0..k):
seq(seq(T(n, k), k=0..n), n=0..10);

A[n_, k_] := A[n, k] = If[n==0, 1,
    Sum[A[n-j, k]*Binomial[j+k-1, k-1], {j, 1, n}]];
T[n_, k_] := Sum[A[n, k-i]*(-1)^i*Binomial[k, i], {i, 0, k}];
Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 08 2017, translated from Maple *)

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

A323300 Number of ways to fill a matrix with the parts of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 2, 2, 4, 1, 6, 1, 4, 4, 3, 1, 6, 1, 6, 4, 4, 1, 12, 2, 4, 2, 6, 1, 12, 1, 2, 4, 4, 4, 18, 1, 4, 4, 12, 1, 12, 1, 6, 6, 4, 1, 10, 2, 6, 4, 6, 1, 12, 4, 12, 4, 4, 1, 36, 1, 4, 6, 4, 4, 12, 1, 6, 4, 12, 1, 20, 1, 4, 6, 6, 4, 12, 1, 10, 3, 4

Offset: 1

Gus Wiseman, Jan 12 2019

nonn

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

			The a(24) = 12 matrices whose entries are (2,1,1,1):
  [1 1 1 2] [1 1 2 1] [1 2 1 1] [2 1 1 1]
.
  [1 1] [1 1] [1 2] [2 1]
  [1 2] [2 1] [1 1] [1 1]
.
  [1] [1] [1] [2]
  [1] [1] [2] [1]
  [1] [2] [1] [1]
  [2] [1] [1] [1]

Positions of 1's are one and prime numbers A008578.
Positions of 2's are primes to prime powers A053810.
Cf. A000005, A001222, A008480, A056239, A063989, A112798, A120733.
Cf. A323295, A323305, A323307, A323351.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
Array[Length[ptnmats[#]]&,100]

a(n) = A008480(n) * A000005(A001222(n)).

cp's OEIS Frontend

A007716 Number of polynomial symmetric functions of matrix of order n under separate row and column permutations.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A101370 Number of zero-one matrices with n ones and no zero rows or columns.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

GAP

Mathematica

PARI

Formula

A114736 Number of planar partitions of n where parts strictly decrease along each row and column.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Extensions

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.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A117433 Number of planar partitions of n with all part sizes distinct.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A321719 Number of non-normal semi-magic squares with sum of entries equal to n.

Views

Author