[论文解读] Asymptotic enumeration of integer matrices with constant row and column sums

E. Rodney Canfield, Brendan D. McKay|arXiv (Cornell University)|Mar 20, 2007

Advanced Combinatorial Mathematics参考文献 17被引用 9

一句话总结

本文建立了当 $m,n \to \infty$ 且平均元素 $\lambda = s/n = t/m$ 不太小时，具有指定常数行和与列和的 $m \times n$ 整数矩阵数量的精确渐近公式。通过鞍点法与多面体几何，证明了数量 $M(m,s;n,t)$ 渐近为 $G(m,s;n,t) \cdot \exp(1/2 + O(n^{-b}))$，其中 $G$ 是基于多项式系数的组合近似，并给出了伯克霍夫多面体的新渐近体积公式。

ABSTRACT

Let s,t,m,n be positive integers such that sm=tn. Let M(m,s;n,t) be the number of m x n matrices over {0,1,2,...} with each row summing to s and each column summing to t. Equivalently, M(m,s;n,t) counts 2-way contingency tables of order m x n such that the row marginal sums are all s and the column marginal sums are all t. A third equivalent description is that M(m,s;n,t) is the number of semiregular labelled bipartite multigraphs with m vertices of degree s and n vertices of degree t. When m=n and s=t such matrices are also referred to as n x n magic squares with line sums equal to t. We prove a precise asymptotic formula for M(m,s;n,t) which is valid over a range of (m,s;n,t) in which m,n become infinite while remaining approximately equal and the average entry is not too small. This range includes the case where m/n, n/m, s/n and t/m are bounded from below.

研究动机与目标

• 推导 $m \times n$ 整数矩阵数量的精确渐近公式，其行和固定为 $s$，列和固定为 $t$，且满足 $ms = nt$。
• 将以往仅适用于小规模或稀疏情况的渐近结果，扩展至包含稠密矩阵在内的广泛参数范围。
• 为伯克霍夫多面体（即双随机矩阵的集合）提供首个渐近体积公式。
• 严格验证古德的启发式近似 $G(m,s;n,t)$ 在广阔参数范围内作为 $M(m,s;n,t)$ 的近似值的准确性。

提出的方法

• 对具有固定行和与列和的整数矩阵的多元生成函数应用鞍点法。
• 通过拉普拉斯型近似，定义一个区域 $\mathcal{R}'$，使得积分的主要贡献来自该区域。
• 在 $\lambda \geq n^5$ 的稠密情形下，利用运输多面体的Ehrhart多项式分析格点计数，其中 $\lambda = s/n = t/m$。
• 通过斯坦利的Ehrhart多项式理论，建立 $M(m,s;n,t)$ 与稀释运输多面体 $q\mathcal{P}$ 中整数点数量之间的联系。
• 证明不同 $\lambda$ 值下计数的比值收敛于某个 2 的幂，从而实现参数范围内的插值。
• 结合稀疏情形（$\lambda = O(n^5)$）与稠密情形（$\lambda \geq n^5$）的结果，证明完整的渐近公式。

实验结果

研究问题

• RQ1 当 $m,n \to \infty$ 时，具有常数行和 $s$ 与列和 $t$ 的 $m \times n$ 整数矩阵数量的渐近增长速率是什么？
• RQ2 古德的启发式近似 $G(m,s;n,t)$ 在不同参数范围内对 $M(m,s;n,t)$ 的估计有多准确？
• RQ3 能否利用此框架对伯克霍夫多面体的体积进行渐近估计？
• RQ4 在何种条件下，近似 $G(m,s;n,t)$ 与真实计数 $M(m,s;n,t)$ 的误差在常数因子范围内？
• RQ5 运输多面体及其Ehrhart多项式的结构如何促进列联表渐近计数的分析？

主要发现

• 数量 $M(m,s;n,t)$ 渐近为 $G(m,s;n,t) \cdot \exp(1/2 + O(n^{-b}))$，其中 $G(m,s;n,t)$ 是涉及二项式系数的组合表达式。
• 该公式在条件 $\frac{(1+2\lambda)^2}{4\lambda(1+\lambda)} \left(1 + \frac{5m}{6n} + \frac{5n}{6m}\right) \leq a \log n$ 下成立，其中常数 $a,b > 0$ 满足 $a+b < 1/2$，该条件确保了参数范围包含 $m/n$ 和 $s/n$ 远离零的情况。
• 当 $\lambda = O(n^5)$ 时，通过多变量积分表示的鞍点分析推导出渐近公式。
• 当 $\lambda \geq n^5$ 时，利用运输多面体的Ehrhart多项式将公式扩展，表明不同 $\lambda$ 值下计数的比值收敛于某个 2 的幂。
• 结果表明，古德的近似 $G(m,s;n,t)$ 极其精确，在广阔参数范围内与 $M(m,s;n,t)$ 的差异至多为常数因子。
• 作为推论，首次获得了伯克霍夫多面体的渐近体积公式，该公式在 $m=n$ 且 $s=t$ 时成立。