短区间上的Erdös-Kac定理
Weighted Erdös-Kac Theorem in Short Intervals
DOI: 10.12677/pm.2025.154116, PDF, HTML, XML,   
作者: 刘新颖:青岛大学数学与统计学院,山东 青岛
关键词: Euler函数Erdös-Kac定理短区间Euler Function Erdös-Kac Theorem Short Intervals
摘要: φ( n ) 是Euler函数,本文给出了短区间上算术序列的Erds-Kac型定理,其中该算术序列与 φ( n ) 分布有关。
Abstract: Assuming φ( n ) is an Euler function, this article provides the Erds-Kac type theorem for arithmetic sequences on short intervals, where the arithmetic sequence is related to the Euler function’s distribution.
文章引用:刘新颖. 短区间上的Erdös-Kac定理[J]. 理论数学, 2025, 15(4): 127-137. https://doi.org/10.12677/pm.2025.154116

1. 引言

Erdös-Kac定理在概率数论中有非常重要的地位。它为数论和概率论之间架起了桥梁,展示了整数的素因子数量具有统计性质,且趋近于泊松分布。这个定理不仅为数学理论提供了新的视角,还对现代计算数学和密码学等领域产生了深远的影响。该定理证明了对于任给的非负整数n,它的不同素因子的个数在 { n:nx } 上的分布近似于期望为 log 2 x:=loglogx ,方差为 ( log 2 x ) 1/2 的高斯分布,即对于任给的 λ

1 x nx ω( n ) log 2 xλ ( log 2 x ) 1/2 1Φ( λ ) x( Φ( λ ):= 1 2π λ e τ 2 /2 dτ )

其中 ω( n ) n的不同素因子的个数, Φ( λ ) 是正态分布函数。

随之,数学家们研究了加权Erdös-Kac定理。Elliott在文章[1]中给出了一个以 d ( n ) α 为权的Erdös-Kac型定理,其中 d( n ) 是除数函数。刘奎和吴杰[2]将Elliott的结果推广到短区间的情况,即

1 D α ( x,y ) x<nx+y ω( n ) log 2 xλ ( log 2 x ) 1/2 d ( n ) α Φ( λ ) D α ( x,y ):= x<nx+y d ( n ) α

对于 x x 7/ 12 +ε yx 是一致的,其中隐藏的常数仅取决于 α 和. ε 。刘晓莉和杨志善[3]建立了一个在短区间上以 a Κ l ( n 2 )( l + ) 加权的Erdös-Kac型,其中K是高斯域, a Κ ( n ) [ i ] 中范数为n的非零整理想的个数。王丹[4]又给出了短区间某些算术序列上以 d ( n ) α 加权的Erdös-Kac型定理。

Selberg-Delange方法是研究Erdös-Kac定理的重要工具之一。在1954年到1971年之间,Selberg [5]和Delange [6]利用与算术函数相关的Dirichlet级数的解析性质发展了一种相当普遍的方法,这种方法现在被称为Selberg-Delange方法。

f( n ) 是一个算术函数,用下式表示对应的Dirichlet级数:

F( s ):= n1 f( n ) n s s>1

对于 s>1 ,假设 F( s ) 有如下分解形式

F( s )=G( s;z )ζ ( s ) z

其中 ζ( s ) 是Riemann ζ( s ) -函数, G( s;z ) 是亚纯函数, z 。在上述假设下,Selberg和Delange利用解析方法给出了算术函数 f( n ) 均值估计

S f ( x ):= nx f( n ) 3

的一个渐近公式。Tenenbaum在文章[7]中对这一理论有较为详细的阐述。

近年来,Selberg-Delange方法被多位数学家进行了适当的改进,在不同的方向和背景下有许多应用。本文主要利用Labihi和Raoujand [8]的推广Selberg-Delange方法,建立了短区间上新的Erds-Kac型定理。

欧拉函数 φ( n ) 是数论中的一个基本函数,它表示小于或等于n的正整数中与n互质的数的个数,即那些与n的最大公约数为1的整数个数,即对于正整数n,我们有

φ( n )=| { 1kn:gcd( k,n )=1 } |

其中 gcd( k,n ) 表示kn的最大公约数。当n很大时,欧拉函数 φ( n ) 可以通过以下近似来得到:

