Improved Iteration Complexities for Overconstrained p-Norm Regression
Abstract
In this paper we obtain improved iteration complexities for solving p regression. We provide methods which given any full-rank A ∈ Rn × d with n ≥ d, b ∈ Rn, and p ≥ 2 solve x ∈ Rd \|A x - b\|p to high precision in time dominated by that of solving Op(dp-23p-2) linear systems in A D A for positive diagonal matrices D. This improves upon the previous best iteration complexity of Op(np-23p-2) (Adil, Kyng, Peng, Sachdeva 2019). As a corollary, we obtain an O(d1/3ε-2/3) iteration complexity for approximate ∞ regression. Further, for q ∈ (1, 2] and dual norm q = p/(p-1) we provide an algorithm that solves q regression in O(dp-22p-2) iterations. To obtain this result we analyze row reweightings (closely inspired by p-norm Lewis weights) which allow a closer connection between 2 and p regression. We provide adaptations of two different iterative optimization frameworks which leverage this connection and yield our results. The first framework is based on iterative refinement and multiplicative weights based width reduction and the second framework is based on highly smooth acceleration. Both approaches yield Op(dp-23p-2) iteration methods but the second has a polynomial dependence on p (as opposed to the exponential dependence of the first algorithm) and provides a new alternative to the previous state-of-the-art methods for p regression for large p.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.