370 likes | 388 Views
This paper presents a nearly optimal deterministic algorithm for volume estimation in the oracle model. It introduces a new construction for "easy-to-enumerate" thin-covering lattices. The algorithm reduces volume estimation to counting lattice points in a well-calibrated lattice.
E N D
Faster Deterministic Volume Estimation in the Oracle ModelviaThin Lattice Coverings Daniel Dadush Centrum Wiskunde & Informatica
Volume Estimation Problem Given convex body and factor , compute s.t. . Answers: is in ? given by a membership oracle.
Volume Estimation Problem Given convex body and factor , compute s.t. . Assume appropriate sandwiching guarantees.
Volume Estimation: Known Bounds Bárány, Füredi 86+88: Any deterministic-approximation uses at least membership queries for . Dyer, Frieze, Kannan 91: Randomized polynomial time-approximation algorithm.
Volume Estimation: Known Bounds Bárány, Füredi 86+88: Any deterministic-approximation uses at least membership queries for . Dyer, Frieze, Kannan 91: Randomized polynomial time-approximation algorithm. Open Problem: Deterministic PTAS for explicit polytopes?
Volume Estimation: Known Bounds Bárány, Füredi 86+88: Any deterministic-approximation uses at least membership queries for . Dyer, Frieze, Kannan 91: Randomized polynomial time-approximation algorithm. More modest question: Can we match the oracle model lower bounds?
Main Result Theorem[D. 14+]: Deterministic -approximation using -time and -space. Close toBárány & Füredilower bound.
Deterministic -approximations symmetric asymmetric
Lattices A lattice is all integral combinations of a basis .
Lattices A lattice is all integral combinations of a basis . The determinant of is.
Lattices A lattice is all integral combinations of a basis . Every tiling domain has volume .
High Level Idea Reduce approximating to counting lattice points in a ``well-calibrated’’ lattice . volume is
High Level Algorithm Input:convex body, approx factor .
High Level Algorithm Compute point about which is ``approximately symmetric’’: Treat above as optimization problem.
High Level Algorithm Compute point about which is ``approximately symmetric’’: Build using an “evolving net” over , evaluatingthe objective using volume algorithm for symmetric bodies.
High Level Algorithm Build “thin” lattice which covers space w.r.t. : b. a.
High Level Algorithm Build “thin” lattice which covers space w.r.t. : c.Enumeration over requires-space & -time.
High Level Algorithm Compute by enumeration. Return .
Algorithm Recap Find approximate center of symmetry . Build easy to enumerate thin covering lattice . Enumerate lattice points in blowup of .
Algorithm Recap Find approximate center of symmetry . Build easy to enumerate thin covering lattice . Enumerate lattice points in blowup of .
How to build Want easy to enumerate thin covering lattice .
How to build A good start:Pick basis of to be axes of an M-ellipsoid of .
Milman’s Ellipsoid is an M-Ellipsoid of if translates of cover and vice versa.
Milman’s Ellipsoid translates of suffice to cover .
Milman’s Ellipsoid translates of suffice to cover .
Milman’s Ellipsoid Theorem[D., Vempala 13]:Can construct M-ellipsoid in deterministic-time and -space.
How to build A good start:Pick basis of to be axes of an M-ellipsoid of . Get:& “easy to enumerate” .
How to build A good start:Pick basis of to be axes of an M-ellipsoid of . Don’t get: covering property.
How to build A good start:Pick basis of to be axes of an M-ellipsoid of . Idea: Find lattice “close” to satisfying covering property.
How to build Sparsify : Take “random” sublattice to remove vectors in .
How to build Sparsify : Take “random” sublattice to remove vectors in .
How to build Densify Greedily add cosets of that miss until no longer possible [Rogers 50].
How to build Densify Greedily add cosets of that miss until no longer possible [Rogers 50].
How to build Densify Greedily add cosets of that miss until no longer possible [Rogers 50].
How to build Return as final lattice.
Summary Nearly optimal deterministic algorithm for volume estimation in the oracle model. New construction of “easy to enumerate” thin-covering lattices. Open Questions Ignoring efficiency, can we match query lower bound? Other applications of thin-covering lattices?
-Net Barrier: Lemma: For any convex body and s.t. we have . Theorem [Barvinok 12]: Take small enough, then for a symmetric convex body , there exists a polytope with at most vertices. Question: Can we compute a as above using only queries?