TR2025-098

Recursive McCormick Linearization of Multilinear Programs


    •  Raghunathan, A., Cardonha, C., Bergman, D., Nohra, C.J., "Recursive McCormick Linearization of Multilinear Programs", INFORMS J Computing, June 2025.
      BibTeX TR2025-098 PDF
      • @article{Raghunathan2025jun,
      • author = {Raghunathan, Arvind and Cardonha, Carlos and Bergman, David and Nohra, Carlos J.},
      • title = {{Recursive McCormick Linearization of Multilinear Programs}},
      • journal = {INFORMS J Computing},
      • year = 2025,
      • month = jun,
      • url = {https://www.merl.com/publications/TR2025-098}
      • }
  • MERL Contact:
  • Research Area:

    Optimization

Abstract:

Linear programming (LP) relaxations are widely employed in exact solution methods for multilinear programs (MLPs). These relaxations can be obtained by using Recursive McCormick Linearizations (RMLs), where an MLP is linearized by iteratively substituting bilinear products with artificial variables and constraints. This article introduces a systematic approach to identifying RMLs. We focus on identifying RMLs with a small number of artificial variables and strong LP bounds. We present a novel mechanism for representing all the possible RMLs, which we use to design an exact mixed-integer programming (MIP) formulation to identify minimum-size RMLs; this problem is NP-hard in general, but we show that it is fixed-parameter tractable if each monomial is composed of at most three variables. Moreover, we explore the structural properties of our formulation to derive an exact MIP model that identifies RMLs of a given size with the best-possible LP relaxation bound. We test our algorithms by conducting numerical experiments on a large collection of MLPs. Numerical results indicate that the RMLs obtained with our algorithms can be significantly smaller than those derived from heuristic or greedy approaches, leading, in many cases, to tighter LP relaxation bounds. Moreover, our linearization strategies can be used to reformulate MLPs as QCPs, which can then be efficiently solved using state-of-the-art solvers for QCPs. This QCP-based solution approach is highly beneficial for hard MLP instances.