Enumerating contingency tables via random permanents
Abstract
Given m positive integers R=(ri), n positive integers C=(cj) such that sum ri = sum cj =N, and mn non-negative weights W=(wij), we consider the total weight T=T(R, C; W) of non-negative integer matrices (contingency tables) D=(dij) with the row sums ri, column sums cj, and the weight of D equal to prod wijdij. We present a randomized algorithm of a polynomial in N complexity which computes a number T'=T'(R,C; W) such that T' < T < alpha(R, C) T' where alpha(R,C) = minprod ri! ri-ri, prod cj! cj-cj NN/N!. In many cases, ln T' provides an asymptotically accurate estimate of ln T. The idea of the algorithm is to express T as the expectation of the permanent of an N x N random matrix with exponentially distributed entries and approximate the expectation by the integral T' of an efficiently computable log-concave function on Rmn. Applications to counting integer flows in graphs are also discussed.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.