Solving Matrix Games with Near-Optimal Matvec Complexity
Ishani Karmarkar, Liam O'Carroll, Aaron Sidford
TL;DR
The paper tackles computing $oldsymbol{ extepsilon}$-approximate Nash equilibria in two-player matrix games with bilinear payoffs under simple strategy sets, focusing on matvec-efficient algorithms. It introduces a prox multi-point outer-loop framework with a dyadic, multi-center approach and an amortized analysis that bounds total movement across iterations, enabling a $ ilde{O}(oldsymbol{ extepsilon}^{-2/3})$ matvec complexity for both $ ext{ell}_1 ext{-ell}_1$ and $ ext{ell}_2 ext{-ell}_1$ settings. The approach combines a dynamic bisection procedure for regularization, a smooth-until-proven-guilty inner solver, and a matrix-approximation path to learn $A$ efficiently, achieving near-optimal performance up to polylogarithmic factors. The framework unifies and extends prior methods, providing a scalable path toward solving broad classes of structured convex-concave problems with limited access to the matrix via matvecs, and it includes detailed robustness and stability analyses that support its practical applicability. Overall, the work significantly advances the matvec efficiency frontier for bilinear matrix games and offers techniques likely applicable to broader monotone-operator problems in optimization and game theory.
Abstract
We study the problem of computing an $ε$-approximate Nash equilibrium of a two-player, bilinear game with a bounded payoff matrix $A \in \mathbb{R}^{m \times n}$, when the players' strategies are constrained to lie in simple sets. We provide algorithms which solve this problem in $\tilde{O}(ε^{-2/3})$ matrix-vector multiplies (matvecs) in two well-studied cases: $\ell_1$-$\ell_1$ (or zero-sum) games, where the players' strategies are both in the probability simplex, and $\ell_2$-$\ell_1$ games (encompassing hard-margin SVMs), where the players' strategies are in the unit Euclidean ball and probability simplex respectively. These results improve upon the previous state-of-the-art complexities of $\tilde{O}(ε^{-8/9})$ for $\ell_1$-$\ell_1$ and $\tilde{O}(ε^{-7/9})$ for $\ell_2$-$\ell_1$ due to [KOS '25]. In both settings our results are nearly-optimal as they match lower bounds of [KS '25] up to polylogarithmic factors.
