[논문 리뷰] Computing the Independence Polynomial: from the Tree Threshold down to the Roots

Nicholas J. A. Harvey, Piyush Srivastava외 1명|arXiv (Cornell University)|2016. 08. 07.

Markov Chains and Monte Carlo Methods참고 문헌 44인용 수 4

한 줄 요약

이 논문은 원점 중심의 루트 없는 복소다면체 내에서 복소수 및 음수 정점 활성도를 갖는 다변수 독립 다항식 $ Z(\mathbf{z}) $ 를 결정론적으로 근사하는 알고리즘을 제시한다. 비균일한 복소수 매개변수를 사용하는 새로운 다변수 상관관계 감쇠 기법을 도입하여, $(1+\epsilon)$-근사가 시간 $(n/\epsilon\alpha)^{O(\log d / \sqrt{\alpha})}$ 내에 달성되며, 기존의 계산 가능한 영역을 통합하고 지수 항에서 $1/\sqrt{\alpha}$ 의 의존성의 날카로움을 증명한다.

ABSTRACT

We study an algorithm for approximating the multivariate independence polynomial $Z(\mathbf{z})$, with negative and complex arguments, an object that has strong connections to combinatorics and to statistical physics. In particular, the independence polynomial with negative arguments, $Z(-\mathbf{p})$, determines the Shearer region, the maximal region of probabilities to which the Lovasz Local Lemma (LLL) can be extended (Shearer 1985). In statistical physics, complex zeros of the independence polynomial relate to existence of phase transitions. Our main result is a deterministic algorithm to compute approximately the independence polynomial in any root-free complex polydisc centered at the origin. Our algorithm is essentially the same as Weitz's algorithm for positive parameters up to the tree uniqueness threshold, and the core of our analysis is a novel multivariate form of the correlation decay technique, which can handle non-uniform complex parameters. In particular, in the univariate real setting our work implies that Weitz's algorithm works in an interval between two critical points $(λ'_c(d), λ_c(d))$, and outside of this interval an approximation of $Z(\mathbf{z})$ is known to be NP-hard. As an application, we give a sub-exponential time algorithm for testing approximate membership in the Shearer region. We also give a new rounding based deterministic algorithm for Shearer's lemma (an extension of the LLL), which, however, runs in sub-exponential time. On the hardness side, we prove that evaluating $Z(\mathbf{z})$ at an arbitrary point in Shearer's region, and testing membership in Shearer's region, are #P-hard problems. We also establish the best possible dependence of the exponent of the run time of Weitz's correlation decay technique in the negative regime on the distance to the boundary of the Shearer region.

연구 동기 및 목표

• 복소수 및 음수 정점 활성도를 포함한, 계산 가능한 영역 전역에서 독립 다항식 $ Z(\mathbf{z}) $ 를 근사하는 통합 알고리즘을 개발하는 것.
• 웨이츠의 상관관계 감쇠 프레임워크를 비균일한 복소수 매개변수를 다룰 수 있도록 확장하여, 양의 실수 영역을 넘어서는 근사가 가능하도록 하는 것.
• 특히 임계 임계점 근처에서 계산 난이도가 발생하는 경우에 대해, $ Z(\mathbf{z}) $ 근사의 날카로운 실행 시간 범위를 확립하는 것.
• 셰이어의 영역에 속하는지 테스트하는 데 효율적인 알고리즘과 슬랙 $ \alpha $ 가 있는 알고리즘적 라보스 로컬 렘마를 제공하는 것.
• 셰이어의 영역 내 임의의 점에서 $ Z(\mathbf{z}) $ 를 근사하거나, 셰이어의 영역에 속하는지 테스트하는 것이 #P-난이도 문제라는 것을 증명하는 것.

제안 방법

• 핵심 방법은 비균일한 복소수 정점 활성도를 다룰 수 있는 새로운 다변수 상관관계 감쇠 형태를 사용하며, 웨이츠의 접근법을 양의 실수를 넘어서 확장한다.
• 알고리즘은 루트 없는 복소다면체 내에서 작동하여, 모든 $ |z'_i| \leq (1+\alpha)|z_i| $ 에 대해 $ Z(\mathbf{z}') \neq 0 $ 임을 보장함으로써 해석성과 안정성을 확보한다.
• 핵심 구성 요소 중 하나는 오차 민감도 매개변수로, 계산 트리에서 상관관계 감쇠의 정도를 제한하여 정확한 값으로의 수렴을 보장한다.
• 분석은 복소해석학 부등식과 재귀식 $ f(x) = 1 - (1 - \lambda x)^d $ 에 대한 재귀적 경계를 사용하며, 고정점으로의 수렴을 추적한다.
• 실행 시간은 $ \epsilon $-정확도를 확보하기 위해 필요한 계산 트리의 레벨 수를 제한함으로써 유도되며, 이는 $ \alpha $ 와 $ d $ 에 따라 달라진다.
• 난이도 결과는 기존의 #P-난이도 문제로의 감소를 통해 확립되며, 셰이어의 영역 내에서도 근사가 어렵다는 것을 보여준다.

실험 결과

연구 질문

• RQ1 웨이츠의 상관관계 감쇠 방법은 복소수 및 음수 정점 활성도를 갖는 독립 다항식에 대해 확장 가능한가?
• RQ2 독립 다항식의 임계점 근처에서 실행 시간의 최적 의존성은 슬랙 매개변수 $ \alpha $ 에 대해 어떻게 되는가?
• RQ3 $ Z(\mathbf{z}) $ 가 근사적으로 계산 가능한 영역이 루트 없는 영역(즉, 0이 아닌 영역)과 일치하는가?
• RQ4 다항식 평가와 루트 없는 영역을 활용하여 알고리즘적 라보스 로컬 렘마를 효율적으로 구현할 수 있는가?
• RQ5 셰이어의 영역 내에서의 소속성 테스트 또는 셰이어의 영역 내에서 $ Z(\mathbf{z}) $ 를 근사하는 것은 계산적으로 어려운가?

주요 결과

• 알고리즘은 원점 중심의 루트 없는 다면체 내에서 임의의 복소수 $ \mathbf{z} $ 에 대해 $(1+\epsilon)$-근사가 시간 $ (n/\epsilon\alpha)^{O(\log d / \sqrt{\alpha})} $ 내에 달성된다.
• 지수 항에서의 $ 1/\sqrt{\alpha} $ 의 의존성은, 임계점 근처에서 재귀식 $ f^l(0) $ 의 수렴 속도에 대한 하한을 통해 최적임을 입증한다.
• 일변수 실수 $ \lambda $ 에 대해 알고리즘은 구간 $ (-\lambda_c'(d), \lambda_c(d)) $ 내에서 작동하며, 이는 근사가 NP-난이도로 변하는 계산 임계점과 정확히 일치한다.
• 소속성 테스트를 다중적 오차 $ 1+\alpha $ 내에서 수행하는 알고리즘은 시간 $ (n/\alpha)^{O(\sqrt{n/\alpha} \log d)} $ 내에 작동한다.
• 이를 바탕으로 $ n $ 개의 사건과 슬랙 $ \alpha $ 를 갖는 셰이어의 렘마에 대한 결정론적 알고리즘은 시간 $ (nm/\alpha)^{O(\sqrt{m/\alpha} \log d)} $ 내에 작동한다.
• 셰이어의 영역 내에서 $ Z(\mathbf{z}) $ 를 평가하거나 셰이어의 영역에 속하는지 테스트하는 것은 모두 #P-난이도 문제이다.