[论文解读] Computing Generalized Rank Invariant for 2-Parameter Persistence Modules via Zigzag Persistence and Its Applications

Tamal K. Dey, Woojin Kim+1 more|arXiv (Cornell University)|Nov 30, 2021

Topological and Geometric Data Analysis被引用 2

一句话总结

本文提出了一种新颖的方法，通过将问题约化为在 Z² 中有限区间边界上的 zigzag 持久性，来计算二维参数持久性模的广义秩不变量。关键贡献是一条定理，表明区间 I 上的广义秩等于限制在 ∂I 上的 zigzag 模的条形码重数，从而实现了 O(t^ω) 时间内的高效计算，其中 t 是单形的数量，ω ∈ [2, 2.373) 是矩阵乘法指数。这导致了测试区间可分解性以及计算可分解模条形码的改进算法。

ABSTRACT

The notion of generalized rank invariant in the context of multiparameter persistence has become an important ingredient for defining interesting homological structures such as generalized persistence diagrams. Naturally, computing these rank invariants efficiently is a prelude to computing any of these derived structures efficiently. We show that the generalized rank over a finite interval $I$ of a $\mathbb{Z}^2$-indexed persistence module $M$ is equal to the generalized rank of the zigzag module that is induced on a certain path in $I$ tracing mostly its boundary. Hence, we can compute the generalized rank over $I$ by computing the barcode of the zigzag module obtained by restricting the bifiltration inducing $M$ to that path. If the bifiltration and $I$ have at most $t$ simplices and points respectively, this computation takes $O(t^ω)$ time where $ω\in[2,2.373)$ is the exponent of matrix multiplication. Among others, we apply this result to obtain an improved algorithm for the following problem. Given a bifiltration inducing a module $M$, determine whether $M$ is interval decomposable and, if so, compute all intervals supporting its summands.

研究动机与目标

• 开发一种高效算法，用于计算 Z² 索引持久性模的广义秩不变量。
• 解决确定二维参数持久性模是否为区间可分解的计算挑战，若可分解，则计算其区间直和分量。
• 通过利用 Z² 中区间边界的 zigzag 持久性，降低广义秩计算的复杂度。
• 改进现有算法，用于二维参数持久性模的区间可分解性测试与条形码计算。

提出的方法

• 本文建立了理论等价性：Z² 模 M 在有限区间 I 上的广义秩等于在 I 的边界 ∂I 上诱导的 zigzag 模的广义秩。
• 它沿着 ∂I 构建了一条 zigzag 路径，以追踪 Z² 中区间 I 的边界，同时保持与秩计算相关的偏序结构。
• 该方法将二维参数持久性问题约化为计算 zigzag 模的条形码，该问题可使用已知算法高效求解。
• 该算法利用了从限制在 ∂I 上的 zigzag 模的条形码中恢复 I 上秩函数的事实，通过极限到余极限映射和 Möbius 反演实现。
• 计算复杂度被限制在 O(t^ω) 时间内，其中 t 是双过滤中单形的数量，ω 是矩阵乘法指数。
• 该方法使得设计一种新算法 INTERVAL 成为可能，该算法可在 O(t^{ω+2}) 时间内计算任意有限区间可分解 Z² 模的条形码。

实验结果

研究问题

• RQ1 能否通过将问题约化为 zigzag 持久性问题，高效计算二维参数持久性模的广义秩不变量？
• RQ2 能否仅通过区间的边界 ∂I 来刻画有限区间 I ⊆ Z² 上的广义秩？
• RQ3 能否通过在区间边界上使用 zigzag 持久性，更高效地测试 Z² 模的区间可分解性？
• RQ4 使用该新方法计算区间可分解 Z² 模的条形码的计算复杂度是多少？
• RQ5 该方法能否推广到更高参数持久性模（d > 2）？
• RQ6 能否通过用高效的 zigzag 计算替代指数级枚举，改进现有区间可分解性算法的复杂度？

主要发现

• Z² 模 M 在有限区间 I 上的广义秩等于限制在 I 的边界 ∂I 上的 zigzag 模的条形码中完整条形的重数。
• 广义秩不变量的计算被约化为计算 zigzag 模的条形码，该计算可在 O(t^ω) 时间内完成，其中 t 是单形数量，ω ∈ [2, 2.373) 是矩阵乘法指数。
• 算法 INTERVAL 可在 O(t^{ω+2}) 时间内计算任意有限区间可分解 Z² 模的条形码，优于先前的方法。
• 算法 ISINTERVALDECOMP 通过利用基于边界的 zigzag 约化，可在 O(t^{ω+2}) 时间内测试 Z² 模的区间可分解性。
• 该方法避免了先前方法（如 Asashiba 等人的算法）所需的指数级区间枚举。
• 理论基础由一个证明支持：通过在 ∂I 上使用 zigzag 持久性计算的广义秩，与在 I 上的广义秩一致，该证明利用了极限到余极限映射和 Möbius 反演。