[論文レビュー] 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⁵)時間で整数四面体を列挙する効率的なアルゴリズムを開発すること。
• {1,…,d}の要素をもつ4×4対称距離行列の数に対する正確な閉形式の公式を導出すること。ただし、三角不等式を満たすものに限る。
• [0,1]の要素をもつランダムな4×4距離行列がユークリッド四面体に対応する幾何学的確率P₃を計算すること。
• 同型を避ける生成のための高速な正規化チェックを、最大6つの整数比較で実装すること。

提案手法

• 論文は、暗黙の列挙を用いた順序付き生成を用いて、同型を除いたすべての整数四面体を体系的に生成する。
• 対称群S₄の作用の下での非ラベル付き行列の数を数えるために、重み付き形のコーシー＝フロベニウスの補題を適用する。
• 同型性のチェックを、最大6つの整数比較で行うことができる決定木アルゴリズムを用いて、正規化ラベル付けの条件を強制する。
• 辺の長さをソートして辞書式順序で四面体を生成することで、同型の重複が生じないことを保証する。
• 対称性と組合せ的列挙に基づいて導出された閉形式式を用いて、有効な距離行列の数を計算する。
• 幾何学的確率P₃は、離散近似に基づき、lim_{d→∞} 4·α(3,d)/d⁵の極限として推定される。

実験結果

リサーチクエスチョン

• RQ1 d ≤ 1000の範囲で、直径がd未満の非同型の整数四面体の正確な数は何か？
• RQ2 同型の重複を避けるO(d⁵)時間ですべての整数四面体を列挙できる効率的なアルゴリズムは存在するか？
• RQ3 {1,…,d}の要素をもつ4×4対称距離行列のうち、三角不等式を満たし、直径がd未満であるものの正確な数は何か？
• RQ4 [0,1]の要素をもつランダムな(4×4)距離行列がユークリッド四面体に対応する漸近的幾何学的確率P₃は何か？
• RQ5 整数四面体の正規化チェックを、最大6つの整数比較で実装できるか？

主な発見

• 直径がd未満の非同型の整数四面体の数は、すべてのd ≤ 1000についてO(d⁵)アルゴリズムを用いて正確に計算された。
• {1,…,d}の要素をもつ4×4対称距離行列の数で、三角不等式を満たすものの正確な公式は以下の通りである： $$ \hat{\alpha}_{\leq}(d,3) = \begin{cases} \frac{17d^6 + 425d^4 + 1628d^2}{2880} & \text{if } d \equiv 0 \pmod{2}, \\ \frac{17d^6 + 425d^4 + 1763d^2 + 675}{2880} & \text{if } 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{if } d \equiv 0 \pmod{2}, \\ \frac{34d^5 - 85d^4 + 680d^3 - 908d^2 + 1722d - 483}{960} & \text{if } d \equiv 1 \pmod{2}. \end{cases} $$
• [0,1]の要素をもつランダムな4×4距離行列がユークリッド四面体に対応する幾何学的確率P₃は、0.090 ≤ P₃ ≤ 0.111を満たす。
• 正規化チェックアルゴリズムは最大6つの整数比較で実行可能であり、順序付き生成プロセスでの使用に最適化されている。
• d = 800の場合、α(800,3)の計算にたった3.3時間で完了し、単純なループと比較して本アルゴリズムの効率性が示された。