素数之美1：所有素数之积

在之前的欧拉数学中，我们计算过所有素数的倒数之和，得出素数的倒数之和是发散的，从而这也是一个关于素数个数为无穷的证明。在本篇文章中，我们尝试计算所有素数之积，通过一个简单的技巧，得到素数之积的一个上限（以后我们也会计算下限），从而也得到$\pi(n)$的一个上限公式。更重要的，该估计是初等地证明Bertrand假设（说的是n与2n之间定有一个素数）的重要基础之一。本文内容部分参考自《数学天书中的证明》和《解析和概率数论导引》。

素数之积

笔者已经说过，数论的神奇之处就是它总是出人意料地把数学的不同领域联系了起来。读者很快就可以看到，本文的证明和组合数学有重要联系（但仅仅是简单的联系）。关于素数之积，我们有以下结论：

不超过$n$的所有素数之积小于$4^{n-1}$。

如果用$p$来表示素数，那么可以记为$\prod\limits_{p\leq n}p < 4^{n-1}$。事实上，由更深入的技巧可以大大改进该估计，因为由素数定理知$\pi(n)\sim \frac{n}{\ln n}$，从而$\prod\limits_{p\leq n}p \sim e^n$，也就是说，对于任意$\varepsilon > 0$，存在$N(\varepsilon) > 0$，使得对于任意的$n > N(\varepsilon)$，都有$\prod\limits_{p\leq n}p < (e+\varepsilon)^n$。但是，本文所采用的方法和估计都是比较简单的，所以只能够得到上限$4^n$，但之后读者就可以发现，对于证明Betrand假设来说，这个估计足够了。

下面我们开始证明它。可以检验，对于n=2该结论是成立的。而一个很显然的事实是，如果该估计对于奇数n都成立，那么它对偶数n+1必然成立。所以我们只证明$n=2m+1$时成立。我们采用数学归纳法，假设$n=1,2,\dots,2m$时，该结论都成立，于是我们将$n=2m+1$时分拆成两部分
$$\prod_{p\leq 2m+1}p=\prod_{p\leq m+1}p\prod_{m+1 < p\leq 2m+1}p < 4^m\prod_{m+1 < p\leq 2m+1}p$$
最后一个不等号是由归纳假设得到的。现在还差一部分，即m+1到2m+1之间的素数之积，我们考虑二项式系数
$$C_{m+1}^{2m+1}=\binom{2m+1}{m+1}=\frac{(2m+1)!}{(m+1)! m!}$$
由定义知$C_{m+1}^{2m+1}$是一个整数，并且它能够被大于m+1、小于等于2m+1的所有素数整除，所以显然
$$\prod_{m+1 < p\leq 2m+1}p\leq C_{m+1}^{2m+1}$$
而对$C_{m+1}^{2m+1}$进行估计就会发现
$$C_{m+1}^{2m+1} \leq 2^{2m}=4^m$$
这是因为$C_{m+1}^{2m+1}=C_{m}^{2m+1}$，且
$$C_{m+1}^{2m+1}+C_{m}^{2m+1}\leq \sum_{k=0}^{2m+1}C_{k}^{2m+1}=2^{2m+1}$$
指数4是最优的，如果不改变估计的方法，从该二项式系数中无法得到更好的指数。为了认识到这一点，只需要用String公式对$C_{m+1}^{2m+1}$进行渐进近似计算。

有了对$C_{m+1}^{2m+1}$的估计，那么就有
$$\prod_{p\leq 2m+1}p< 4^m\prod_{m+1 < p\leq 2m+1}p\leq 4^{2m}$$
也就是说，结论对于$n=2m+1$成立。从而由归纳假设知，该结论对于所有的正整数都成立。

素数个数

有了这个估计，我们还可以给出$\pi(n)$的一个上限，只需要注意到
$$t^{\pi(n)-\pi(t)} < \prod_{p\leq n}p < 4^{n-1}$$
$t$是任意常数，两边取对数得
$$\pi(n) < \frac{(n-1)\ln 4}{\ln t}+\pi(t)\leq \frac{n\ln 4}{\ln t}+t$$
在数论中，为了改进估算结果，通常是这样的思路，引入一些未知参数，比如上面的t，然后通过调整未知参数来得到最优估计。比如上面的不等式，我们设$f(t)=\frac{(n-1)\ln 4}{\ln t}+t$，求它的最小值：
$$f'(t)=-\frac{n\ln 4}{(\ln t)^2 t}+1$$
令导数为0，得到$n\ln 4=(\ln t)^2 t$，这是一道超越方程，没有办法得到显式精确解，通过迭代可以给出近似解，改成
$$t_{k+1}=\frac{n\ln 4}{(\ln t_k)^2}$$
取$t_0=n$，迭代一次得到$t_1=\frac{n\ln 4}{(\ln n)^2}$，为了不得到复杂的解，我们的迭代到此结束，并且简单起见，省去因子$\ln 4$，即简单地取$t=n/(\ln n)^2$，代入并整理得到
$$\pi(n) < \left(\ln 4+\frac{8\ln\ln n}{\ln n}\right)\frac{n}{\ln n}$$
该估计弱于素数定理很多，但它是基于最简单的计数的结果，足可满意。

转载到请包括本文地址：https://kexue.fm/archives/2789

更详细的转载事宜请参考：《科学空间FAQ》