Let $p$ be a prime $\equiv 3$ mod 4, $p>3$, and suppose that 10 has the order $(p-1)/2$ mod p. Then $1/p$ has a decimal period of length $(p-1)/2$. We express the frequency of each digit $0,\ldots,9$ in this period in terms of the class numbers of two imaginary quadratic number fields. We also exhibit certain analogues of this result, so for the case that 10 is a primitive root mod $p$ and for octal digits of $1/p$.
论文ID : 2510.07873标题 : On the decimal digits of 1 / p 1/p 1/ p 作者 : Kurt Girstmair (University of Innsbruck)分类 : math.NT (Number Theory)发表时间 : 2025年10月28日 (arXiv v4)论文链接 : https://arxiv.org/abs/2510.07873 设p p p 是一个满足p ≡ 3 ( m o d 4 ) p \equiv 3 \pmod{4} p ≡ 3 ( mod 4 ) 且p > 3 p > 3 p > 3 的素数,假设10在模p p p 下的阶为( p − 1 ) / 2 (p-1)/2 ( p − 1 ) /2 。则1 / p 1/p 1/ p 的十进制周期长度为( p − 1 ) / 2 (p-1)/2 ( p − 1 ) /2 。本文用两个虚二次数域的类数来表达该周期中每个数字0 , … , 9 0,\ldots,9 0 , … , 9 的频率。文章还展示了该结果的类似情形,包括10是模p p p 的原根的情况以及1 / p 1/p 1/ p 的八进制数字的情况。
问题背景 : 自Girstmair在1994年首次建立1 / p 1/p 1/ p 的数字与类数之间的联系以来,这一领域产生了众多研究。现有研究主要关注1 / p 1/p 1/ p 或更一般的m / p m/p m / p (其中m m m 是模p p p 的二次剩余)的十进制数字分布。问题重要性 : 理解素数倒数的数字分布不仅具有纯数学意义,还与代数数论中的类数理论密切相关。这种联系为计算大素数的数字频率提供了有效方法。现有局限性 : 已知结果只能给出数字的平均值,即∑ j = 1 ( p − 1 ) / 2 a j = 9 ( p − 1 ) / 4 + 9 h 1 / 2 \sum_{j=1}^{(p-1)/2} a_j = 9(p-1)/4 + 9h_1/2 ∑ j = 1 ( p − 1 ) /2 a j = 9 ( p − 1 ) /4 + 9 h 1 /2 ,其中h 1 h_1 h 1 是虚二次数域Q ( − p ) \mathbb{Q}(\sqrt{-p}) Q ( − p ) 的类数,但无法给出每个数字的具体频率。研究动机 : 本文旨在精确确定周期中每个数字0 , … , 9 0,\ldots,9 0 , … , 9 的频率,并用两个虚二次数域的类数来表达这些频率。主要定理 : 提出了Theorem 1,给出了在特定条件下每个数字频率的精确公式,涉及两个类数h 1 h_1 h 1 和h 2 h_2 h 2 。扩展结果 : 提供了10是原根情况下的类似结果(Theorem 2)和八进制数字的对应结果(Theorem 3)。计算优势 : 为大素数(如10 15 10^{15} 1 0 15 量级)的数字频率计算提供了有效方法,避免了直接计算的困难。理论深化 : 首次在文献中给出涉及两个类数的数字分布结果。给定素数p ≡ 3 ( m o d 4 ) p \equiv 3 \pmod{4} p ≡ 3 ( mod 4 ) ,p > 3 p > 3 p > 3 ,且10在模p p p 下的阶为( p − 1 ) / 2 (p-1)/2 ( p − 1 ) /2 ,确定m / p m/p m / p (m m m 为二次剩余)的十进制展开式中每个数字的频率。
对于m / p m/p m / p 的十进制展开:
m / p = ∑ j = 1 ∞ a j 10 − j m/p = \sum_{j=1}^{\infty} a_j 10^{-j} m / p = ∑ j = 1 ∞ a j 1 0 − j
其中a j ∈ { 0 , … , 9 } a_j \in \{0,\ldots,9\} a j ∈ { 0 , … , 9 } ,序列( a 1 , … , a ( p − 1 ) / 2 ) (a_1,\ldots,a_{(p-1)/2}) ( a 1 , … , a ( p − 1 ) /2 ) 构成周期。
Lemma 1 : 设θ b ( k ) = b ( k ) p − ( b k ) p p \theta_b(k) = \frac{b(k)_p - (bk)_p}{p} θ b ( k ) = p b ( k ) p − ( bk ) p ,则
a j = θ 10 ( m ⋅ 10 j − 1 ) a_j = \theta_{10}(m \cdot 10^{j-1}) a j = θ 10 ( m ⋅ 1 0 j − 1 )
Lemma 2 : 基于Berndt的结果,给出了二次剩余和非剩余在特定区间内的分布差异。
Theorem 1 : 在上述设置下,对于k ∈ { 0 , … , 9 } k \in \{0,\ldots,9\} k ∈ { 0 , … , 9 } :
n k = 1 2 ( ⌊ ( k + 1 ) p 10 ⌋ − ⌊ k p 10 ⌋ + δ k ) n_k = \frac{1}{2}\left(\left\lfloor\frac{(k+1)p}{10}\right\rfloor - \left\lfloor\frac{kp}{10}\right\rfloor + \delta_k\right) n k = 2 1 ( ⌊ 10 ( k + 1 ) p ⌋ − ⌊ 10 k p ⌋ + δ k )
n 9 − k = n k − δ k , k = 0 , … , 4 n_{9-k} = n_k - \delta_k, \quad k = 0,\ldots,4 n 9 − k = n k − δ k , k = 0 , … , 4
其中δ k \delta_k δ k 的值依赖于p m o d 8 p \bmod 8 p mod 8 :
若p ≡ 3 ( m o d 8 ) p \equiv 3 \pmod{8} p ≡ 3 ( mod 8 ) :δ 0 = 0 \delta_0 = 0 δ 0 = 0 δ 1 = δ 4 = 3 h 1 / 2 − h 2 / 4 \delta_1 = \delta_4 = 3h_1/2 - h_2/4 δ 1 = δ 4 = 3 h 1 /2 − h 2 /4 δ 2 = − δ 1 \delta_2 = -\delta_1 δ 2 = − δ 1 δ 3 = 3 h 1 / 2 + h 2 / 4 \delta_3 = 3h_1/2 + h_2/4 δ 3 = 3 h 1 /2 + h 2 /4 若p ≡ 7 ( m o d 8 ) p \equiv 7 \pmod{8} p ≡ 7 ( mod 8 ) :δ 0 = h 1 − h 2 / 2 \delta_0 = h_1 - h_2/2 δ 0 = h 1 − h 2 /2 δ 1 = h 2 / 4 \delta_1 = h_2/4 δ 1 = h 2 /4 δ 2 = 3 h 2 / 4 \delta_2 = 3h_2/4 δ 2 = 3 h 2 /4 δ 3 = δ 4 = − δ 1 \delta_3 = \delta_4 = -\delta_1 δ 3 = δ 4 = − δ 1 区间分析法 : 通过建立数字与特定区间内整数的对应关系,将数字频率问题转化为区间内二次剩余计数问题。类数表示 : 利用Berndt的深刻结果,将二次剩余的分布差异用虚二次数域的类数表达。对称性利用 : 巧妙利用p − l p-l p − l 与l l l 的二次剩余性质对偶关系,建立n k n_k n k 与n 9 − k n_{9-k} n 9 − k 的关系。论文提供了几个具体的数值验证:
p = 67 p = 67 p = 67 : h 1 = 1 h_1 = 1 h 1 = 1 , h 2 = 18 h_2 = 18 h 2 = 18 1 / 67 = 0. 014925373134328358208955223880597 ‾ 1/67 = 0.\overline{014925373134328358208955223880597} 1/67 = 0. 014925373134328358208955223880597 理论预测的频率与实际计算完全吻合 p = 79 p = 79 p = 79 : 展示了更一般情况下的应用涉及多个分数1 / 79 1/79 1/79 , 9 / 79 9/79 9/79 , 2 / 79 2/79 2/79 的周期分析 p = 47 p = 47 p = 47 : 原根情况的验证在广义黎曼假设下,满足条件的素数密度约为A / 2 = 0.186977 … A/2 = 0.186977\ldots A /2 = 0.186977 … (其中A A A 是Artin常数) 约37%的素数p ≡ 3 ( m o d 4 ) p \equiv 3 \pmod{4} p ≡ 3 ( mod 4 ) 满足10的阶为( p − 1 ) / 2 (p-1)/2 ( p − 1 ) /2 的条件 精确匹配 : 所有数值验证都完美验证了理论公式的正确性计算效率 : 对于大素数,使用类数计算频率比直接计算快得多分布规律 : 揭示了数字频率的内在规律,如p ≡ 3 ( m o d 8 ) p \equiv 3 \pmod{8} p ≡ 3 ( mod 8 ) 时n 3 > n 6 n_3 > n_6 n 3 > n 6 原根情况 : Theorem 2处理了10是原根的情形,给出了奇偶位置数字的不同分布八进制情况 : Theorem 3将结果推广到八进制展开,涉及Q ( − 2 p ) \mathbb{Q}(\sqrt{-2p}) Q ( − 2 p ) 的类数论文详细回顾了该领域的发展历程:
Girstmair (1994)首次建立联系 后续研究包括Hirabayashi, Murty & Thangadurai, Mizuno等的工作 本文是首个涉及两个类数的结果 成功建立了1 / p 1/p 1/ p 数字频率与两个虚二次数域类数的精确关系 提供了计算大素数数字分布的有效方法 揭示了数字分布的深层数论结构 仅适用于特定条件的素数(p ≡ 3 ( m o d 4 ) p \equiv 3 \pmod{4} p ≡ 3 ( mod 4 ) 且10的阶为( p − 1 ) / 2 (p-1)/2 ( p − 1 ) /2 ) 对于p ≡ 1 ( m o d 4 ) p \equiv 1 \pmod{4} p ≡ 1 ( mod 4 ) 的情况,虽有公式但难以用类数解释 依赖于类数的高效计算 扩展到更一般的素数情况 研究其他进制下的类似结果 探索与椭圆曲线等其他代数结构的联系 理论深度 : 将初等的数字分布问题与深刻的代数数论联系起来计算实用性 : 为大素数的数字分布计算提供了实际可行的方法结果完整性 : 不仅给出主要结果,还提供了多个扩展和类比验证充分 : 通过多个具体例子验证了理论的正确性适用范围限制 : 结果仅适用于满足特定条件的素数证明依赖性 : 关键步骤依赖于Berndt的深刻结果,本文贡献相对初等推广困难 : 向其他情况的推广面临技术障碍学术价值 : 为数字分布与类数理论的联系提供了新的视角计算意义 : 在大素数计算中具有实际应用价值启发性 : 可能启发更多相关方向的研究大素数的数字分布快速计算 代数数论中类数的应用研究 数字序列分析和密码学相关应用 论文引用了16篇相关文献,主要包括:
B.C. Berndt的经典二次剩余定理 作者自己的早期工作 近年来该领域的主要进展 类数计算的相关方法 这篇论文在数论领域具有重要意义,它不仅解决了一个具体的数字分布问题,更重要的是揭示了初等数论与代数数论之间的深刻联系,为相关研究提供了新的工具和视角。