中文

通过随机积和式枚举列联表

组合数学 2007-05-23 v2 概率论

摘要

给定mm个正整数R=(ri)R=(r_i)nn个正整数C=(cj)C=(c_j)满足ri=cj=N\sum r_i = \sum c_j =N,以及mnmn个非负权W=(wij)W=(w_{ij}),我们考虑具有行和rir_i、列和cjc_j且权为wijdij\prod w_{ij}^{d_{ij}}的非负整数矩阵(列联表)D=(dij)D=(d_{ij})的总权T=T(R,C;W)T=T(R, C; W)。我们给出一个关于NN的多项式复杂度的随机算法,计算数T=T(R,C;W)T'=T'(R,C; W)使得T<T<α(R,C)TT' < T < \alpha(R, C) T',其中α(R,C)=min{ri!riri,cj!cjcj}NN/N!\alpha(R,C) = \min\{\prod r_i! r_i^{-r_i}, \prod c_j! c_j^{-c_j}\} N^N/N!。在许多情形下,lnT\ln T'给出lnT\ln T的渐近准确估计。该算法的思想是将TT表示为一个N×NN \times N具指数分布元素的随机矩阵的积和式的期望,并用RmnR^{mn}上可有效计算的对数凹函数之积分TT'逼近该期望。也讨论了在图中计数整数流的应用。

关键词

引用

@article{arxiv.math/0511596,
  title  = {Enumerating contingency tables via random permanents},
  author = {Alexander Barvinok},
  journal= {arXiv preprint arXiv:math/0511596},
  year   = {2007}
}

备注

19 pages, bounds are sharpened, references are added