Linear Reweighted Regularization Algorithms for Graph Matching Problem
Rongxuan Li
TL;DR
This work addresses graph matching, a challenging NP-hard problem, by relaxing permutation constraints to a convex set and introducing a linear reweighted regularization that preserves convexity. The core idea minimizes $f(X) = ||A X - X B||_F^2$ over a convex domain while adding a linear term $h_W(X)$ with iteratively updated weights, solved via a projected gradient method and accelerated projection onto the doubly stochastic set. A convergence theory shows that, under suitable conditions on the iterate proximity and regularization strength, the method converges to the graph matching solution; a practical variant updates weights and solves a sequence of convex subproblems, yielding sparse, high-quality solutions. The approach is demonstrated on synthetic networks and shapes, highlighting efficient sparsity-inducing behavior and potential for scalable graph matching in network alignment and computer vision tasks. Overall, the method provides a convexity-preserving, computationally efficient framework for sparse graph matching with broad applicability in pattern recognition and geometry processing.
Abstract
The graph matching problem is a significant special case of the Quadratic Assignment Problem, with extensive applications in pattern recognition, computer vision, protein alignments and related fields. As the problem is NP-hard, relaxation and regularization techniques are frequently employed to improve tractability. However, most existing regularization terms are nonconvex, posing optimization challenges. In this paper, we propose a linear reweighted regularizer framework for solving the relaxed graph matching problem, preserving the convexity of the formulation. By solving a sequence of relaxed problems with the linear reweighted regularization term, one can obtain a sparse solution that, under certain conditions, theoretically aligns with the original graph matching problem's solution. Furthermore, we present a practical version of the algorithm by incorporating the projected gradient method. The proposed framework is applied to synthetic instances, demonstrating promising numerical results.
