Carl Størmer and his Numbers

基本信息

• 论文ID: 2511.03030
• 标题: Carl Størmer and his Numbers
• 作者: Matthew Kroesche, Lance L. Littlejohn, Graeme Reinhart
• 分类: math.HO (数学史) math.NT (数论)
• 发表时间: 2025年11月4日
• 论文链接: https://arxiv.org/abs/2511.03030

摘要

在费马二平方定理的许多证明中，二次同余方程 $x^2 \equiv -1 \bmod p$ 的最小正剩余解 $x_0$ 起着关键作用，其中 $p$ 是素数且 $p \equiv 1 \bmod 4$。这样的 $x_0$ 被称为Størmer数，以挪威数学家和天文学家Carl Størmer（1874-1957）命名。本文建立了 $x_0 \in \mathbb{N}$ 成为某个素数 $p \equiv 1 \bmod 4$ 的Størmer数的充分必要条件。Størmer对这些数的研究主要源于他对将 $\pi$ 表示为Gregory-MacLaurin级数 $\arctan(1/x)$ 某些值的有限线性组合的恒等式研究。由于1900年时已知的 $\pi$ 位数不到600位，逼近 $\pi$ 是一个重要课题。Størmer在1896年发现的一个恒等式被Kanada及其团队在2002年用于计算出1.24万亿位 $\pi$。

研究背景与动机

问题背景

1. 费马二平方定理的重要性：该定理被G.H. Hardy称为"算术中最精美的定理之一"，它指出素数 $p$ 能表示为两个正整数的平方和当且仅当 $p \equiv 1 \bmod 4$。
2. 二次同余的核心作用：大多数费马二平方定理的证明都涉及研究二次同余方程 $x^2 \equiv -1 \bmod p$。例如，Hermite和Serret在1848年独立给出的证明假设已知最小剩余解 $x_0$，然后使用欧几里得算法产生求解正整数 $a, b$ 满足 $p = a^2 + b^2$ 的算法。
3. 历史意义：虽然Girard在1632年陈述了该定理，Fermat在1640年宣布了证明，但第一个有记录的存在性证明是Euler在1749年给出的，唯一性证明则由Gauss在1801年完成。

研究动机

1. 理论需求：当素数 $p$ 很大时，找到 $x_0$ 在实践中是困难的，需要系统的判定条件。
2. 历史价值：Carl Størmer的工作连接了数论和 $\pi$ 的逼近，具有重要的数学史意义。
3. 应用价值：Størmer的恒等式至今仍在高精度 $\pi$ 计算中发挥作用。

核心贡献

1. 建立了Størmer数的完整刻画：给出了正整数 $x_0$ 成为某个素数的Størmer数的充分必要条件（定理2.3）
2. 证明了Størmer函数的单射性：证明了每个Størmer数对应唯一的素数（定理2.1）
3. 提供了Størmer数自然密度的启发性证明：支持Everest和Harman猜想，即Størmer数的自然密度为 $\ln 2$
4. 系统阐述了Størmer数与Gregory数的联系：展示了Størmer的 $\pi$ 逼近理论及其现代应用
5. 提供了丰富的计算实例：通过具体例子说明了理论的实际应用

方法详解

任务定义

输入：正整数 $x_0$输出：判断 $x_0$ 是否为某个素数 $p \equiv 1 \bmod 4$ 的Størmer数，如果是，找到对应的素数 $p$约束：$p$ 必须是形如 $4n+1$ 的素数，且 $1 < x_0 \leq (p-1)/2$，$x_0^2 \equiv -1 \bmod p$

核心理论框架

定义2.1（Størmer数）

设 $p = 4n+1$ 是素数，正整数 $x_0$ 满足：

• (i) $1 < x_0 \leq \frac{p-1}{2}$
• (ii) $x_0^2 \equiv -1 \bmod p$

则称 $x_0$ 为 $p$ 的Størmer数，记作 $S(p) = x_0$。

关键定理

定理2.1（单射性）：Størmer函数 $S: P_{4n+1} \to \mathbb{N}$ 是单射的。

证明思路：假设 $S(p_1) = S(p_2)$ 但 $p_1 \neq p_2$，不失一般性设 $p_1 > p_2$。由于 $p_j | (x_0^2 + 1)$，得到 $p_1 p_2 | (x_0^2 + 1)$，即： $x_0^2 + 1 \geq p_1 p_2$

另一方面，由于 $x_0 \leq (p_2-1)/2$ 和 $p_2 \leq p_1 - 2$，可以推导出： $x_0^2 + 1 < p_1 p_2$

这产生矛盾，因此 $p_1 = p_2$。

定理2.2（素数性质）：若 $x_0^2 + 1 = 2^r p_1^{r_1} \cdots p_m^{r_m}$，则每个 $p_j \equiv 1 \bmod 4$。

定理2.3（充分必要条件）：设 $n = x_0^2 + 1$ 的素因数分解为 $2^r p_1^{r_1} \cdots p_m^{r_m}$，其中 $p_1 < \cdots < p_m$。则存在素数 $p$ 使得 $x_0$ 为其Størmer数当且仅当 $2x_0 + 1 \leq p_m$。若条件满足，则 $p = p_m$ 且 $S(p_m) = x_0$。

技术创新点

1. 统一的判定框架：将Størmer数的判定问题转化为对 $x_0^2 + 1$ 的素因数分解的分析
2. 构造性证明：不仅给出判定条件，还提供了找到对应素数的具体方法
3. 连接数论与分析：建立了Størmer数与Gregory级数、$\pi$ 的逼近之间的深层联系

实验设置

数据展示

论文提供了两个重要的数据表：

表1：素数与对应Størmer数的有序对