φ( n )n p|n ( 1 1 p )

其中 p|n 表示所有n的素因子。根据这个公式可以看出,欧拉函数的值大致与n的素因子分布有关,特别是素因子越多,欧拉函数的值相对越小。

定义

U( x,y ):= x<φ( n )x+y 1

其中 φ( n ) 是Euler函数。本文的主要结果如下:

定理1.1. ε>0 ,则对任意实数 λ ,有

1 U( x,y ) x<φ( n )x+y ω( n ) log 2 xλ ( log 2 x ) 1/2 1 =Φ( λ )+O( 1 log 2 x )

x x 19/ 24 +ε yx 一致成立。特别地,误差项是最优的。

为了证明误差项的最优性,我们需要建立一个在短区间上的Landau素数定理。对于 k + ,定义

U k ( x,y ):= x<φ( n )x+y ω( n )=k 1

对上述和式我们有以下的结果:

定理1.2. ε>0 ,有

U k ( x,y )= y logx ( log 2 x ) k1 ( k1 )! ×Y( l 1 ( 1 ), l 2 ( 1 ), k1 log 2 x )+O( ( log 2 x ) 2 klogx + k1 ( log 2 x ) 2 )

x3 x 19/ 24 +ε yx 1k log 2 x 一致成立,其中

Y( l 1 ( 1 ), l 2 ( 1 ),z )= 1 τ( z+1 ) p ( 1 1 p ) z ( 1+ z p ( 1 p 1 ) 2 )

且隐藏常数与 ε>0 有关。

2. 一些引理

设函数 f: * 是一个算术函数, g: * ( 0, ) 满足

lim n g( n )= limsup logn logg( n ) =κ

以下均设 s=σ+it 。Dirichlet级数

F( s )= 1 f( n )g ( n ) s

σ σ c 处收敛,在 σ σ a 处绝对收敛,其中 σ a σ c +κ 。下面给出两个定义:

定义1 ([8]) 如果以下条件成立,我们说横坐标有限的收敛Dirichlet级数

F( s )= 1 f( n )g ( n ) s

具有性质 Ρ( A,M, M 1 , M 2 ,δ, l 1 , l 2 )

(1) g( n ) 是实值乘法函数,使得 g:[ 1, ) g( n ) 随着 n 趋于无穷大,且存在 κ0 的常数,使得

lim n sup logn logg( n ) =κ

(2) 在 σ>1 时,可以将 F( s ) 分解为

F( s )=G( s, l 1 ( s ), l 2 ( s ) )ζ ( s ) l 1 ( s ) ζ ( 2s ) l 2 ( s )

(3) 函数 G( s, l 1 ( s ), l 2 ( s ) ) σ>1/2 的区域内有解析延拓,满足

| G( s, l 1 ( s ), l 2 ( s ) ) |M ( 3+| t | ) max{ δ( 1σ ),0 } log A ( 3+| t | )

(4) 函数 l 1 ( s ) l 2 ( s ) σ>1/2 的区域内是解析的,并且在 σ[ 1/2 ,2+ε ] 满足

| l 1 ( s ) | M 1 | l 2 ( s ) | M 2

定义2 ([8]) 设Dirichlet级数

F( s )= 1 f( n )g ( n ) s

具有性质 Ρ( A,M, M 1 , M 2 ,δ, l 1 , l 2 ) 。如果存在正实数序列 ( f + ˜ ( n ) ) n1 ,使得对所有 n1 ,有 | f( n ) | f + ˜ ( n ) 其中Dirichlet级数

F + ˜ ( s ):= 1 f + ˜ ( n )g ( n ) s

Ρ( A,M, M 1 , M 2 ,δ, l ˜ 1 , l ˜ 2 ) 型,那么我们称 F( s ) Τ( A,M, M 1 , M 2 ,δ, l 1 , l 2 , l ˜ 1 , l ˜ 2 ) 型。

首先,我们给出Labihi和Raoujand关于短区间Selberg-Delange方法的一般结果,该结果在定理1.1的证明中起着关键作用。

引理1 ([4]) 设Dirichlet级数

F( s )= 1 f( n )g ( n ) s

Τ( A,M, M 1 , M 2 ,δ, l 1 , l 2 , l ˜ 1 , l ˜ 2 ) 型。对于任意的 ε>0 ,我们有

