Mining for diamonds—Matrix generation algorithms for binary quadratically constrained quadratic problems

• Published in: Computers & operations research Vol. 142; p. 105735
• Main Authors: Bettiol, Enrico, Bomze, Immanuel, Létocart, Lucas, Rinaldi, Francesco, Traversi, Emiliano
• Format: Journal Article
• Language: English
• Published: New York Elsevier Ltd 01.06.2022; Pergamon Press Inc; Elsevier
• Subjects: Algorithms; Binary quadratic problem; Boolean; Boolean algebra; Boolean Quadric Polytope; Computer Science; Dantzig–Wolfe Reformulation; Decomposition; Diamonds; Discrete Mathematics; Equivalence; Operations Research; Polytopes; Semidefinite programming; Sparse problem; Sparsity; 90C20; Dantzig–Wolfe Reformulation; Binary quadratic problem; 90C26; 49M27; Sparse problem; 90C09; Boolean Quadric Polytope
• ISSN: 0305-0548; 1873-765X; 1873-765X; 0305-0548
• DOI: 10.1016/j.cor.2022.105735

In this paper, we consider binary quadratically constrained quadratic problems and propose a new approach to generate stronger bounds than the ones obtained using the Semidefinite Programming relaxation. The new relaxation is based on the Boolean Quadric Polytope and is solved via a Dantzig–Wolfe Reformulation in matrix space. For block-decomposable problems, we extend the relaxation and analyze the theoretical properties of this novel approach. If overlapping size of blocks is at most two (i.e., when the sparsity graph of any pair of intersecting blocks contains either a cut node or an induced diamond graph), we establish equivalence to the one based on the Boolean Quadric Polytope. We prove that this equivalence does not hold if the sparsity graph is not chordal and we conjecture that equivalence holds for any block structure with a chordal sparsity graph. The tailored decomposition algorithm in the matrix space is used for efficiently bounding sparsely structured problems. Preliminary numerical results show that the proposed approach yields very good bounds in reasonable time. •New relaxation based upon Boolean Quadric Polytope.•Solved via a Dantzig–Wolfe Reformulation in matrix space.•Tighter bounds than SDP-based ones.•Efficient bounds for sparse structures.