Improving the performance of standard solvers for quadratic 0-1 programs by a tight convex reformulation: The QCR method

  • Authors:

  • Alain Billionnet;Sourour Elloumi;Marie-Christine Plateau

  • Affiliations:

  • Laboratoire CEDRIC, Institut d'Informatique d'Entreprise, 18 allée Jean Rostand, F-91025 Evry, France;Laboratoire CEDRIC, Conservatoire National des Arts et Métiers, 292 rue Saint Martin, F-75141 Paris, France;Laboratoire CEDRIC, Conservatoire National des Arts et Métiers, 292 rue Saint Martin, F-75141 Paris, France

  • Venue:
  • Discrete Applied Mathematics
  • Year:
  • 2009

Quantified Score

Hi-index 0.05

Visualization

Abstract

Let (QP) be a 0-1 quadratic program which consists in minimizing a quadratic function subject to linear equality constraints. In this paper, we present QCR, a general method to reformulate (QP) into an equivalent 0-1 program with a convex quadratic objective function. The reformulated problem can then be efficiently solved by a classical branch-and-bound algorithm, based on continuous relaxation. This idea is already present in the literature and used in standard solvers such as CPLEX. Our objective in this work was to find a convex reformulation whose continuous relaxation bound is, moreover, as tight as possible. From this point of view, we show that QCR is optimal in a certain sense. State-of-the-art reformulation methods mainly operate a perturbation of the diagonal terms and are valid for any {0,1} vector. The innovation of QCR comes from the fact that the reformulation also uses the equality constraints and is valid on the feasible solution domain only. Hence, the superiority of QCR holds by construction. However, reformulation by QCR requires the solution of a semidefinite program which can be costly from the running time point of view. We carry out a computational experience on three different combinatorial optimization problems showing that the costly computational time of reformulation by QCR can however result in a drastically more efficient branch-and-bound phase. Moreover, our new approach is competitive with very specific methods applied to particular optimization problems.