[论文解读] Factoring Permutation Matrices Into a Product of Tridiagonal Matrices

Michael Daniel Samson, Martianus Frederic Ezerman|arXiv (Cornell University)|Jul 20, 2010

graph theory and CDMA systems参考文献 4被引用 3

一句话总结

本文证明了吉尔伯特·斯特朗格的猜想：任何带宽为 $ w $ 的置换矩阵均可通过贪心并行冒泡排序算法，分解为最多 $ 2w - 1 $ 个带宽为 1 的置换矩阵。该方法通过相邻行交换系统性地减少逆序对，确保在少于 $ 2w $ 步内收敛，并在 $ n = 9 $ 以内通过理论与实证验证，同时引入了贪心冒泡矩阵以优化分解过程。

ABSTRACT

Gilbert Strang posited that a permutation matrix of bandwidth $w$ can be written as a product of $N < 2w$ permutation matrices of bandwidth 1. A proof employing a greedy ``parallel bubblesort'' algorithm on the rows of the permutation matrix is detailed and further points of interest are elaborated.

研究动机与目标

• 证明吉尔伯特·斯特朗格的猜想：带宽为 $ w $ 的置换矩阵可表示为最多 $ 2w - 1 $ 个带宽为 1 的置换矩阵的乘积。
• 开发并形式化一种贪心的“并行冒泡排序”算法，通过迭代方式减少置换矩阵中的逆序对。
• 引入并分析贪心冒泡矩阵（GBM）的概念，作为构建最优或近似最优分解的工具。
• 建立所需因子数的理论界，特别是证明 $ \mathrm{fact}(P) < 2w $，并在 $ n \leq 9 $ 的范围内通过实证加以验证。

提出的方法

• 本文采用贪心并行冒泡排序策略，通过仅在形成逆序对时交换相邻行，识别并解决置换矩阵中的逆序对。
• 定义了“减少”与“超越”交换，用于逆序块内部，以确保每一步均减少总逆序对数，且不会重新引入先前的逆序对。
• 引入“贪心冒泡矩阵” $ G $ 的概念，使得 $ GP $ 与 $ P $ 不共享任何共同的逆序对，从而确保每次迭代均实现最大进度。
• 在对称群的凯莱图上使用广度优先生成树算法，计算最小分解长度 $ \overline{\mathrm{fact}}(P) $，其中连接集由带宽为 1 的置换定义。
• 该算法被形式化为机会性超越矩阵算法，通过结构化、迭代的方式处理逆序块，并应用交换以保持进度。
• 理论分析利用“分块”（块对角分量）与行/列符号模式的概念，对因子化过程中的置换矩阵进行分类与管理。

实验结果

研究问题

• RQ1 每个带宽为 $ w $ 的置换矩阵是否都能如斯特朗格所猜想的那样，分解为少于 $ 2w $ 个带宽为 1 的置换矩阵？
• RQ2 应用于置换矩阵行的贪心并行冒泡排序算法，是否总能在少于 $ 2w $ 步内收敛到单位矩阵？
• RQ3 贪心冒泡矩阵在确保最小或近似最小分解长度方面起什么作用？
• RQ4 置换矩阵的结构性质——如分块、逆序块以及行/列符号——如何影响分解过程？
• RQ5 是否存在一种系统性的算法方法，用于计算或近似 $ \overline{\mathrm{fact}}(P) $，即给定 $ P $ 所需的最少带宽为 1 的因子数？

主要发现

• 本文证明了斯特朗格的猜想：任意 $ n \times n $ 带宽为 $ w $ 的置换矩阵可分解为最多 $ 2w - 1 $ 个带宽为 1 的置换矩阵。
• 贪心并行冒泡排序算法确保在最多 $ 2w - 1 $ 步内收敛到单位矩阵，因为每次迭代严格减少逆序对数。
• 对所有 $ n \leq 9 $ 的 $ n \times n $ 置换矩阵进行实证测试，结果表明 $ \mathrm{fact}(P) \leq \overline{\mathrm{fact}}(P) + \lfloor n/3 \rfloor $，支持了贪心分解的可行性。
• 带宽 $ w \leq 1 $ 的 $ n \times n $ 置换矩阵的数量为第 $ n $ 个斐波那契数 $ F_n $，为分解复杂性提供了组合基准。
• 贪心冒泡矩阵 $ G $ 的概念确保 $ GP $ 与 $ P $ 不共享任何共同的逆序对，从而实现高效且无冗余的分解步骤。
• 实证数据支持猜想 8：所有测试的带宽为 $ w $ 的置换矩阵均可分解为少于 $ 2w $ 个贪心冒泡矩阵，表明一般界 $ \mathrm{fact}(P) < 2w $ 成立。