n1 x<g( n )x+y f( n ) =y ( logx ) l 1 ( 1 )1 { λ( l 1 ( 1 ), l 2 ( 1 ) )+O( log 2 x logx ) }

该结论对于 x2 x θ+ε yx θ:=1 1 ( 24/5 )+2δ | l 1 ( s ) | M 1 | l 2 ( s ) | M 2 成立,其中,

λ( l 1 ( 1 ), l 2 ( 1 ) ):= G( 1, l 1 ( 1 ), l 2 ( 1 ) )ζ ( 2 ) l 2 ( 1 ) Γ( l 1 ( 1 ) ) G( s, l 1 ( s ), l 2 ( s ) ):=F( s )ζ ( s ) l 1 ( s ) ζ ( 2s ) l 2 ( s )

O项中的隐含常数仅取决于 A , M 1 , M 2 , δ , ε

对于定理1.1的证明,我们使用加性函数理论中的高斯误差律和Berry-Esseen不等式。

引理2 ([9]) 设 f( n ) 是一个加法函数,满足 f( m 1 m 2 )=f( m 1 )+f( m 2 ) ,其中 ( m 1 , m 2 )=1 ,假设

F( n )= p<n f ( p ) 2 p

收敛,那么

f( m )< p<m f( p ) p +ω 2F( m )

的密率对于任给的实数 ω ,都等于

1 π ω exp( y 2 )dy

假设 F( x ) 是一个满足 F( )=0 F( )=1 的分布函数,我们将 F( x ) 的特征函数定义为

f( τ ):= e iτx dF( x )

引理3 ([7]) 设FG是分布函数,分别具有特征函数fg。设G是可微的, G 上有界。那么对于所有的 T>0

FG 16 G T +6 T T | f( τ )g( τ ) τ |dτ

其中对于任何实值函数 H( x ) ,定义 H := sup λ | H( λ ) |

φ( n ) 显然满足 lim n φ( n )= limsup logn logφ( n ) =1 ,为了证明定理1.2,我们取 g( n )=φ( n ) ,故可以给出关于 z ω( n ) 的如下渐近公式。

引理4 B>0 | z |B ε>0 ,那么

n1 x<φ( n )x+y z ω( n ) =y ( logx ) z1 { zY( l 1 ( 1 ), l 2 ( 1 ),z )+O( log 2 x logx ) }

x2 x 19/ 24 +ε yx 一致成立,其中,

特别地,

U( x,y )=y{ Y( l 1 ( 1 ), l 2 ( 1 ),1 )+O( log 2 x logx ) }

x2 x 19/ 24 +ε yx 一致成立,其中

Y( l 1 ( 1 ), l 2 ( 1 ),1 )= Π p ( 1 1 p )( 1+ 1 p ( 1 p 1 ) 2 )

证明 因为函数 n z ω( n ) 是可乘的,对于 s>1 ,我们可以得到如下式子:

F z ( s ):= n1 z ω( n ) φ ( n ) s =ζ ( s ) l 1 ( s ) ζ ( 2s ) l 2 ( s ) G( s, l 1 ( s ), l 2 ( s ) )

其中, l 1 ( s )=z l 2 ( s )= ( z z 2 )/2 | z |B

G( s, l 1 ( s ), l 2 ( s ) ):= Π p ( 1 1 p s ) z ( 1 1 p 2s ) ( z z 2 )/2 ( 1+ z p s ( 1 p s ) ( 1 p 1 ) s )

函数 G( s, l 1 ( s ), l 2 ( s ) ) { s:s>1/3 } 上是解析的,并且对于所有 ε>0 G( s, l 1 ( s ), l 2 ( s ) ) ε 1

{ s:s>1/3 +ε }{ s:s>1/2 }

上是解析的,所以对于 1/2 σ2+ε F z ( s ) 是一类Dirichlet级数,满足

Τ( 0,M( B ),B,1/8 ,0,z, ( z z 2 )/2 , l ˜ 1 , l ˜ 2 )

利用引理1,我们得到了 z ω( n ) 的均值渐近公式。

3. 定理1.2的证明

定义

U k ( x,y ):= ω( n )=k x<φ( n )x+y 1

我们可以得到