(5,2), (13,5), (17,4), (29,12), (37,6), (41,9), (53,23)
(61,11), (73,27), (89,14), (97,22), (101,10), (109,33), (113,15)
...

表2：按递增顺序排列的前若干个Størmer数

1, 2, 4, 5, 6, 9, 10, 11, 12, 14, 15, 16, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 33, 34, 35, 36, 37, 39, 40, 42, ...

验证方法

论文通过具体例子验证理论结果：

例子1：$x_0 = 3$

• $x_0^2 + 1 = 10 = 2 \times 5$
• $2x_0 + 1 = 7 > 5$
• 因此3不是Størmer数

例子2：$x_0 = 15$

• $x_0^2 + 1 = 226 = 2 \times 113$
• $2x_0 + 1 = 31 \leq 113$
• 因此 $S(113) = 15$

实验结果

自然密度分析

论文给出了Størmer数自然密度的数值证据：

这些数据支持Everest-Harman猜想：$\lim_{n \to \infty} \frac{|\{S(p) | p \leq n\}|}{n} = \ln 2 \approx 0.693$

启发性密度证明

论文提供了一个基于概率论的启发性证明：

对于 $x_0 \in \mathbb{N}$，如果 $x_0 = S(p)$，则 $2x_0 + 1 \leq p \leq x_0^2 + 1$。假设每个整数成为Størmer数的"概率"为 $2/(p-1)$，则 $x_0$ 成为Størmer数的概率约为：

$\sum_{2x_0+1 \leq p \leq x_0^2+1, p \equiv 1 \bmod 4} \frac{2}{p-1} \sim \sum_{2x_0 \leq p \leq x_0^2} \frac{1}{p} \sim \ln 2$

使用Merten估计得到最终结果。

Gregory数与π的逼近

理论基础

Størmer建立了Størmer数与Gregory数 $t_x := \arctan(1/x)$ 之间的联系：

定理5.1（Størmer定理）：如果 $x \in \mathbb{N}$ 不是Størmer数，则 $t_x$ 可以唯一表示为 $t_n$ 的有限线性组合，其中每个 $n$ 都是Størmer数且 $n < x$。

历史应用

1. Machin公式（1706）：$\frac{\pi}{4} = 4\arctan\frac{1}{5} - \arctan\frac{1}{239}$
2. Euler公式：$\frac{\pi}{4} = 5\arctan\frac{1}{7} + 2\arctan\frac{79}{3}$
3. Størmer公式（1896）：$\frac{\pi}{4} = 44\arctan\frac{1}{57} + 7\arctan\frac{1}{239} - 12\arctan\frac{1}{682} + 24\arctan\frac{1}{12943}$

现代意义

Størmer的1896年恒等式被Kanada团队在2002年用于计算1.24万亿位π，展现了经典数学理论的持久价值。

相关工作

历史发展

1. 费马二平方定理：从Girard（1632）到Fermat（1640）到Euler（1749）的发展历程
2. Smith方法（1855）：使用欧几里得算法和连分数理论的优雅存在性证明
3. Hermite-Serret算法（1848）：基于已知Størmer数构造平方和分解的算法

相关研究

1. Conway-Guy定义：略有不同的Størmer数定义，包含n=1的情况
2. Lehmer工作：研究arccot函数的恒等式
3. Todd定理：关于arctan函数不可约性的刻画

结论与讨论

主要结论

1. 完整刻画：建立了Størmer数的充分必要判定条件，解决了一个长期存在的问题
2. 理论统一：将Størmer数、费马二平方定理、Gregory级数和π的逼近统一在一个理论框架内
3. 历史价值：系统梳理了Carl Størmer的数学贡献及其现代意义

局限性

1. 密度证明：自然密度的证明仍是启发性的，缺乏严格的数学证明
2. 计算复杂性：对于大的 $x_0$，判定其是否为Størmer数需要分解 $x_0^2 + 1$，这在计算上可能是困难的
3. 应用范围：主要局限于理论数学，实际应用相对有限

未来方向

1. 严格密度证明：寻求Everest-Harman猜想的严格数学证明
2. 算法优化：开发更高效的Størmer数判定和计算算法
3. 推广研究：探索类似结构在其他数学分支中的应用

深度评价

优点

1. 理论完整性：提供了Størmer数的完整数学刻画，填补了理论空白
2. 历史价值：系统介绍了Carl Størmer的生平和贡献，具有重要的数学史意义
3. 连接性强：巧妙地连接了数论、分析学和π的计算等多个数学分支
4. 实例丰富：通过大量具体例子说明理论的应用，增强了可读性
5. 写作清晰：论文结构清晰，从基本定义到深入应用层层递进

不足

1. 计算效率：对于实际应用，素因数分解的计算复杂性可能成为瓶颈
2. 密度证明：自然密度的证明缺乏数学严谨性，仅为启发性论证
3. 应用局限：虽然理论完整，但现代实用价值相对有限
4. 创新程度：主要是对已有理论的系统化整理，原创性贡献相对较少

影响力

1. 理论贡献：为数论研究提供了新的工具和视角
2. 教育价值：适合作为数论和数学史的教学材料
3. 启发意义：可能激发对相关数学结构的进一步研究
4. 可复现性：理论结果清晰，计算实例易于验证

适用场景

1. 数论研究：为研究二次同余和素数性质提供新方法
2. 数学教育：作为连接不同数学分支的优秀教学案例
3. 数学史研究：为理解19-20世纪数学发展提供重要资料
4. 算法设计：为相关计算问题提供理论基础

参考文献

论文引用了21篇重要文献，涵盖了从经典数论教科书到现代计算数学的广泛范围，体现了研究的深度和广度。关键参考文献包括Andrews的《数论》、Conway和Guy的《数之书》、以及Størmer的原始论文等。