[论文解读] The fast intersection transform with applications to counting paths

Andreas Björklund, Thore Husfeldt+2 more|ArXiv.org|Sep 15, 2008

Advanced Graph Theory Research参考文献 25被引用 14

一句话总结

本文提出了快速交集变换（fast intersection transform）算法，该算法在输入和输出函数支撑集的下闭包时间内，计算子集格上的线性变换。该方法可在时间 $ O^*( \exp(n cdot H(\ell/(2n)))) $ 内计数有向图中指定长度与权重的路径，突破了长路径计数的经典 $ O^*(2^n) $ 上限，在哈密顿路径与TSP的计数场景中具有应用价值。

ABSTRACT

We present an algorithm for evaluating a linear ``intersection transform'' of a function defined on the lattice of subsets of an $n$-element set. In particular, the algorithm constructs an arithmetic circuit for evaluating the transform in ``down-closure time'' relative to the support of the function and the evaluation domain. As an application, we develop an algorithm that, given as input a digraph with $n$ vertices and bounded integer weights at the edges, counts paths by weight and given length $0\leq\ell\leq n-1$ in time $O^*(\exp(n\cdot H(\ell/(2n))))$, where $H(p)=-p\log p-(1-p)\log(1-p)$, and the notation $O^*(\cdot)$ suppresses a factor polynomial in $n$.

研究动机与目标

• 开发一种高效算法，用于计算n元集的子集格上的线性交集变换。
• 将有向图中按长度与权重计数路径的时间复杂度降低至经典 $ O^*(2^n) $ 上限以下，尤其针对长路径。
• 为子集查询问题与组合优化中的精确算法提供新工具。
• 证明在聚焦于长路径时，路径计数可快于 $ O^*(2^n) $，并利用基于熵的时间界进行表征。

提出的方法

• 该算法构建一个算术电路，计算所有 $ Y \in \mathcal{G} $ 和 $ j \in \{0,\dots,n\} $ 的交集变换 $ f\iota_j(Y) = \sum_{X \subseteq U, |X \cap Y| = j} f(X) $。
• 利用输入函数支撑集的下闭包 $ \downarrow\mathcal{F} $ 和输出域的下闭包 $ \downarrow\mathcal{G} $，将计算限制在相关子集范围内。
• 核心技术将Yates算法适配于稀疏输入，实现时间复杂度 $ O^*(|\downarrow\mathcal{F}| + |\downarrow\mathcal{G}|) $。
• 在路径计数中，方法将路径在中点处分割，计算子路径的生成多项式，并利用交集变换组合，以计数恰好共享一个顶点的路径。
• 算法使用将边权重编码在变量中的生成多项式，并对大小为 $ \lfloor \ell/2 \rfloor + 1 $ 和 $ \lceil \ell/2 \rceil + 1 $ 的子集求和。
• 最终路径计数通过中间顶点 $ a $ 的求和计算，利用交集变换将子路径计数组合，其中 $ |S \cap T| = 1 $。

实验结果

研究问题

• RQ1 能否在长路径场景下打破有向图中路径计数的 $ O^*(2^n) $ 经典上界？
• RQ2 是否存在一种高效算法，用于计算子集格上的线性变换，且能利用输入与输出支撑集的稀疏性？
• RQ3 交集变换能否在 $ \ell \ll n $ 时，以亚指数时间 $ o(2^n) $ 计数给定长度与权重的路径？
• RQ4 使用如 $ H(\ell/(2n)) $ 这类熵函数是否能对稀疏情况下的路径计数复杂度提供紧致表征？

主要发现

• 快速交集变换在时间 $ O^*(|\downarrow\mathcal{F}| + |\downarrow\mathcal{G}|) $ 内计算交集变换，其中 $ \mathcal{F} $ 与 $ \mathcal{G} $ 分别为输入与输出函数的支撑集。
• 对于在加权有向图中从 $ s $ 到 $ t $ 计数长度为 $ \ell $ 的路径，该算法运行时间为 $ O^*(\exp(n \cdot H(\ell/(2n)))) $，其中 $ H $ 为二项熵函数。
• 当 $ \ell = 0.5n $ 时，运行时间为 $ O(1.7548^n) $，显著低于经典的 $ O^*(2^n) $ 上限。
• 当 $ \ell = 0.9999n $ 时，运行时间为 $ O(1.999999999^n) $，接近但未达到 $ O^*(2^n) $。
• 该算法可在相同时间界内计数指定长度的环，并通过自归约性，以相同时间复杂度列出此类路径。
• 该方法表明，使用当前工具以 $ O((2 - \epsilon)^n) $ 时间构造所有长度 $ \ell \gg n/2 $ 的路径，对任意 $ \epsilon > 0 $ 独立于 $ n $，是不可行的，因为这将对TSP问题构成突破。