n阶矩阵中填入0,1,使得每行1的个数相等,每列1的个数两两不同的填法总数

由于每行1的个数相等，不妨设每行1的个数为 r ,在 n阶矩阵 A的 n 列中任取 $i_1,i_2,---,i_k$ 这 k 列,如果这 k 列的每一列都有 m 个1,现在来计算满足这一条件的填法总数.

考察这 k 列和 A 每一行交叉的 k 个位置,如果这 k 个位置有 t 个1,那么使这 k 个位置有 t 个1的取法总数是

$C_t^k$ ,每一种取法的0和1在 k 个交叉位置形成一个位模式,如果把所有取法按其对应的位模式的字典序编号为 $1,2,---,C_t^k$ ，则每一种取法在 n 行中出现的次数(非负整数)依次记为 ${\alpha }_{t1},{\alpha }_{t2},---,{\alpha }_{t{C}_t^k}$ ,对于每一种取法，由于在 k 个交叉位置有 t 个1,故同一行剩下 $n-k$ 个位置还可以放置 $r-t$ 个1，方法总数为 $C_{r-t}^{n-k}$ 故由乘法原理，在不考虑 ${\alpha }_{t1},{\alpha }_{t2},---,{\alpha }_{t{C}_t^k}$ 在 n 行中的具体分布的情况下所有的 $C_t^k$ 种取法对应的矩阵 A 的填法总数是 $(C_{r-t}^{n-k}{)}^{{\alpha }_{t1}+{\alpha }_{t2}+---+{\alpha }_{t{C}_t^k}}$ 再用乘法原理对 t 计算连乘积得到不考虑 ${\alpha }_{t1},{\alpha }_{t2},---,{\alpha }_{t{C}_t^k}$ 在 n 行中的具体分布的情况下所有可能的 t 对应的所有可能得取法对应的矩阵 A 的填法总数

${\prod }_{t=0}^{min(r,k)}(C_{r-t}^{n-k}{)}^{{\alpha }_{t1}+{\alpha }_{t2}+---+{\alpha }_{t{C}_t^k}}$

注意我们有约束关系 ${\sum }_{t=0}^k{\alpha }_{t1}+{\alpha }_{t2}+---+{\alpha }_{t{C}_t^k}=n$

以及 ${\sum }_{t=1}^k{\alpha }_{t1}^{(i)}+{\alpha }_{t2}^{(i)}+---+{\alpha }_{t{C}_{t-1}^{k-1}}^{(i)}=m$ $1<=i<=k$

(把 k 个交叉位置从左至右第 i 个位置上的数取1,剩下 $k-1$ 个位置取 $t-1$ 个1对应的每一种取法在

${\alpha }_{t1},{\alpha }_{t2},---,{\alpha }_{t{C}_t^k}$ 中对应的数依次记为 ${\alpha }_{t1}^{(i)},{\alpha }_{t2}^{(i)},---,{\alpha }_{t{C}_{t-1}^{k-1}}^{(i)}$ 约束关系 ${\sum }_{t=1}^k{\alpha }_{t1}^{(i)}+{\alpha }_{t2}^{(i)}+---+{\alpha }_{t{C}_{t-1}^{k-1}}^{(i)}=m$ 即指所有可能的 t 对应的所有可能得取法对应的矩阵 A 的填法中

$i_1,i_2,---,i_k$ 中的第 i 列的1的个数为m)

如果再考虑 ${\alpha }_{t1},{\alpha }_{t2},---,{\alpha }_{t{C}_t^k}$ ($0<=t<=k$) 在 n 行中的具体分布,由组合数公式所有的分布可能数

为 $\frac{n!}{{\prod }_{t=0}^k{\alpha }_{t1}!{\alpha }_{t2}!---{\alpha }_{t{C}_t^k}!}$ 从而由乘法原理所有可能的 t 对应的所有可能得取法对应的矩阵 A 的填法总数是

$\frac{n!}{{\prod }_{t=0}^k{\alpha }_{t1}!{\alpha }_{t2}!---{\alpha }_{t{C}_t^k}!}{\prod }_{t=0}^{min(r,k)}(C_{r-t}^{n-k}{)}^{{\alpha }_{t1}+{\alpha }_{t2}+---+{\alpha }_{t{C}_t^k}}$

对满足约束条件的所有可能的 ${\alpha }_{t1},{\alpha }_{t2},---,{\alpha }_{t{C}_t^k}$ 求和得到矩阵 A的 $i_1,i_2,---,i_k$ 这 k 列的每一列都有 m 个1的所有可能的填法总数 $f(n,m,k,r)=$

${\sum }_{{\sum }_{t=0}^k{\alpha }_{t1}+{\alpha }_{t2}+---+{\alpha }_{t{C}_t^k}=n,{\sum }_{t=1}^k{\alpha }_{t1}^{(i)}+{\alpha }_{t2}^{(i)}+---+{\alpha }_{t{C}_{t-1}^{k-1}}^{(i)}=m,1<=i<=k}\frac{n!}{{\prod }_{t=0}^k{\alpha }_{t1}!{\alpha }_{t2}!---{\alpha }_{t{C}_t^k}!}{\prod }_{t=0}^{min(r,k)}(C_{r-t}^{n-k}{)}^{{\alpha }_{t1}+{\alpha }_{t2}+---+{\alpha }_{t{C}_t^k}}$

下面令 $B_n$ 表示{$1,2,---,n$}的所有二元子集构成的集合, $B_{n,k}$ 表示 $B_n$ 的所有 k 元子集构成的集合

