Convex volume approximation
In the analysis of algorithms, several authors have studied the computation of the volume of high-dimensional convex bodies, a problem that can also be used to model many other problems in combinatorial enumeration. Often these works use a black box model of computation in which the input is given by a subroutine for testing whether a point is inside or outside of the convex body, rather than by an explicit listing of the vertices or faces of a convex polytope. It is known that, in this model, no deterministic algorithm can achieve an accurate approximation,[1][2] and even for an explicit listing of faces or vertices the problem is #P-hard.[3] However, a joint work by Martin Dyer, Alan M. Frieze and Ravindran Kannan provided a randomized polynomial time approximation scheme for the problem, providing a sharp contrast between the capabilities of randomized and deterministic algorithms.[4]
The main result of the paper is a randomized algorithm for finding an approximation to the volume of a convex body in -dimensional Euclidean space by assuming the existence of a membership oracle. The algorithm takes time bounded by a polynomial in , the dimension of and . The algorithm combines two ideas:
- By using a Markov chain Monte Carlo (MCMC) method, it is possible to generate points that are nearly uniformly randomly distributed within a given convex body. The basic scheme of the algorithm is a nearly uniform sampling from within by placing a grid consisting of -dimensional cubes and doing a random walk over these cubes. By using the theory of rapidly mixing Markov chains, they show that it takes a polynomial time for the random walk to settle down to being a nearly uniform distribution.[4]
- By using rejection sampling, it is possible to compare the volumes of two convex bodies, one nested within another, when their volumes are within a small factor of each other. The basic idea is to generate random points within the outer of the two bodies, and to count how often those points are also within the inner body.
The given convex body can be approximated by a sequence of nested bodies, eventually reaching one of known volume (a hypersphere), with this approach used to estimate the factor by which the volume changes at each step of this sequence. Multiplying these factors gives the approximate volume of the original body.
This work earned its authors the 1991 Fulkerson Prize.[5]
Improvements
[edit]Although the time for this algorithm is polynomial, it has a high exponent. Subsequent authors improved the running time of this method by providing more quickly mixing Markov chains for the same problem.[6][7][8][9]
Generalizations
[edit]The polynomial-time approximability result has been generalized to more complex structures such as the union and intersection of objects.[10] This relates to Klee's measure problem.
References
[edit]- ^ Elekes, G. (1986), "A geometric inequality and the complexity of computing volume", Discrete and Computational Geometry, 1 (4): 289–292, doi:10.1007/BF02187701, MR 0866364
- ^ Bárány, Imre; Füredi, Zoltán (1987), "Computing the volume is difficult", Discrete and Computational Geometry, 2 (4): 319–326, doi:10.1007/BF02187886, hdl:1813/8572, MR 0911186
- ^ Dyer, Martin; Frieze, Alan (1988), "On the complexity of computing the volume of a polyhedron", SIAM Journal on Computing, 17 (5): 967–974, doi:10.1137/0217060, MR 0961051
- ^ a b Dyer, Martin; Frieze, Alan; Kannan, Ravi (1991), "A random polynomial-time algorithm for approximating the volume of convex bodies", Journal of the ACM, 38 (1): 1–17, doi:10.1145/102782.102783, MR 1095916, S2CID 13268711
- ^ Fulkerson Prize winners, American Mathematical Society, retrieved 2017-08-03.
- ^ Applegate, David; Kannan, Ravi (1991), "Sampling and Integration of Near Log-concave Functions", Proceedings of the Twenty-Third Annual ACM Symposium on Theory of Computing (STOC '91), New York, NY, USA: ACM, pp. 156–163, doi:10.1145/103418.103439, ISBN 978-0-89791-397-3, S2CID 15432190
- ^ Kannan, Ravi; Lovász, László; Simonovits, Miklós (1997), "Random walks and an volume algorithm for convex bodies", Random Structures & Algorithms, 11 (1): 1–50, doi:10.1002/(SICI)1098-2418(199708)11:1<1::AID-RSA1>3.0.CO;2-X, MR 1608200
- ^ Lovász, L.; Simonovits, M. (1993), "Random walks in a convex body and an improved volume algorithm", Random Structures & Algorithms, 4 (4): 359–412, doi:10.1002/rsa.3240040402, MR 1238906
- ^ Lovász, L.; Vempala, Santosh (2006), "Simulated annealing in convex bodies and an volume algorithm", Journal of Computer and System Sciences, 72 (2): 392–417, doi:10.1016/j.jcss.2005.08.004, MR 2205290
- ^ Bringmann, Karl; Friedrich, Tobias (2010-08-01). "Approximating the volume of unions and intersections of high-dimensional geometric objects". Computational Geometry. 43 (6): 601–610. arXiv:0809.0835. doi:10.1016/j.comgeo.2010.03.004. hdl:11858/00-001M-0000-000F-1603-4. ISSN 0925-7721. S2CID 5930593.