x<φ( n )x+y z ω( n ) = x<φ( n )x+y ω( n )=k z k = k x<φ( n )x+y z k = k U k ( x,y ) z k

通过柯西积分公式,

U k ( x,y )= 1 2πi | z |=r ( x<φ( n )x+y z ω( n ) ) dz z k+1 (1)

其中, r= k log 2 x 。通过引理4,

n1 x<φ( n )x+y z ω( n ) =y ( logx ) z1 { zY( l 1 ( 1 ), l 2 ( 1 ),z )+O( log 2 x logx ) }

代入(1),我们可以得到

U k ( x,y )= 1 2πi | z |=r ( x<φ( n )x+y z ω( n ) ) dz z k+1 = 1 2πi | z |=r zy ( logx ) z1 Y( l 1 ( 1 ), l 2 ( 1 ),z ) 1 z k+1 dz + 1 2πi | z |=r y ( logx ) z1 O( log 2 x logx ) 1 z k+1 dz = y logx 1 2πi | z |=r ( logx ) z Y( l 1 ( 1 ), l 2 ( 1 ),z ) z k dz +O( y log 2 x ( logx ) 2 | z |=r ( logx ) z | z | k+1 | dz | ) = y logx I k ( x,r )+O( y log 2 x ( logx ) 2 | z |=r ( logx ) z | z | k+1 | dz | ),

其中, I k ( x,r )= 1 2πi | z |=r ( logx ) z Y( l 1 ( 1 ), l 2 ( 1 ),z ) z k dz

对于 U k ( x,y ) 的误差项,利用变量替换 t=k( 1cosθ ) ,我们有以下估计:

| z |=r ( logx ) z | z | k+1 | dz | = | z |=r ( logx ) rcosθ r k+1 rdθ = | z |=r ( logx ) rcosθ r k dθ = ( log 2 x k ) k 0 2π ( logx ) kcosθ log 2 x dθ = ( log 2 x k ) k 0 2π e kcosθ dθ ( log 2 x k ) k ( 0 π 2 e kcosθ dθ +1 ) ( log 2 x k ) k ( e k k 0 k e t t 1 2 dt+1 ) ( log 2 x ) k k! .

接下来将分别计算 I k ( x,r ) k=1 k2 的值。

k=1 时,

I 1 ( x,r )= 1 2πi | z |=r e z log 2 x Y( l 1 ( 1 ), l 2 ( 1 ),z ) z dz =Y( l 1 ( 1 ), l 2 ( 1 ),0 )=1

所以, U 1 ( x,y )= y logx { 1+O( ( log 2 x ) 2 logx ) }

对于 k2 ,通过柯西积分公式我们可以得到 I k ( x,r )= I k ( x, r 0 ) r 0 = k1 log 2 x

Y( l 1 ( 1 ), l 2 ( 1 ),z ) z= r 0 处泰勒展开可以得到

Y( l 1 ( 1 ), l 2 ( 1 ),z )=Y( l 1 ( 1 ), l 2 ( 1 ), r 0 )+ Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 )( z r 0 )+ r 0 z Y ( l 1 ( 1 ), l 2 ( 1 ),u )( zu )du =Y( l 1 ( 1 ), l 2 ( 1 ), r 0 )+ Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 )( z r 0 ) + ( z r 0 ) 2 0 1 ( 1t ) Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 +t( z r 0 ) )dt .

那么,

I k ( x,r )= Y( l 1 ( 1 ), l 2 ( 1 ), r 0 ) 2πi | z |= r 0 e z log 2 x z k dz + Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 ) 2πi | z |= r 0 e z log 2 x ( z r 0 ) z k dz + 1 2πi | z |= r 0 ( logx ) 2 ( z r 0 ) 2 0 1 ( 1t ) Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 +t( z r 0 ) )dt z k dz .

接下来我们依次估计这三部分。通过柯西积分公式,第一部分为

Y( l 1 ( 1 ), l 2 ( 1 ), r 0 ) 2πi | z |= r 0 e z log 2 x z k dz = ( log 2 x ) k1 ( k1 )! Y( l 1 ( 1 ), l 2 ( 1 ), r 0 )

