[论文解读] Counting all equilateral triangles in {0,1,2,...,n}^3

Eugen J. Ionaşcu|ArXiv.org|Jan 3, 2007

Mathematics and Applications参考文献 3被引用 4

一句话总结

本文提出一种计算方法，用于枚举所有顶点位于整数格点 {0,1,...,n}³ 中的等边三角形，通过从几何与数论约束推导出的丢番图方程的参数化解。作者在 Maple 中实现该算法，计算了 n=1 到 55 时的 ET(n)，并推测 limₙ→∞ ln ET(n)/ln(n+1) 存在，表明 ET(n) 的增长速率呈多项式阶 (n+1)^C，其中 C 为某常数。

ABSTRACT

We describe a procedure of counting all equilateral triangles in the three dimensional space whose coordinates are allowed only in the set $\{0,1,...,n\}$. This sequence is denoted here by ET(n) and it has the entry A102698 in "The On-Line Encyclopedia of Integer Sequences". The procedure is implemented in Maple and its main idea is based on the results in \cite{eji}. Using this we calculated the values ET(n) for n=1..55 which are included here. Some facts and conjectures about this sequence are stated. The main of them is that $\ds \lim_{n o \infty} \frac{\ln ET(n)}{\ln n+1}$ exists.

研究动机与目标

• 系统计算所有顶点受限于整数格点 {0,1,...,n}³ 的等边三角形数量。
• 基于 ℤ³ 中等边三角形的数论参数化，开发一种计算程序。
• 计算 n=1 到 55 时的 ET(n)，并分析该序列的增长行为。
• 研究 ET(n) 的渐近增长速率，并推测 limₙ→∞ ln ET(n)/ln(n+1) = C 是否存在。

提出的方法

• 利用定理 1.1 将 ℤ³ 中等边三角形的边长表征为 √[2(m²−mn+n²)]，其中 m,n 为整数。
• 应用定理 1.2，生成所有满足 a²+b²+c²=3d² 的整数解 (a,b,c,d)，这些解定义了包含此类三角形的平面的法向量。
• 利用定理 1.3，通过整数参数 m,n 及方程 2(a²+b²)=s²+3r² 的解 (r,s)，对给定平面内的所有三角形进行参数化。
• 在 Maple 中实现递归过程，计算在立方体 [0,n]³ 内最小三角形的所有有效平移与旋转数量。
• 利用几何不变量如 transl(T)、inters(T) 和 intersch(T)，通过确保放置位置互不重叠，避免重复计数。
• 将所有有效 (a,b,c,d) 与 (r,s) 对的贡献相加，以 n 的多项式项之和形式计算 ET(n)。

实验结果

研究问题

• RQ1 对于给定的 n，{0,1,...,n}³ 中顶点位于三维整数格点的等边三角形的确切数量是多少？
• RQ2 如何利用数论参数化系统地生成所有此类三角形？
• RQ3 当 n→∞ 时，ET(n)（此类三角形的数量）的渐近增长速率如何？
• RQ4 极限 limₙ→∞ ln ET(n)/ln(n+1) 是否存在？若存在，其值是多少？

主要发现

• 作者计算了 n=1 到 55 时的 ET(n)，其值从 n=1 时的 8 增加到 n=55 时的 488,902,856。
• 序列 ET(n) 的增长速度快于任意三次多项式，但似乎被 (n+1)^5 所限制，表明对所有 n 有 ET(n) ≤ (n+1)^5。
• 作者观察到序列 aₙ = ln ET(n)/ln(n+1) 是递增的，且在数值上被 5 上限所限制，支持了该极限存在的猜想。
• ET(n) 的增长被建模为 ζ = n−1 的多项式贡献之和，每个新 n 值增加形如 uₙ(ζ−n+1)^3 + vₙ(ζ−n+1)^2 + wₙ(ζ−n+1) + sₙ 的项。
• 发现了一个反例，推翻了参数化中的简化假设，但该方法对所有满足 d ≤ 4096 且 gcd(a,b,c)=1 的奇数 d 仍保持有效。
• 作者推测 limₙ→∞ ln ET(n)/ln(n+1) = C 存在，表明 ET(n) 呈精确的多项式增长速率，阶为 (n+1)^C，其中 C 为某常数。