符号定义
沿用前一节中定义的符号, 可知:
S(A,P,z)=#{a∈A:p∤a,∀p∈P∩[2,z)},(2.1)
很明显当 P(z) 表示全体 P∩[2,z) 元素的乘积时 (2.1) 就可以化成:
S(A,P,z)=a∈A(a,P(z))=1∑1.
利用 Möbius 反演可知:
S(A,P,z)=d∣P(z)∑μ(d)∣Ad∣,(2.2)
因此在估计 S(A,P,z) 时我们需要对 ∣Ad∣ 的大小有一定的了解. 一般来讲我们会要求其满足如下性质:
存在 X>0 和积性函数 g(d) 使得: ∣Ad∣=g(d)X+r(d)其中当 p 为素数时 0<g(p)<1.
在本讲义的后续章节中我们都会默认 A 满足这条性质.
将这个渐近公式回代至 (2.2) 中, 便有:
S(A,P,z)=XV(z)+d∣P(z)∑μ(d)r(d),
其中 V(z) 表示下列乘积:
V(z)=d∣P(z)∑μ(d)g(d)=p∣P(z)∏(1−g(p)).
利用这一点, 我们就得到了著名的 Eratosthenes 筛法:
当 Ad 满足命题 2.1.1 时, 总存在 −1≤θ≤1 使得: S(A,P,z)=XV(z)+θd∣P(z)∑∣r(d)∣.
然而 z 较大时 Eratosthenes 筛法中产生的误差项通常以指数形式增长, 所以这个时候我们就需要引入一些高级工具来处理这个问题.
上下界的构造
我们连续使用两次 Buchstab 迭代式, 可知:
S(A,P,z)=∣A∣−p1<z∑∣Ap1∣+p1>p2p1<z∑p2<z∑S(Ap1p2,P,p2).(2.3)
通过在 (2.3) 上动手脚, 我们就可以构造出上下界.
下界筛
如果我们设 z1≤z 并要求 (2.3) 中的 p2<z1, 则有:
S(A,P,z)≥∣A∣−p1<z∑∣Ap1∣+p1>p2p1<z∑p2<z1∑S(Ap1p2,P,p2),
反复使用该公式就能发现当 z≥z1≥z2≥⋯≥zn 时总有:
S(A,P,z)≥∣A∣−p1<z∑∣Ap1∣+p1>p2p1<z∑p2<z1∑∣Ap1p2∣−p1>p2>p3p1<z∑p2<z1∑p3<z1∑∣Ap1p2p3∣+⋯+p1>p2>⋯>p2np1<z∑p2<z1∑⋯p2n<zn∑∣Ap1p2⋯p2n∣−p1>p2>⋯>p2n>p2n+1p1<z∑p2<z1∑⋯p2n<zn∑p2n+1<zn∑∣Ap1p2⋯p2np2n+1∣=d∣P(z)d∈D−∑μ(d)∣Ad∣,
其中 D− 表示满足下列条件的 d:
d=p1p2⋯pm,p1>p2>⋯>pm,p2k<zk,1≤k≤2m.(2.4)
上界筛
对 (2.3) 再做一次迭代, 就可以发现当 z1≤z 时总有:
S(A,P,z)≤∣A∣−p1<z∑∣Ap1∣+p1>p2p1<z∑p2<z∑∣Ap1p2∣−p1>p2>p3p1<z∑p2<z∑p3<z1∑S(Ap1p2p3,P,p3).
迭代使用这个不等式, 便得:
S(A,P,z)≤∣A∣−p1<z∑∣Ap1∣+p1>p2p1<z∑p2<z∑∣Ap1p2∣−p1>p2>p3p1<z∑p2<z∑p3<z1∑∣Ap1p2p3∣+⋯−p1>p2>⋯>p2n−1p1<z∑p2<z∑⋯p2n+1<zn∑∣Ap1p2⋯p2n+1∣+p1>p2>⋯>p2n>p2n+2p1<z∑p2<z∑⋯p2n+1<zn∑p2n+2<zn∑∣Ap1p2⋯p2np2n+1∣=d∣P(z)d∈D+∑μ(d)∣Ad∣,
其中 d∈D+ 当且仅当:
d=p1p2⋯pm,p1>p2>⋯>pm,p2k+1<zk,1≤k≤2m−1.(2.5)
上下界的统一
为了方便起见, 我们定义 χ0 为集合 D− 的特征函数并让 χ1 为 D+ 的特征函数, 则根据 (2.4) 和 (2.5) 可知当 v=0,1 时 χv(d)=1 当且仅当:
d=p1⋯pm,p1>⋯>pm,p2k+v<zk,1≤k≤min(n,2m−v),(2.6)
因此当我们定义:
Sv(A,P,z)=d∣P(z)∑μ(d)χv(d)∣Ad∣(2.7)
时总有:
S0(A,P,z)≤S(A,P,z)≤S1(A,P,z).(2.8)
至此我们的任务就变成了估计 Sv(A,P,z) 的大小了.
事实上将 (2.8) 的左侧下标更换成任意非负偶数、右侧下标更换成任意正奇数时依然成立.
Sv(A,P,z) 的估计
当 ∣Ad∣ 满足命题 2.1.1 时, 可知存在 −1≤θ≤1 使:
Sv(A,P,z)=XQv(z)d∣P(z)∑μ(d)χv(d)g(d)+θRv(z)d∣P(z)∑χv(d)∣r(d)∣,(2.9)
其中通过精确地处理, Qv(z) 将成为 Sv(A,P,z) 的主项而通过放缩 Rv(z) 会成为 Sv(A,P,z) 的余项.
Qv(z) 的初步处理
从 Qv(z) 的构造来看, 我们期待它最终能近似 V(z), 所以我们不妨先做一次分离, 得:
Qv(z)=V(z)−d∣P(z)∑μ(d)[1−χv(d)]g(d).
结合 (2.8), 可知:
d∣P(z)∑μ(d)[1−χv(d)]g(d)=1≤k≤n∑p1>⋯>p2k+v≥zkzk>p2k+v+1>⋯>pmp2ℓ+v<zℓ,1≤ℓ≤k∑μ(p1⋯pm)g(p1⋯pm)=1≤k≤n∑p1>⋯>p2k+v≥zkp2ℓ+v<zℓ,1≤ℓ≤k∑μ(p1⋯p2k+v)g(p1⋯p2k+v)δ∣P(zk)∑μ(δ)g(δ)=(−1)v1≤k≤n∑V(zk)p1>⋯>p2k+v≥zkp2ℓ+v<zℓ,1≤ℓ≤k∑g(p1⋯p2k+v)=(−1)vV(z)Δv(z),
其中 Δv(z) 满足:
Δv(z)≤1≤k≤n∑V(z)V(zk)(2k+v)!1(zk≤p<z∑g(p))2k+v.
由于当 0≤y<1 时:
y≤y+2y2+3y3+⋯=log(1−y)−1,
所以当 Lk=log[V(zk)/V(z)] 时总有:
Δv(z)≤1≤m≤n∑(2m+v)!eLmLm2m+v.(2.10)
条件 Ω(κ) 与 zm 的设置
为了继续估计 Δv(z), 我们有必要对 g(d) 做更多的要求.
存在常数 A>0 使得 2≤w<z 时总有: w≤p<z∏(1−g(p))−1≤(logwlogz)κ(1+logzA).
因此有:
Lm≤κloglogzmlogz+logzmA.(2.11)
为了让 Lm 最终能够有一个不受制于 z 的上界, Brun 构造了这种形式的 zm:
logzm=e−mΛlogz,1≤m≤n,(2.12)
其中 Λ 是一个固定的正实数. 另一方面我们要求 n 足够大使得 zn≤2<zn−1, 换言之:
e(n−1)Λ<log2logz≤enΛ.(2.13)
Qv(z) 的展开式
将 (2.12) 代入到 (2.11) 中, 便有:
Lm≤mΛκ+logzAemΛ=mΛκ(1+mΛemΛ−1κlogzA)+logzA,
由于
yey−1=m≥1∑m!ym−1
是单调递增函数, 所以根据 (2.13), 可知:
Lm≤mΛκ(1+nΛenΛ−1κlogzA)+logzA<mΛκ(1+nΛe(n−1)ΛκlogzAeΛ)+logzA<mΛκ(1+loglogzA1)+logzA.
这意味着对于所有的 ε>0, 均存在 Z>0 使得当 z>Z 时总有:
Lm<mτΛκ(1+ε)+logzA.(2.14)
现在将 (2.14) 代入到 (2.10), 即得:
Δv(z)<1≤m≤n∑(2m+v)!emτ(mτ+logzA)2m+veA/logz=1≤m≤n∑(2m+v)!emτ(mτ)2m+v(1+mτlogzA)2m+veA/logz<1≤m≤n∑(2m+v)!emτ(mτ)2m+vexp{m2m+vτlogzA+logzA}≤1≤m≤n∑(2m+v)!emτ(mτ)2m+v{1+O(logz1)}.
我们不难看出当 rm 表示上述和式的第 m+1 项除以第 m 项的商时, 有:
rm∼4m2eτ{mτ(m+1)τ}2m[(m+1)τ]2∼4τ2e2+τ,
因此当 41τ2e2+τ<1 时可知当 z 充分大时便有:
Δv(z)<1≤m≤n∑(2m+v)!emτ(mτ)2m+v+O(logz1)<m≥1∑(2m+v)!emτ(mτ)2m+v.
综上所述我们就得到了结论:
对于每个充分大的 z, 总存在 −1≤θ≤1 使得: V(z)Qv(z)=1+θm≥1∑(2m+v)!emτ(mτ)2m+v,其中 v 为非负整数, 正实数 τ 满足: 0<4τ2e2+τ<1.
接下来就只需要估计 Rv(z) 的上界了.
Rv(z) 的上界估计
为了灵活性, 我们并不打算像 g(d) 一样给 r(d) 增添限制条件. 所以在本节中我们只会估计这个和式:
Mv(z;f)=d∣P(z)∑χv(d)f(d),(2.15)
其中 f(d) 为一个积性函数. 因此根据 (2.6), 可知 v≥0 时总有:
Mv(z;f)≤(1+p<z∑f(p))1+v1≤m<n∏(1+p<zm∑f(p))2.(2.16)
在很多应用场景中 f(p) 通常是有界的, 所以此时我们可以追加条件, 即存在 B>0 使 y≥2 时总有:
1+p<y∑f(p)≤logyBy.
将这个条件代入到 (2.16) 中, 便可发现 z 充分大时:
Mv(z;f)<z1+v1≤m<n∏(logzmBzm)2<z1+v+2(e−Λ+e−2Λ+⋯)1≤m<n∏(logzBemΛ)2=z1+v+eΛ−12en(n−1)Λ1≤m<n∏(logzB)2=z1+v+eΛ−12(logzBenΛ/2)2(n−1).
现在根据 (2.13), 可知:
logzenΛ/2<logzlog2eΛ/2,
另一方面根据 (2.14) 我们知道存在 ε′>0 使得:
eΛ−12=eτ/κ−12+ε′.
将这些结论结合起来, 我们就得到了用于估计 Rv(z) 的辅助工具:
倘若积性函数 f(d) 满足: p<z∑f(p)=O(logzz),则对于所有的 ε>0, 当 z 充分大时总有: d∣P(z)∑χv(d)f(d)<z1+v+eτ/κ−12+ε.
结论
现将引理 2.3.2 和 (2.9) 相结合, 就有:
若特征函数 χv(d)=1 当且仅当 d 满足 (2.6), 且Rv(z)=d∣P(z)∑χv(d)∣r(d)∣,ηv(τ)=m≥1∑(2m+v)!emτ(mτ)2m+v,则当 A 满足命题 2.1.1 和 Ω(κ), 正数 τ 满足 0<τ2e2+τ<4, v1 为正奇数, v2 为非负偶数时S(A,P,z)<XV(z)[1+ηv1(τ)]+Rv2(z),S(A,P,z)>XV(z)[1−ηv2(τ)]−Rv2(z)对所有充分大的 z 均成立.