[論文レビュー] Concave Quadratic Cuts for Mixed-Integer Quadratic Problems

Jaehyun Park, Stephen Boyd|arXiv (Cornell University)|Oct 21, 2015

Optimization and Packing Problems参考文献 123被引用数 4

ひとこと要約

本稿では、整数格子上での有効な凹二次不等式を用いて、半定値プログラミング（SDP）緩和を強化する、混合整数二次プログラミング（MIQPs）のための新規カット平面法を提案する。これらのカットを追加することにより、非凸なQCQPsにおける最適値の下界が向上し、整数最小二乗問題およびマックスカット問題において顕著なギャップ低減が実証された。特に、ゲーマンズ＝ウィリアムソン境界を上回る27%までの最適性ギャップ改善が達成された。

ABSTRACT

The technique of semidefinite programming (SDP) relaxation can be used to obtain a nontrivial bound on the optimal value of a nonconvex quadratically constrained quadratic program (QCQP). We explore concave quadratic inequalities that hold for any vector in the integer lattice ${\bf Z}^n$, and show that adding these inequalities to a mixed-integer nonconvex QCQP can improve the SDP-based bound on the optimal value. This scheme is tested using several numerical problem instances of the max-cut problem and the integer least squares problem.

研究の動機と目的

• 非凸な混合整数二次制約付き二次プログラミング（MIQCQPs）の最適値の下界を向上させること。
• 整数格子 Z^n 上で有効な凹二次不等式に基づくカット平面法を考案すること。
• これらのカットを用いて、半定値プログラミング（SDP）緩和を強化し、NP困難な組合せ問題における境界品質を向上させること。
• 整数最小二乗問題およびマックスカット問題といったベンチマーク問題に対して、タイトな境界が不可欠であるため、その評価を実施すること。

提案手法

• 整数ベクトルの構造に基づき、すべての x ∈ Z^n に対して成り立つ凹二次不等式を導出する。
• A が負半定値である形式 x^T A x + b^T x + c ≤ 0 の有効な不等式を、すべての整数ベクトルが満たすように定式化する。
• これらの不等式をMIQCQPのSDP緩和にカット平面として統合し、問題の実行可能領域を変更せずに境界を改善する。
• ||a||_2 = √3 を満たすランダムなカットベクトル a の生成を用いて、有効なカットを同定・追加し、緩和をタイトにする。
• 制約の再定式化により、整数最小二乗問題およびマックスカット問題を含む、QCQPsの最小化および最大化形に本手法を適用する。
• 標準的なSDPソルバを用いて、追加カットを含むSDPを解き、新しい境界をゲーマンズ＝ウィリアムソン境界などの古典的SDP緩和と比較する。

実験結果

リサーチクエスチョン

• RQ1 整数格子上での有効な凹二次不等式を用いて、非凸なMIQCQPsのSDP緩和を強化できるか？
• RQ2 ゲーマンズ＝ウィリアムソン境界などの古典的SDP境界と比較して、これらのカットは最適性ギャップをどの程度低減できるか？
• RQ3 制約構造が異なるにもかかわらず、本手法は整数最小二乗問題およびマックスカット問題の両方でよりタイトな境界を達成できるか？
• RQ4 境界の向上に伴い、追加するカットの数と計算コストのトレードオフはどのようなものか？
• RQ5 本手法と高階のラッセルレ階層の間に構造的関係は存在するか？

主な発見

• 整数最小二乗問題において、提案手法によるカットは、古典的SDP境界と比較して最適性ギャップを最大27%まで低減した。
• マックスカット問題において、本手法は平均して12〜27%の最適性ギャップ低減を達成し、元のギャップと削減後のギャップの比 α は0.73から0.88の範囲に収まった。
• ||a||_2 = √3 を満たすカットは顕著な改善をもたらしたが、||a||_2 = √2 を満たすカットは測定可能な利益をもたらさなかった。
• テストした10個のマックスカットインスタンスすべてにおいて、本手法はゲーマンズ＝ウィリアムソン境界を上回るよりタイトな上界を生成した。特にG54300およびG54600で最も顕著な改善が観察された。
• 約1,000個のカットを追加した場合、平均的な解法時間は約1秒から約15〜23秒に増加し、境界品質と計算コストの明確なトレードオフが示された。
• 本手法は、不確定な二次カットがブール問題に対しても有効であることを示した。これは、本研究で検証されたケースを超えた幅広い応用可能性を示唆している。