[论文解读] Approximately Counting Independent Sets of a Given Size in Bounded-Degree Graphs

Ewan Davies, Will Perkins|arXiv (Cornell University)|Jan 1, 2021

Markov Chains and Monte Carlo Methods被引用 2

一句话总结

本文確立了在有界度圖中近似計數和抽樣指定大小的獨立集的計算閾值。它識別出與 ∆-正則樹上硬核模型唯一性閾值相關的臨界密度 αc(∆)，表明當密度低於 αc(∆) 時存在多項式時間近似方案，而當密度高於 αc(∆) 時則存在難度，除非 NP = RP。

ABSTRACT

We determine the computational complexity of approximately counting and sampling independent sets of a given size in bounded-degree graphs. That is, we identify a critical density α_c(Δ) and provide (i) for α < α_c(Δ) randomized polynomial-time algorithms for approximately sampling and counting independent sets of given size at most α n in n-vertex graphs of maximum degree Δ; and (ii) a proof that unless NP=RP, no such algorithms exist for α > α_c(Δ). The critical density is the occupancy fraction of hard core model on the clique K_{Δ+1} at the uniqueness threshold on the infinite Δ-regular tree, giving α_c(Δ) ~ e/(1+e)1/(Δ) as Δ → ∞.

研究动机与目标

• 確定在有界度圖中近似計數和抽樣指定大小獨立集的計算複雜度。
• 識別出一個明確的閾值 αc(∆)，用以區分此問題的可解與不可解實例。
• 將先前僅知於硬核模型分區函數 ZG(λ) 的計算相變範疇，推廣至固定大小獨立集計數問題。
• 透過系綜等價性與佔有分數分析，建立演算法複雜度與極值組合數學之間的聯繫。
• 將此框架推廣至反鐵磁性 2 自旋系統，包括伊辛模型，並證明類似的計算閾值存在。

提出的方法

• 利用統計物理中正則系綜與巨觀系綜的等價性，將硬核模型的抽樣與固定大小獨立集的均勻抽樣聯繫起來。
• 利用最大度 ∆ 的圖上 Glauber 動力學的快速混合性（當 fugacity λ < λc(∆) 時），實現硬核分佈的高效抽樣。
• 應用 Peters 和 Regts 所得的獨立多項式 ZG(λ) 的零自由區域，以確保近似過程中的集中性並避免奇點。
• 運用 Michelen 和 Sahasrabudhe 所得的硬核模型中佔用頂點數量的中心極限定理，以控制獲得指定大小獨立集的機率。
• 透過構造一個 ∆-正則且無三角形的圖 K，滿足 αK(λc(∆)) < α，將硬核分區函數近似問題歸約至固定大小獨立集計數問題，從而建立難度結果。
• 針對伊辛模型，採用類似技術：Glauber 動力學的快速混合性、複平面上的零自由區域，以及局部中心極限定理，以將 + 自旋數量與目標大小聯繫起來。

实验结果

研究问题

• RQ1 是否存在一個計算閾值，用於近似計數有界度圖中指定大小 k 的獨立集，類似於硬核模型分區函數已知的閾值？
• RQ2 對於最大度 ∆ 的 n 個頂點圖中大小 k ≤ αn 的獨立集計數問題，精確的臨界密度 αc(∆) 是什麼？
• RQ3 統計力學中巨觀系綜與正則系綜的等價性是否可用於設計固定大小獨立集計數的高效演算法？
• RQ4 固定大小獨立集計數的計算閾值是否與無限 ∆-正則樹上硬核模型的唯一性閾值一致？
• RQ5 此框架能否推廣至其他反鐵磁性 2 自旋系統（如伊辛模型），以證明類似的計算閾值存在？

主要发现

• 臨界密度 αc(∆) 等於在無限 ∆-正則樹上唯一性閾值處，硬核模型在完全圖 K∆+1 上的佔有分數，且當 ∆ → ∞ 時，αc(∆) ∼ e/(1+e)。
• 當 α < αc(∆) 時，在最大度 ∆ 的 n 個頂點圖中，對於大小 k ≤ αn 的獨立集計數，存在一個全多項式時間隨機近似方案（FPRAS）。
• 當 α > αc(∆) 時，除非 NP = RP，否則不存在多項式時間演算法能近似此類圖中大小 k ≤ αn 的獨立集數量。
• 難度結果是透過將硬核分區函數 ZG(λ) 在 λ > λc(∆) 時的近似問題歸約而建立，此問題在相同假設下已被證明是難解的。
• 該方法適用於反鐵磁性 2 自旋系統，包括伊辛模型，在此情況下存在類似的閾值 αinf(B, λc, ∆)，用以決定具有固定數量 + 自旋的狀態計數的可解性。
• 基於極值組合數學與系綜等價性論證，推測臨界閾值 αc(∆) 由完全圖 K∆+1 實現。