同样地,第二部分为

Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 ) 2πi | z |= r 0 e z log 2 x ( z r 0 ) z k dz = Y ( l 1 ( 1 ), l 2 ( 1 ), r 0 )( ( log 2 x ) k2 ( k2 )! r 0 ( log 2 x ) k1 ( k1 )! )=0

对于 0<t1 | z |= r 0 ,我们有

| r 0 +t(z r 0 ) |=| r 0 (1t)+tz | r 0 (1t)+t| z |= r 0

因为在 | z |B Y( l 1 ( 1 ), l 2 ( 1 ),z ) 是解析的,所以存在一个正常数 C α ,使得

| Y ( l 1 ( 1 ), l 2 ( 1 ),z ) | C α

因此,那么第三部分利用变量替换 t=( k1 )( 1cosθ ) ,有

| z |= r 0 ( logx ) z | z r 0 | 2 | z | k | dz | = 0 2π e z log 2 x r 0 2 | e iθ 1 | 2 r 0 k r 0 dθ = r 0 ( k3 ) 0 2π e log 2 x k1 log 2 x cosθ | e iθ 1 | 2 dθ = r 0 ( k3 ) 0 2π e ( k1 )cosθ | e iθ 1 | 2 dθ r 0 ( k3 ) ( 0 π 2 e ( k1 )cosθ ( 1cosθ )dθ+π ) r 0 ( k3 ) e k1 ( k1 ) 3 2 ( 0 k1 e t t 1 2 dt+π ) ( log 2 x ) k3 ( k2 )! .

那么我们得到

U k ( x,y )= y ( log 2 x ) k1 ( k1 )!logx ( Y( l 1 ( 1 ), l 2 ( 1 ), k1 log 2 x ) )+O( y ( log 2 x ) k+1 k! ( logx ) 2 + y ( log 2 x ) k3 ( k2 )!logx )

4. 定理1.1的证明

定义

F x,y ( λ )= 1 U( x,y ) x<φ( n )x+y ω( n ) log 2 xλ ( log 2 x ) 1/2 1

φ x,y ( τ ) F x,y ( λ ) 的特征函数,

φ x,y ( τ ):= e iτλ d F x,y ( λ ) = 1 U( x,y ) x<φ( n )x+y exp{ iτ ω( n ) log 2 x log 2 x } = e iτT U( x,y ) x<φ( n )x+y e i( τ/T )ω( n ) ,

其中, T= log 2 x 。通过Berry-Esseen不等式,

F x,y ϕ 16 2π T +6 T T | φ x,y ( τ ) e τ 2 2 τ |dτ ,

接下来我们需要去证明 T T | φ x,y ( τ ) e τ 2 2 τ |dτ 1 T 。将 z= e it 引入引理4可得,

n1 x<φ( n )x+y e itω( n ) =y ( logx ) e it 1 { e it Y( l 1 ( 1 ), l 2 ( 1 ), e it )+O( log 2 x logx ) }

z=1 代入引理4可得,

U( x,y )= n1 x<φ( n )x+y 1 =y{ Y( l 1 ( 1 ), l 2 ( 1 ),1 )+O( log 2 x logx ) }

所以,

1 U( x,y ) x<φ( n )x+y e itω( n ) = ( logx ) e it 1 { A( e it )+O( log 2 x logx ) }

其中, A( z )= zY( l 1 ( 1 ), l 2 ( 1 ),z ) Y( l 1 ( 1 ), l 2 ( 1 ),1 ) 并且 A( 1 )=1 。令 t=τ/T

φ x,y ( τ )= ( logx ) e i τ T 1 e iτT { A( e i τ T )+O( log 2 x logx ) }

| t |1 时, cost12 ( t π ) 2 ,我们有

| ( logx ) e i τ T 1 e iτT |=| e T 2 ( e i τ T 1 ) e i τ T T 2 |= e T 2 ( cos τ T 1 ) e 2 ( τ π ) 2

因此,

± T 1 3 ±T | φ x,y ( τ ) e τ 2 2 τ |dτ ± T 1 3 ±T e 2 ( τ π ) 2 dτ 1 T

( logx ) 1 <| τ | T 1 3 时,由泰勒展开式,

A( e i τ T )=1+O( τ T ) e i τ T 1=i τ T 1 2 ( τ T ) 2 +O( ( τ T ) 3 )

那么,

