Approximability of Symmetric Bimatrix Games and Related Experiments

In this work we present a simple quadratic formulation for the problem of computing Nash equilibria in symmetric bimatrix games, inspired by the well-known formulation of Mangasarian and Stone [26]. We exploit our formulation to shed light to the approximability of NE points. First we observe that any KKT point of this formulation (and indeed, any quadratic program) is also a stationary point, and vice versa. We then prove that any KKT point of the proposed formulation (is not necessarily itself, but) indicates a ( 0, in time polynomial in the size of the game and quasi-linear in 1/δ, exploiting Ye's algorithm for approximating KKT points of QPs [34]. This is (to our knowledge) the first polynomial time algorithm that constructs e-NE points for symmetric bimatrix games for any e close to 1/3. We extend our main result to (asymmetric) win lose games, as well as to games with maximum aggregate payoff either at most 1, or at least 5/3. To achieve this, we use a generalization of the Brown & von Neumann symmetrization technique [6] to the case of non-zero-sum games, which we prove that is approximation preserving. Finally, we present our experimental analysis of the proposed approximation and other quite interesting approximations for NE points in symmetric bimatrix games.