[论文解读] Review of a Quantum Algorithm for Betti Numbers

Sam Gunn, Niels Kornerup|arXiv (Cornell University)|Jun 18, 2019

Topological and Geometric Data Analysis参考文献 10被引用 9

一句话总结

该论文提出了一种简化的量子算法，通过在组合拉普拉斯算子的分块对角矩阵表示上使用相位估计算法，计算单纯复形的贝蒂数。该方法相较于先前工作提高了效率，但一般情况下仍为非多项式时间复杂度，仅在特定条件下（如贝蒂数较大且谱间隙有界）才可能实现多项式时间运行。

ABSTRACT

We looked into the algorithm for calculating Betti numbers presented by Lloyd, Garnerone, and Zanardi (LGZ). We present a new algorithm in the same spirit as LGZ with the intent of clarifying quantum algorithms for computing Betti numbers. Our algorithm is simpler and slightly more efficient than that presented by LGZ. We present a thorough analysis of our algorithm, pointing out reasons that both our algorithm and that presented by LGZ do not run in polynomial time for most inputs. However, the algorithms do run in polynomial time for calculating an approximation of the Betti number to polynomial multiplicative error, when applied to some class of graphs for which the Betti number is exponentially large.

研究动机与目标

• 简化并改进Lloyd等人（2014年）最初提出的计算贝蒂数的量子算法。
• 阐明从维特里斯-里普斯滤波导出的单纯复形背景下，拓扑不变量的量子计算。
• 分析该算法的计算复杂度，并确定其可实现多项式时间运行的条件。
• 评估将该方法扩展至持久同调的可行性，后者在实际应用中更具相关性。

提出的方法

• 该算法使用格罗弗算法的一种变体，对图中所有k-单纯形进行均匀叠加态制备，其中标记项的数量未知或被估计。
• 将k-单纯形编码为n量子比特态，其汉明权重为k，从而减少量子比特需求，并通过近似计数和振幅放大实现高效的态制备。
• 该方法在分块对角厄米特矩阵B上使用相位估计算法，其平方即为每个k-单纯形维度的拉普拉斯算子Δk。
• 在k-单纯形的希尔伯特空间上应用修改后的格罗弗扩散算符，使用态制备酉算符Pk及其共轭。
• 通过量子计数估计Δk的核的维数，利用相位估计算法结果作为预言机，以计数零特征值。
• 关键组成部分是界定Δk的最小非零特征值λmin，其决定了整体门复杂度和运行时间。

实验结果

研究问题

• RQ1 与原始LGZ方案相比，能否设计出更简单、更高效的计算贝蒂数的量子算法？
• RQ2 在何种条件下，所提出的量子算法可实现计算贝蒂数的多项式时间运行？
• RQ3 是否可能将该算法扩展以计算持久贝蒂数，后者在拓扑数据分析中更具相关性？
• RQ4 拉普拉斯算子的谱间隙λmin在决定算法运行时间和可行性方面起什么作用？
• RQ5 为何原始LGZ算法和改进版本在大多数输入下均无法实现多项式时间复杂度？

主要发现

• 所提出的算法比原始LGZ算法更简单且效率稍高，其门复杂度为Õ(√(βk|Sk|)(nk√(C(n,k)/|Sk|) + n²k/λmin))。
• 由于依赖于√(βk|Sk|)和最小非零特征值λmin的倒数，该算法在大多数输入下无法实现多项式时间运行。
• 仅当βk呈指数级大且λmin远离零时，才可能实现多项式时间运行，而这些条件在实践中极为罕见。
• 原始LGZ算法所声称的复杂度基于对U_k中投影算子使用方式的错误假设，作者认为这些算子是不必要的。
• 本文表明，精确计算βk需要δ ≈ 1/βk，导致复杂度接近所提算法，而该算法仍为非多项式时间，除非βk较小且k为常数。
• 将该方法扩展至持久贝蒂数并非直接可行，因为持久同调不存在类似的拉普拉斯算子结构，且标准经典方法难以直接映射到量子设置。