科学空间|Scientific Spaces

BoJone很喜欢Python，也很喜欢数论，所以就喜欢利用Python玩数论了。平时也喜欢自己动手写一些数论函数，毕竟Python支持大整数高精度运算，这点是非常好的；但是，在很多实际应用中，还是希望能有一个现成的数论函数库来调用。之前尝试过数学研发网的HugeCalc库，但是由于各种不熟悉不了了之。后来论坛上的无心老兄推荐了PARI/GP，小试一下，居然在Python上成功调用了。以后再也不用担心Python上的数论计算问题了，呵呵~

点击阅读全文...

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

素数之积

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

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

点击阅读全文...

有了上一篇文章的$\prod\limits_{p\leq n}p < 4^{n-1}$的基础，我们其实已经很接近Bertrand假设的证明了。Bertrand假设的证明基于对二项式系数$C_n^{2n}$的素因子次数的细致考察，而在本篇文章中，我们先得到一个关于素数之积的下限公式，然后由此证明一个比Bertrand假设稍微弱一点的假设。最后，则通过一个简单的技巧，将我们的证明推动至Bertrand假设。

二项式系数的素因子

首先，我们考察$n!$中的素因子$p$的次数，结果是被称为Legendre定理的公式：

$n$中素因子$p$的次数恰好为$\sum\limits_{k\geq 1}\left\lfloor\frac{n}{p^k}\right\rfloor$。

证明很简单，因为$n!=1\times 2\times 3\times 4\times \dots \times n$，每隔$p$就有一个$p$的倍数，每隔$p^2$就有一个$p^2$的倍数，每隔$p^3$就有一个$p^3$的倍数，每增加一次幂，将多贡献一个$p$因子，所以把每个间隔数叠加即可。注意该和虽然写成无穷形式，但是非零项是有限的。

点击阅读全文...

在小学的时候，数学老师就教我们除法运算：

被除数 = 除数 × 商 + 余数

其中，余数要小于除数。不过，我们也许未曾想到过，这一运算的成立，几乎是自然数$\mathbb{N}$所有算术（数论）运算性质成立的基础！在代数中，上面的运算等式称为带余除法（division algorithm）。如果在一个整环中成立带余除法，那么该整环几乎就拥有了所有理想的性质，比如唯一分解性，也就是我们说的算术基本定理。这样的一个整环，被称为唯一分解整环（Unique factorization domain）。

欧几里得整环

Euklid-von-Alexandria_1

唯一分解定理说的是在一个整环之中，所有的元素都可以分解为该整环的某些“素元素”之积，并且在不考虑元素相乘的顺序和相差单位数的意义之下，分解形式是唯一的。我们通常说的自然数就成立唯一分解定理，比如$60=2^2\times 3\times 5$，这种分解是唯一的，这看起来相当显然，但实际上唯一分解定理相当不显然。首先，并不是所有的整数环都成立唯一分解定理的，我们考虑所有偶数组成的环$2\mathbb{Z}$，要注意，在$2\mathbb{Z}$中，2、6、10、30都是素数，因为它们无法分解成两个偶数的乘积了，但是$60=6\times 10=2\times 30$，存在两种不同的分解，因此在这样的数环中，唯一分解定理就不成立了。

点击阅读全文...

现在可以开始$n=3$的证明了。在实整数范围内n=3的证明看起来相当复杂，而且跟n=4的证明似乎没有相通之处。然而，如果我们在$\mathbb{Z}[\omega]$中考虑$x^3+y^3+z^3=0$无解的证明，就会跟n=4时有很多类似的地方，而且事实上证明比n=4时简单（要注意在实整数范围内的证明，n=4比n=3简单。费马完成了n=4的证明，但是没完成n=3的证明。）。我想，正是这样的类似之处，才让当初还没有完成证明的数学家拉梅就自信他从这条路可以完成费马大定理的证明。（不过，这自信却是失败的案例：拉梅的路不能完全走通，而沿着这条路走得更远的当属库默，但即便这样，库默也没有证明费马大定理。）

证明跟$n=4$的第二个证明是类似的。我们先往方程中添加一个单位数，然后证明无论单位数是什么，方程在$\mathbb{Z}[\omega]$中都无解。这是一个很妙的技巧，让我们证明了更多的方程无解，但是却用到了更少的步骤。事实上，存在着只证明$x^3+y^3+z^3=0$无解的证明，但需要非常仔细地分析里边的单位数情况，这是相当麻烦的。本证明是我参考了Fermats last theorem blogspot上的证明，然后结合本系列n=4的第二个证明，简化而来，主要是减少了对单位数的仔细分析。

点击阅读全文...

是时候了！

前面用好几篇文章为费马大定理的证明铺设了道路，当然，相当于完整的费马大定理证明来说，这几篇文章只不过是沧海一粟而已。不过，它们已经足够用来完成费马大定理在n=4时的证明了。我们很快会看到高斯整数为n=4所带来的简洁的证明，而这让我们坚信，这一道路可以走得更远。

不定方程$x^4+y^4=z^2$在$\mathbb{Z}[i]$中没有全不为0的解。

点击阅读全文...

在上一篇文章中，笔者提到似乎证明n=4时必须要证明$x^4+y^4=z^2$无解而不能只证明$x^4+y^4=z^4$无解。不过，在今天中午研究的时候，笔者发现了另外一个n=4的证明，它同样是在$\mathbb{Z}[i]$中，但是，证明的则是指数全是4的形式，但是，又不单单是$x^4+y^4=z^4$的形式，而是$\varepsilon x^4+y^4=z^4$，$\varepsilon$是单位数。这个证明过程，我觉得应该更接近n等于其他奇素数时的证明，遂补充了这篇文章，供大家参考。读者可以对比着上一篇文章进行比较阅读。

引理

用$\varepsilon_1,\varepsilon_2,\varepsilon_3,\varepsilon$表示$\mathbb{Z}[i]$中的单位数，下面先证明

如果方程$\varepsilon_1 x'^4 +\varepsilon_2 y'^4+\varepsilon_3 z'^4=0$在$\mathbb{Z}[i]$中有全不为0的解，那么在经过适当的化简和整理之后，方程必有形式$\varepsilon x^4+y^4=z^4$，其中$(x,y,z)$是$(x',y',z')$的某个置换，$\xi^2|x$。

点击阅读全文...

Gotthold_Eisenstein

是时候向n=3进军了，为了证明这个情况，我们需要一个新的数环：艾森斯坦整数（Eisenstein Integer）。艾森斯坦是德国著名数学家，同时代的高斯曾经评价：“只有三个划时代的数学家：阿基米德，牛顿和艾森斯坦。”足见艾森斯坦的成就斐然。事实上，阅读费马大定理的研究史，同时也是在阅读数学名人录——没有超高的数学，几乎不可能在费马大定理中有所建树。

基本定义

跟高斯整数一样，艾森斯坦整数也是复整数的一种，其中，高斯整数是以1和$i$为基，$i$其实是一个四次单位根，也就是$x^4-1=0$的一个非实数根，因此高斯整数也叫做四次分圆整数；而艾森斯坦整数以1和$\omega$为基，$\omega$是三次单位根，也就是$x^3-1=0$的一个非实数根。任意一个艾森斯坦整数都可以记为$a+b\omega,\,a,b\in\mathbb{Z}$，艾森斯坦整数环记为$\mathbb{Z}[\omega]$，也称为三次分圆整数环。

点击阅读全文...