Sampling binary contingency tables with a greedy start

  • Authors:
  • Ivona Bezáková;Nayantara Bhatnagar;Eric Vigoda

  • Affiliations:
  • University of Chicago, Chicago, IL;Georgia Institute of Technology, Atlanta, GA;Georgia Institute of Technology, Atlanta, GA

  • Venue:
  • SODA '06 Proceedings of the seventeenth annual ACM-SIAM symposium on Discrete algorithm
  • Year:
  • 2006

Quantified Score

Hi-index 0.00

Visualization

Abstract

We study the problem of counting and randomly sampling binary contingency tables. For given row and column sums, we are interested in approximately counting (or sampling) 0/1 n x m matrices with the specified row/column sums. We present a simulated annealing algorithm with running time O((nm)2D3dmaxlog5(n+m)) for any row/column sums where D is the number of non-zero entries and dmax is the maximum row/column sum. This is the first algorithm to directly solve binary contingency tables for all row/column sums. Previous work reduced the problem to the permanent, or restricted attention to row/column sums that are close to regular. The interesting aspect of our simulated annealing algorithm is that it starts at a non-trivial instance, whose solution relies on the existence of short alternating paths in the graph constructed by a particular Greedy algorithm.