令 $A_{ij}$ 表示 A 满足每行1的个数为 r ,且i,j列1的个数不等的填法的集合

由容斥原理 $|{\bigcap }_{1<=i<j<=n}{A}_{ij}|=|S|+{\sum }_{k=1}^{C_2^n}(-1{)}^k{\sum }_{((i_1,j_1),(i_2,j_2),---,(i_k,j_k))\in B_{n,k}}|\overline{A_{i_1{j}_1}}\cap \overline{A_{i_2{j}_2}}\cap ---\cap \overline{A_{i_k{j}_k}}|$

$=|S|+{\sum }_{k=2}^n{\sum }_{1<=i_1<---<i_k<=n}E(i_1,i_2,---,i_k)|B(i_1,i_2,---,i_k)|$

其中 $S$ 表示向矩阵 A 填入0,1使得每行1的个数为 r 的填法集合 $|S|=(C_r^n{)}^n$

$B(i_1,i_2,---,i_k)$ 表示 $i_1,i_2,---,i_k$ 这k列1的个数相等的填法集合

下面来决定系数 $E(i_1,i_2,---i_k)$

令 $D_{k,f}$ 表示{ $i_1,i_2,---,i_k$ }的 $f$ 元子集构成的集合

$D_k^m$ 表示 $D_{k,2}$ 的所有 m 元子集构成的集合 $|D_k^m|=C_m^{C_2^k}$

$B_{i_t}$ ( $1<=t<=k$ ) 表示 $D_k^m$ 中包含包含 $i_t$ 的二元组的元素构成的集合,则由容斥原理

$|B_{i_1}\cap ---B_{i_k}|=|D_k^m|+{\sum }_{s=1}^k(-1{)}^s{\sum }_{(t_1,t_2,---,t_s)\in D_{k,s}}|\overline{B_{t_1}}\cap ---\cap \overline{B_{t_s}}|$

其中 $\overline{B_{t_1}}\cap ---\cap \overline{B_{t_s}}$表示 $D_k^m$ 中不包含包含 $i_t$( $1<=t<=s$ ) 的二元组的元素构成的集合

实际上 $\overline{B_{t_1}}\cap ---\cap \overline{B_{t_s}}$ 就是{ $i_1,i_2,---i_k$ }-{ $t_1,t_2,---,t_s$ }的二元子集构成的集合的全部m元子集构成的集合，其大小即为 $C_m^{C_2^{r-s}}$ 因此 $|\overline{B_{t_1}}\cap ---\cap \overline{B_{t_s}}|=C_m^{C_2^{r-s}}$

故 $|B_{i_1}\cap ---B_{i_k}|=C_m^{C_2^k}+{\sum }_{s=1}^{min(r,k)}(-1{)}^s{C}_s^k{C}_m^{C_2^{r-s}}$

从而 $E(i_1,i_2,---i_k)={\sum }_{m=1}^{C_2^k}(-1{)}^m|B_{i_1}\cap ---B_{i_k}|={\sum }_{m=1}^{C_2^k}(-1{)}^m(C_m^{C_2^k}+{\sum }_{s=1}^{min(r,k)}(-1{)}^s{C}_s^k{C}_m^{C_2^{r-s}})$

由前文结论 $|B(i_1,i_2,---,i_k)|={\sum }_{m=0}^{\left[\frac{nr}{k}\right]}f(n,m,k,r)$

求和上限是因为约束条件 $km<=rn$

于是我们有 $g(n,r)=|{\bigcap }_{1<=i<j<=n}{A}_{ij}|$ $=(C_r^n{)}^n+{\sum }_{k=2}^n{\sum }_{1<=i_1<---<i_k<=n}({\sum }_{m=1}^{C_2^k}(-1{)}^m(C_m^{C_2^k}+{\sum }_{s=1}^{min(r,k)}(-1{)}^s{C}_s^k{C}_m^{C_2^{r-s}})){\sum }_{m=0}^{\left[\frac{nr}{k}\right]}f(n,m,k,r)$

$=(C_r^n{)}^n+$ ${\sum }_{k=2}^n{C}_k^n({\sum }_{m=1}^{C_2^k}(-1{)}^m(C_m^{C_2^k}+{\sum }_{s=1}^{min(r,k)}(-1{)}^s{C}_s^k{C}_m^{C_2^{r-s}})){\sum }_{m=0}^{\left[\frac{nr}{k}\right]}f(n,m,k,r)$

最后,由于填入0,1后 A 的每一列1的个数都不同且每一列1的个数在0和n之间(包括0和n)又因为 A 有n列所以 A 的各列1的个数构成的集合实际上为 $A_{i,n}$ = { $0,1,---,n$ }-{ $i$ } ( $0<=i<=n$ )

于是所有可能的 r 的集合 $E_n$ = { $\frac{{\sum }_{x\in A_{i,n}}x}{n}$ | $\frac{{\sum }_{x\in A_{i,n}}x}{n}$ 为整数, $0<=i<=n$ }

从而n阶矩阵 A 中填入0,1,使得每行1的个数相等,每列1的个数两两不同的填法总数

$h(n)={\sum }_{r\in E_n}g(n,r)$

事实上,通过简单的计算可以发现设 $k$为正整数 $n=2k$ 时 $r=k$ $n=2k+1$ 时 $r=k$ 或 $k+1$

故 $n=2k$ 时 $E_n=${ $k$ } $n=2k+1$ 时 $E_n$ = { $k,k+1$ }