[论文解读] Polynomial-Time Algorithms for Submodular Laplacian Systems

Kaito Fujii, Tasuku Soma+1 more|arXiv (Cornell University)|Mar 29, 2018

Graph theory and applications参考文献 13被引用 8

一句话总结

本文提出了多项式时间算法来求解子模拉普拉斯系统，这是经典拉普拉斯系统向子模变换的推广。通过利用椭球法和参数化最小 $st$-割计算，作者证明了精确系统 $L_F(\bm{x}) \ni \bm{b}$ 及其回归形式均可高效求解，将谱图理论统一扩展至有向图、超图和概率模型框架。

ABSTRACT

Let $G=(V,E)$ be an undirected graph, $L_G\in \mathbb{R}^{V imes V}$ be the associated Laplacian matrix, and $b \in \mathbb{R}^V$ be a vector. Solving the Laplacian system $L_G x = b$ has numerous applications in theoretical computer science, machine learning, and network analysis. Recently, the notion of the Laplacian operator $L_F:\mathbb{R}^V o 2^{\mathbb{R}^V}$ for a submodular transformation $F:2^V o \mathbb{R}_+^E$ was introduced, which can handle undirected graphs, directed graphs, hypergraphs, and joint distributions in a unified manner. In this study, we show that the submodular Laplacian system $L_F( x) i b$ can be solved in polynomial time. Furthermore, we also prove that even when the submodular Laplacian system has no solution, we can solve its regression form in polynomial time. Finally, we discuss potential applications of submodular Laplacian systems in machine learning and network analysis.

研究动机与目标

• 将拉普拉斯系统理论从无向图推广至子模变换。
• 建立子模拉普拉斯系统 $L_F(\bm{x}) \ni \bm{b}$ 在多项式时间内可解的可 tractability。
• 通过高效求解回归形式来处理无精确解存在的情形。
• 在单一子模拉普拉斯框架下统一处理无向图、有向图、超图和联合分布。

提出的方法

• 本文使用椭球法求解子模拉普拉斯系统 $L_F(\bm{x}) \ni \bm{b}$，证明其具有多项式时间可解性。
• 将问题表述为在拟模多面体上的凸优化问题，利用 $F$ 的核的 Birkhoff 表示。
• 对于回归形式，作者通过参数化最小 $st$-割算法最小化包含子模函数 $H$ 和二次项的复合目标函数。
• 该方法构建一个有向图 $G = (U \cup \{s,t\}, A \cup \bar{A})$，其容量通过 $\min\{\alpha - \bm{b}(u), 0\}$ 和 $\min\{-\alpha + \bm{b}(u), 0\}$ 定义，以建模参数化流。
• 应用 GGT 结构处理参数化流，实现所有临界 $\alpha$-值的计算，即最小割发生变化的点。
• 当给定 $\ker(F)$ 的 Birkhoff 表示时，该算法的时间复杂度为 $O(nm^2 \log(n^2/m))$，达到强多项式时间复杂度。

实验结果

研究问题

• RQ1 即使算子是多值且非线性的，子模拉普拉斯系统 $L_F(\bm{x}) \ni \bm{b}$ 是否可在多项式时间内求解？
• RQ2 当无精确解存在时，是否存在一种高效算法求解子模拉普拉斯系统的回归形式？
• RQ3 子模拉普拉斯框架如何用于将谱图理论推广至有向图和超图？
• RQ4 在假设子模变换 $F$ 为归一化的情况下，求解子模拉普拉斯系统的计算复杂度是多少？
• RQ5 能否通过在导出的流网络上使用参数化最小 $st$-割技术高效计算解？

主要发现

• 通过椭球法，子模拉普拉斯系统 $L_F(\bm{x}) \ni \bm{b}$ 可在多项式时间内求解，确立了其可 tractability。
• 即使不存在精确解，通过在拟模多面体上进行凸优化，该系统的回归形式也可在多项式时间内求解。
• 解空间由 $F$ 的核表征，其 Birkhoff 表示可多项式时间计算。
• 当给定 $\ker(F)$ 的 Birkhoff 表示时，该算法达到强多项式时间复杂度 $O(nm^2 \log(n^2/m))$。
• 该框架在单一子模拉普拉斯算子下统一处理了无向图、有向图、超图和概率模型。
• 该方法利用具有 GGT 结构的容量参数化最小 $st$-割算法，实现了所有临界 $\alpha$-值的高效计算。