A Branch-and-Bound Approach for Maximum Low-Diameter Dense Subgraph Problems
Yi Zhou, Chunyu Luo, Zhengren Wang, Zhang-Hua Fu
TL;DR
The paper defines the maximum low-diameter $f(\cdot)$-dense subgraph problem (M$f$DS), a connectivity-respecting generalization of dense-subgraph models. It proposes a practical decomposition-based framework that partitions the input graph into small subproblems, combined with a degeneracy or two-hop degeneracy ordering and a branch-and-bound solver augmented by a sorting-based upper bound (SortBound). A mixed integer programming model (MIP-$f$D) is provided to formalize the problem, and extensive experiments on 139 real-world graphs show that the decomposition-enhanced approach outperforms standalone MIP or pure BnB, solving roughly twice as many instances within one hour. The work demonstrates that exact solutions for large-scale M$f$DS problems are feasible in practice and highlights the potential for parallelization and broader applicability to other diameter-constrained dense-subgraph tasks.
Abstract
A graph with $n$ vertices is an $f(\cdot)$-dense graph if it has at least $f(n)$ edges, $f(\cdot)$ being a well-defined function. The notion $f(\cdot)$-dense graph encompasses various clique models like $γ$-quasi cliques, $k$-defective cliques, and dense cliques, arising in cohesive subgraph extraction applications. However, the $f(\cdot)$-dense graph may be disconnected or weakly connected. To conquer this, we study the problem of finding the largest $f(\cdot)$-dense subgraph with a diameter of at most two in the paper. Specifically, we present a decomposition-based branch-and-bound algorithm to optimally solve this problem. The key feature of the algorithm is a decomposition framework that breaks the graph into $n$ smaller subgraphs, allowing independent searches in each subgraph. We also introduce decomposition strategies including degeneracy and two-hop degeneracy orderings, alongside a branch-and-bound algorithm with a novel sorting-based upper bound to solve each subproblem. Worst-case complexity for each component is provided. Empirical results on 139 real-world graphs under two $f(\cdot)$ functions show our algorithm outperforms the MIP solver and pure branch-and-bound, solving nearly twice as many instances optimally within one hour.
