A321732 Number of nonnegative integer square matrices with sum of entries equal to n, no zero rows or columns, and the same row sums as column sums.

1, 1, 3, 11, 53, 317, 2293, 19435, 188851, 2068417, 25203807, 338117445, 4951449055, 78589443061, 1343810727205, 24626270763109, 481489261372381, 10004230113283129, 220125503239710879, 5113204953106107087, 125037079246130168973

Offset

0,3

Links

Ludovic Schwob, Table of n, a(n) for n = 0..30

Example

The a(3) = 11 matrices:
  [3]
.
  [2 0] [1 1] [1 0] [0 1]
  [0 1] [1 0] [0 2] [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]

Mathematica

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/@#], Total/@prs2mat[#]==Total/@Transpose[prs2mat[#]]]&]], {n, 5}]

Crossrefs

Cf. A000700, A000701, A006052, A007016, A120732, A319056, A320451, A321718, A321719, A321722, A321733, A321734, A321735, A321736, A321739.

Sequence in context: A000255 A318912 A121580 * A224345 A081367 A156171

Adjacent sequences: A321729 A321730 A321731 * A321733 A321734 A321735

Keyword

nonn

Author

Gus Wiseman, Nov 18 2018

Extensions

a(7) onwards from Ludovic Schwob, Apr 03 2024

Status

approved