[论文解读] Enumeration of integral tetrahedra

Sascha Kurz|ArXiv.org|Apr 8, 2008

Analytic Number Theory Research参考文献 14被引用 5

一句话总结

本文提出了一种高效算法，用于枚举直径至多为 d 的非同构整数四面体，时间复杂度为 O(d⁵)。推导出满足三角不等式的对称 4×4 距离矩阵数量的精确公式，并计算了 d ≤ 1000 时的 α(3,d)，确定几何概率 P₃ 落在 0.090 到 0.111 之间。

ABSTRACT

We determine the numbers of integral tetrahedra with diameter $d$ up to isomorphism for all $d\le 1000$ via computer enumeration. Therefore we give an algorithm that enumerates the integral tetrahedra with diameter at most $d$ in $O(d^5)$ time and an algorithm that can check the canonicity of a given integral tetrahedron with at most 6 integer comparisons. For the number of isomorphism classes of integral $4 imes 4$ matrices with diameter $d$ fulfilling the triangle inequalities we derive an exact formula.

研究动机与目标

• 确定 d ≤ 1000 时非同构整数四面体的数量，其直径至多为 d。
• 开发一种在 O(d⁵) 时间内枚举整数四面体的高效算法。
• 推导出 entries 属于 {1,…,d} 的对称 4×4 距离矩阵数量的精确闭式表达式，且满足三角不等式。
• 计算几何概率 P₃，即随机 4×4 距离矩阵（entries 属于 [0,1]）对应欧几里得四面体的概率。
• 通过最多六次整数比较，实现快速的规范性检查，用于无同构生成。

提出的方法

• 本文使用有序生成法（通过隐式枚举）系统地生成所有至多同构的整数四面体。
• 应用加权形式的 Cauchy-Frobenius 引理由对称群 S₄ 作用下的无标签矩阵进行计数。
• 通过决策树算法强制执行规范标记条件，确保在最多六次整数比较内完成同构性检查。
• 算法通过排序边长以字典序生成四面体，确保不产生同构重复。
• 利用由对称性和组合枚举推导出的闭式表达式，计算有效距离矩阵的数量。
• 几何概率 P₃ 基于离散逼近，估计为极限 lim_{d→∞} 4·α(3,d)/d⁵。

实验结果

研究问题

• RQ1 对于 d ≤ 1000，非同构整数四面体（直径至多为 d）的确切数量是多少？
• RQ2 能否设计一种高效算法，在 O(d⁵) 时间内枚举所有整数四面体，同时避免同构重复？
• RQ3 对于 entries 属于 {1,…,d} 且满足三角不等式、直径至多为 d 的对称 4×4 距离矩阵，其确切数量是多少？
• RQ4 随机 (4×4) 距离矩阵（entries 属于 [0,1]）对应欧几里得四面体的渐近几何概率 P₃ 是多少？
• RQ5 能否通过最多六次整数比较，实现整数四面体的规范性检查？

主要发现

• 使用 O(d⁵) 算法，对所有 d ≤ 1000 精确计算出直径至多为 d 的非同构整数四面体数量。
• 推导出 entries 属于 {1,…,d} 的对称 4×4 距离矩阵数量的精确公式，且满足三角不等式： $$ \hat{\alpha}_{\leq}(d,3) = \begin{cases} \frac{17d^6 + 425d^4 + 1628d^2}{2880} & \text{若 } d \equiv 0 \pmod{2}, \\ \frac{17d^6 + 425d^4 + 1763d^2 + 675}{2880} & \text{若 } d \equiv 1 \pmod{2}. \end{cases} $$
• 直径恰好为 d 的此类矩阵数量由下式给出： $$ \hat{\alpha}(d,3) = \begin{cases} \frac{34d^5 - 85d^4 + 680d^3 - 962d^2 + 1776d - 960}{960} & \text{若 } d \equiv 0 \pmod{2}, \\ \frac{34d^5 - 85d^4 + 680d^3 - 908d^2 + 1722d - 483}{960} & \text{若 } d \equiv 1 \pmod{2}. \end{cases} $$
• 随机 4×4 距离矩阵（entries 属于 [0,1]）对应欧几里得四面体的几何概率 P₃ 满足 0.090 ≤ P₃ ≤ 0.111。
• 规范性检查算法最多需要六次整数比较，且已针对有序生成过程进行优化。
• 当 d = 800 时，计算 α(800,3) 仅耗时 3.3 小时，表明该算法相比朴素循环具有显著效率优势。