A Faster Cutting Plane Method and its Implications for Combinatorial and Convex Optimization
Yin Tat Lee, Aaron Sidford, Sam Chiu-wai Wong · 2015
In this paper we improve upon the running time for finding a point in a convex set given a separation oracle. In particular, given a separation oracle for a convex set K ⊂ Rnthat is contained in a box of radius R we show how to either compute a point in K or prove that K does not contain a ball of radius ϵ using an expected O(n log(nR/ϵ)) evaluations of the oracle and additional time O(n3logO(1)(nR/ϵ)). This matches the oracle complexity and improves upon the O(nω+1log(nR/ϵ)) additional time of the previous fastest algorithm achieved over 25 years ago by Vaidya [91] for the current value of the matrix multiplication constant w2log nM · EO + n3logO(1)nM) and O(n3log2n · EO + n4logO(1)n), improving upon the previous best of O((n4· EO + n5)logM) and O(n5· EO + n6) respectively. · Submodular Flow: n = |V|, m = |E|, C is the maximum edge cost in absolute value and U is maximum edge capacity in absolute value. We obtain a faster weakly polynomial running time of O(n2log nCU · EO + n3logO(1) nCU), improving upon the previous best of O(mn5log nU · EO) and O (n4h min {log C, log U}) from 15 years ago by a factor of Õ(n4). We also achieve faster strongly polynomial time algorithms as a consequence of our result on submodular minimization. · Matroid Intersection: n is the size of the ground set, r is the maximum size of independent sets, M is the maximum absolute value of element weight, Trankand Tindare the time for each rank and independence oracle query. We obtain a running time of O((nr log2nTrank+n3logO(1)n) log nM) and O((n2log nTind+n3logO(1)n) log nM), achieving the first quadratic bound on the query complexity for the independence and rank oracles. In the unweighted case, this is the first improvement since 1986 for independence oracle. · Semidefinite Programming: n is the number of constraints, m is the number of dimensions and S is the total number of non-zeros in the constraint matrices. We obtain a running time of O(n(n2+ mω+ S)), improving upon the previous best of Õ(n(nω+ mω+ S)) for the regime S is small.