Counting Integral Points in Polytopes via Numerical Analysis of Contour Integration
Hiroshi Hirai,Ryunosuke Oshiro,Ken’ichiro Tanaka
DOI: https://doi.org/10.1287/moor.2019.0997
IF: 2.215
2020-05-01
Mathematics of Operations Research
Abstract:In this paper, we address the problem of counting integer points in a rational polytope described by <i>P</i>(<b><i>y</i></b>) = {<i><b>x</b></i> ∈ ℝ<sup><i>m</i></sup>: <i>A<b>x</b></i>= <i><b>y, x</b> ≥</i> 0}, where <i>A</i> is an <i>n</i> × <i>m</i> integer matrix and <i><b>y</b></i> is an <i>n</i>-dimensional integer vector. We study the Z transformation approach initiated in works by Brion and Vergne, Beck, and Lasserre and Zeron from the numerical analysis point of view and obtain a new algorithm on this problem. If <i>A</i> is nonnegative, then the number of integer points in <i>P</i>(<i><b>y</b></i>) can be computed in <span class="equationTd inline-formula"><math>O(poly(n,m,‖y‖∞)(‖y‖∞+1)n)</math></span> time and <span class="equationTd inline-formula"><math>O(poly(n,m,‖y‖∞))</math></span> space. This improves, in terms of space complexity, a naive DP algorithm with <span class="equationTd inline-formula"><math>O((‖y‖∞+1)n)</math></span>-size dynamic programming table. Our result is based on the standard error analysis of the numerical contour integration for the inverse Z transform and establishes a new type of an inclusion-exclusion formula for integer points in <i>P</i>(<i><b>y</b></i>).We apply our result to hypergraph <b><i>b</i></b>-matching and obtain a <span class="equationTd inline-formula"><math>O(poly(n,m,‖b‖∞)(‖b‖∞+1)(1−1/k)n)</math></span> time algorithm for counting <b><i>b</i></b>-matchings in a <i>k</i>-partite hypergraph with <i>n</i> vertices and <i>m</i> hyperedges. This result is viewed as a <b><i>b</i></b>-matching generalization of the classical result by Ryser for <i>k</i> = 2 and its multipartite extension by Björklund and Husfeldt.
mathematics, applied,operations research & management science