[論文レビュー] Simulated Annealing is a Polynomial-Time Approximation Scheme for the Minimum Spanning Tree Problem

Benjamin Doerr, Amirhossein Rajabi+1 more|arXiv (Cornell University)|Apr 5, 2022

VLSI and FPGA Design Techniques参考文献 20被引用数 4

ひとこと要約

本稿は、$1 - 1/\ell$ の乗法的冷却スケジュールを適切に選ぶことで、$\ell = (mn \ln m)^{1+1/\varepsilon + o(1)}$ の要因を持つ焼きなまし法が、高確率で多項式時間内に $(1+\varepsilon)$-近似最小全域木を計算できることを証明している。主な結果は、Wegenerが2005年に提起した予想である、焼きなまし法がMST問題に対する多項式時間近似スキーム（PTAS）として機能することを裏付けている。

ABSTRACT

We prove that Simulated Annealing with an appropriate cooling schedule computes arbitrarily tight constant-factor approximations to the minimum spanning tree problem in polynomial time. This result was conjectured by Wegener (2005). More precisely, denoting by $n, m, w_{\max}$, and $w_{\min}$ the number of vertices and edges as well as the maximum and minimum edge weight of the MST instance, we prove that simulated annealing with initial temperature $T_0 \ge w_{\max}$ and multiplicative cooling schedule with factor $1-1/\ell$, where $\ell = \omega (mn\ln(m))$, with probability at least $1-1/m$ computes in time $O(\ell (\ln\ln (\ell) + \ln(T_0/w_{\min}) ))$ a spanning tree with weight at most $1+\kappa$ times the optimum weight, where $1+\kappa = \frac{(1+o(1))\ln(\ell m)}{\ln(\ell) -\ln (mn\ln (m))}$. Consequently, for any $\epsilon>0$, we can choose $\ell$ in such a way that a $(1+\epsilon)$-approximation is found in time $O((mn\ln(n))^{1+1/\epsilon+o(1)}(\ln\ln n + \ln(T_0/w_{\min})))$ with probability at least $1-1/m$. In the special case of so-called $(1+\epsilon)$-separated weights, this algorithm computes an optimal solution (again in time $O( (mn\ln(n))^{1+1/\epsilon+o(1)}(\ln\ln n + \ln(T_0/w_{\min})))$), which is a significant speed-up over Wegener's runtime guarantee of $O(m^{8 + 8/\epsilon})$.

研究の動機と目的

• Wegenerの2005年の予想を解決する。すなわち、焼きなまし法が$(1+\varepsilon)$-近似最小全域木を多項式時間で計算できることを示すこと。
• 現代の解析的手法を用いて、最小全域木問題に対する焼きなまし法の厳密な実行時間解析を提供すること。
• $(1+\varepsilon)$-分離された重みに対して、Wegenerの以前の実行時間上限$O(m^{8+8/\varepsilon})$を改善すること。
• 任意の辺の重みに対して、焼きなまし法を用いた一般の多項式時間近似スキーム（PTAS）を確立すること。
• 温度の低下と初期温度の境界に基づいた実用的な停止基準を提示すること。

提案手法

• 冷却要因$\beta = 1 - 1/\ell$の乗法的冷却スケジュールを使用し、$\ell = (mn \ln m)^{1+1/\varepsilon + o(1)}$ とする。
• メトロポリス基準に基づく温度依存の受理確率を採用する。
• ドリフト解析と集中不等式を用いて、より良いまたは近似的に最適な全域木への遷移確率を分析する。
• 時間経過に伴い、アルゴリズムの状態を最適または近似的に最適な解と比較するためのカップリング論法を適用する。
• 対数的温度低下モデルを用いて、収束に至るまでの反復回数を上限付ける。
• $(1+\varepsilon)$-分離性を活用し、臨界閾値の間の途中の辺の重みが存在しないようにすることで、正確な収束を可能にする。

実験結果

リサーチクエスチョン

• RQ1 乗法的冷却スケジュールを用いた焼きなまし法が、最小全域木問題に対して$(1+\varepsilon)$-近似を達成できるか。
• RQ2 一般の辺の重みに対して、焼きなまし法の実行時間は$n$、$m$、および$1/\varepsilon$に関して多項式的にスケーリングするか。
• RQ3 冷却要因$\beta = 1 - 1/\ell$が収束時間と近似品質に与える影響は何か。
• RQ4 辺の重みが$(1+\varepsilon)$-分離されている場合、最適な最小全域木を確実に見つけられるか。その場合の実行時間は何か。
• RQ5 焼きなまし法が$(1+\varepsilon)$-近似を達成するための最もタイトな実行時間上限は何か。また、それ以前の研究と比較するとどうか。

主な発見

• 冷却要因$\beta = 1 - 1/\ell$、$\ell = (mn \ln m)^{1+1/\varepsilon + o(1)}$ を用いた焼きなまし法は、確率$1 - 1/m$以上で、時間$O((mn \ln n)^{1+1/\varepsilon + o(1)}(\ln \ln n + \ln(T_0/w_{\min})))$ で$(1+\varepsilon)$-近似最小全域木を計算する。
• $(1+\varepsilon)$-分離された重みに対しては、同じ時間上限内で最適な最小全域木が得られ、Wegenerの$O(m^{8+8/\varepsilon})$の上限を顕著に改善する。
• 実行時間は$n$、$m$、および$1/\varepsilon$に関して多項式的であるため、焼きなまし法がMST問題に対する多項式時間近似スキーム（PTAS）であることが確認された。
• 初期温度$T_0$は、重い辺が高確率で含められるように、$w_{\max}$以上でなければならない。
• アルゴリズムの実行時間は$O(\ell \ln(T_0/w_{\min}))$で上限づけられており、$w_{\min}$まで温度が低下するため、これはタイトである。
• 解析により、温度が$w_{\min}$未満になると、新たな辺が追加できなくなることが明らかになった。これにより、自然な停止基準が得られる。