( logx ) e i τ T 1 A( e i τ T ) e i( τT ) = e ( e iτ T 1 ) T 2 iτT A( e iτ T )= e τ 2 2 +O( τ 3 T ) { 1+O( | τ | T ) }= e τ 2 2 { 1+O( | τ |+ | τ | 3 T ) }

所以, φ x,y ( τ )= e τ 2 2 { 1+O( | τ |+ | τ | 3 T ) }{ 1+O( log 2 x logx ) } 。那么,

± 1 logx ± T 1 3 | φ x,y ( τ ) e τ 2 2 τ |dτ 1 logx T 1 3 ( e τ 2 2 1+ τ 2 T + log 2 x τlogx )dτ 1 T + ( log 2 x ) 2 logx 1 T

所以对于 | τ | 1 logx ,有

| τ ω( n ) log 2 x log 2 x | | τ |logx T exp{ iτ ω( n ) log 2 x log 2 x }=1+O( | τ |logx T ) ,

故有, φ x,y ( τ )=1+O( | τ |logx T ) 。又因为 e τ 2 2 =1+O( τ 2 ) ,所以,

1 logx 1 logx | φ x,y ( τ ) e τ 2 2 τ |dτ 1 logx 1 logx ( logx T +| τ | )dτ 1 T

5. 误差

下面我们将证明误差是合理的。令

R λ ( x,y ):= 1 U( x,y ) x<φ( m )x+y ω( n ) log 2 xλ ( log 2 x ) 1/2 1ϕ( λ ) R( x,y ):= sup λ R λ ( x,y )

其中 k=[ log 2 x ] θ=k log 2 x

U k ( x,y ) U( x,y ) = F x,y ( θ log 2 x ) F x,y ( θ 1 2 log 2 x )ϕ( θ log 2 x )ϕ( θ 1 2 log 2 x )+2R( x,y ) = θ 1 2 log 2 x θ log 2 x e τ 2 2 dτ+ 2R( x,y ) 1 2 π log 2 x +2R( x,y ).

另一方面,

U k ( x,y ) U( x,y ) ~ Y( l 1 ( 1 ), l 2 ( 1 ), k1 log 2 x ) Y( l 1 ( 1 ), l 2 ( 1 ),1 ) ( log 2 x ) k1 logx( k1 )! 1 π log 2 x

所以有,

R( x,y ) 1+o( 1 ) 2 π log 2 x 1 4 π log 2 x = 1+o( 1 ) 4 π log 2 x

参考文献

[1] Elliott, P.D.T.A. (2014) Central Limit Theorems for Classical Cusp Forms. The Ramanujan Journal, 36, 81-98.
https://doi.org/10.1007/s11139-013-9516-9
[2] Liu, K. and Wu, J. (2021) Weighted Erdős-Kac Theorem in Short Intervals. The Ramanujan Journal, 55, 1-12.
https://doi.org/10.1007/s11139-020-00343-1
[3] Liu, X.L. and Yang, Z.S. (2020) Weighted Erdős-Kac Type Theorems over Gaussian Field in Short Intervals. Acta Mathematica Hungarica, 162, 465-482.
https://doi.org/10.1007/s10474-020-01087-6
[4] Wang, D. (2024) Weighted Erdős-Kac Type Theorems in Short Intervals. Acta Arithmetica, 212, 255-267.
https://doi.org/10.4064/aa230129-22-8
[5] Sathe, L.G. (1945) On a Congruence Property of the Divisor Function. American Journal of Mathematics, 67, 397.
https://doi.org/10.2307/2371953
[6] Delange, H. (1959) Sur des formules dues à Atle Selberg. Bulletin des Sciences Mathématiques, 2, 101-111.
[7] Tenenbaum, G. (2015). Introduction to Analytic and Probabilistic Number Theory. American Mathematical Society, 163.
https://doi.org/10.1090/gsm/163
[8] Labihi, O. and Raouj, A. (2020) Estimation de certaines sommes courtes. Journal of Number Theory, 206, 231-249.
https://doi.org/10.1016/j.jnt.2019.06.012
[9] Erdös, P. and Kac, M. (1939) On the Gaussian Law of Errors in the Theory of Additive Functions. Proceedings of the National Academy of Sciences, 25, 206-207.
https://doi.org/10.1073/pnas.25.4.206