[论文解读] Computing the Independence Polynomial: from the Tree Threshold down to the Roots

Nicholas J. A. Harvey, Piyush Srivastava+1 more|arXiv (Cornell University)|Aug 7, 2016

Markov Chains and Monte Carlo Methods参考文献 44被引用 4

一句话总结

本文提出了一种确定性算法，用于在原点为中心的无根复数多圆柱内近似计算任意顶点活动为复数或负数时的多元独立多项式 $ Z(\mathbf{z}) $。通过使用一种新颖的多变量相关衰减技术，结合非均匀复参数，该算法在时间 $ (n/\epsilon\alpha)^{O(\log d / \sqrt{\alpha})} $ 内实现了 $ (1+\epsilon) $-近似，统一了已知的可计算区域，并证明了指数中 $ 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}) $。
• 将 Weitz 的相关衰减框架扩展至处理非均匀复参数，从而实现对正实数范围之外的近似。
• 建立近似 $ Z(\mathbf{z}) $ 的紧致时间复杂度边界，特别是在计算困难性出现的临界阈值附近。
• 为测试 Shearer 区域的成员资格以及具有松弛参数 $ \alpha $ 的算法化 Lovász 局部引理提供高效算法。
• 证明在 Shearer 区域内任意点近似 $ Z(\mathbf{z}) $ 或测试其成员资格均为 #P-难问题。

提出的方法

• 核心方法采用一种新颖的多变量相关衰减形式，可处理非均匀复数顶点活动，将 Weitz 的方法从正实数范围扩展至更广领域。
• 该算法运行于一个无根复数多圆柱内，确保对所有 $ |z'_i| \leq (1+\alpha)|z_i| $，有 $ Z(\mathbf{z}') \neq 0 $，从而保证了解析性与稳定性。
• 关键组成部分是误差敏感性参数，用于限制计算树中相关性的衰减，确保收敛至正确值。
• 分析过程使用复分析不等式和递推关系 $ f(x) = 1 - (1 - \lambda x)^d $ 的递归边界，追踪其收敛至不动点的过程。
• 时间复杂度通过限制为达到 $ \epsilon $-精度所需的计算树层级数而推导得出，其依赖于 $ \alpha $ 和 $ d $。
• 通过从已知的 #P-难问题进行归约，建立了下界结果，表明即使在 Shearer 区域内，近似也是计算困难的。

实验结果

研究问题

• RQ1 Weitz 的相关衰减方法能否被扩展以处理独立多项式中复数和负数顶点活动？
• RQ2 在独立多项式临界阈值附近，运行时间对松弛参数 $ \alpha $ 的最优依赖关系是什么？
• RQ3 $ Z(\mathbf{z}) $ 可近似计算的区域是否与非零区域（即无根区域）一致？
• RQ4 能否通过多项式求值和无根区域，高效实现算法化 Lovász 局部引理？
• RQ5 测试 Shearer 区域的成员资格或在 Shearer 区域内近似 $ Z(\mathbf{z}) $ 是否具有计算困难性？

主要发现

• 该算法在时间 $ (n/\epsilon\alpha)^{O(\log d / \sqrt{\alpha})} $ 内实现了对任意复数 $ \mathbf{z} $ 在无根多圆柱内的 $ (1+\epsilon) $-近似。
• 指数中 $ 1/\sqrt{\alpha} $ 的依赖关系是紧致的，这通过在临界点附近递推 $ f^l(0) $ 的收敛速率下界得到证明。
• 对于单变量实数 $ \lambda $，该算法在区间 $ (-\lambda_c'(d), \lambda_c(d)) $ 内有效，与已知的计算阈值一致，即在此区间外近似将变为 NP-难。
• 一个在乘法误差 $ 1+\alpha $ 内测试 Shearer 区域成员资格的算法，其运行时间为 $ (n/\alpha)^{O(\sqrt{n/\alpha} \log d)} $。
• 一个针对 $ n $ 个事件和松弛参数 $ \alpha $ 的 Shearer 引理的确定性算法，其运行时间为 $ (nm/\alpha)^{O(\sqrt{m/\alpha} \log d)} $。
• 在 Shearer 区域内任意点求值 $ Z(\mathbf{z}) $ 和测试其成员资格均为 #P-难问题。