追求完整的无理数：挑战手机计算器的极限

本文来自微信公众号：返朴 （ID：fanpu2019），作者：姬扬

撰文 | 姬扬（中国科学院半导体研究所）

我们都知道，无理数（irrational number）“不讲理”。按整数标准，它们大多身形残缺，但是它们中间有很多数很有理想，总是想追求完整：有的是靠自己，有的是靠同志。

比如说，是无理数，可能也是无理数，但是就是整数了，而且它很“2”。

再比如说，自然常数 e 和圆周率 π 都是无理数，但是只需要一个虚数的帮助，它们就可以变成整数：e^iπ=-1；结果是个负整数，它还不够“2”，否则就会有 e^2iπ=1。

当然，不是所有的努力都能够功德完满。自然常数 e、圆周率 π 和 163 的平方根（），它们也努力过追求完整，可惜就差了那么一点点，很小的一点点。

多小的一点点呢？

它跟整数的差别在于小数点后第 13 位，也就是说，小数点后经过了 12 个 9，才变成了 2。

为什么这样呢？这里面有一套很复杂的数学理论，我也看不懂，但是可以引用一个公式：

它是否正确，我们可以用计算器检验一下。

手机的计算器很强大，特别是它的科学计算器模式。例如它可以轻松地把圆周率 π 的值算到小数点后一百位。

为了便于大家亲自检验，我会介绍详细的操作过程。首先，打开手机自带的计算器（我的手机是华为 P30），它有两个模式 —— 竖版的“标准计算器”和横版的“科学计算器”，可以任意地选择切换。

进入“科学计算器”，就自动切换到横版，上方是两行的显示屏，下方是键盘。按“π”键，π 就出现在显示屏的第一行，同时在第二行显示出它的数值，“3.1415...1971”。这一行可以显示 42 个数字或字符。通过触屏操作，我们可以对数字进行复制粘贴操作。

现在，我们来验证手机自带的 π 究竟算到了多少位。用 π 减去刚才得到的 42 位数字，即“π-3.1415…1971”。

继续同样的操作，可以得到更多的位数。

从菜单里选择“历史纪录”，可以找到以前做的计算，比如说我刚做的把 π 算到 100 位。

算得对不对呢，我从网上找了 π 值小数点后 200 位，如下：

π= 3.1415926535 8979323846 2643383279 5028841971 6939937510 5820974944 5923078164 0628620899 8628034825 3421170679 8214808651 3282306647 0938446095 5058223172 5359408128 4811174502 8410270193 8521105559 6446229489 5493038196 。

跟上面的图片对照一下就可以看到，手机计算器的结果有 115 位小数，全部正确。使用这一方法显然还可以算更多位，但是究竟能算多少位，我没有试验，毕竟用手机操作复制和粘贴太麻烦了。

值得一提的是，这个表现比我的 PC 电脑系统附带的计算器厉害。利用类似的操作，PC 系统的计算器可以得到“π= 3.14...... 6939931148 ”（第 41-50 位应该是 6939937510），从小数点后第 47 位开始就错了。

网上的 π 值就一定正确吗？它可能是正确的，也有可能手机计算器特地保留了 π 的很多位数，也许藏在哪里备用呢。

接下来，让我们算算 —— 想必这个数是不会预存的。结果发现，两个计算器的表现还是不一样。仔细观察可以看出，两部计算器的第一次计算结果应是相同的：计算机的结果最后的 4 位数字是 9871；而手机对应的位置是 9870，其后面还有 9 个数字“892347382”。当用电脑计算器做减法后，得到的结果是“-1.076524...”，前面还有一个负号。考虑到负号，这两个计算结果一直到我们手机显示的倒数第三位，也就是 3（手机计算）和 4（电脑计算），这里确实有差别。我觉得，相比于电脑，我们还是相信手机好了，毕竟它给出的位数更多，而且计算的 π 值也跟网上的结果一样。

电脑计算器计算-12.767…781（的 30 位小数近似值）的结果。

手机计算器计算的结果，有 39 位小数。

最后，我们算一算。先把前面的公式再复习一遍：

我用手机自带的科学计算器做了几个不同的计算，分别是直接计算的值（用前文的方法）和用公式近似计算，结果如下：

①直接计算，得到 18 位整数和 23 位小数：最开始的 12 个小数都是 9，第 13 个是 2。

②计算公式的一级近似，也就是第一项的贡献。

③计算公式的二级近似，也就是前两项的贡献。

④计算公式的三级近似，也就是前三项的贡献。

⑤第二次直接计算，这次算的是减去 18 位整数和小数点后的 12 个 9。

⑥第三次直接计算，这次算的是减去第一次直接计算得到的结果，也就是 18 位整数和小数点后的 23 位小数。

⑦第四次直接计算，减去第一次直接计算和第三次直接计算的结果，得到了小数点以后的 99 位小数。

简单分析可以得出：直接计算（①）最简单也最粗糙，但是已经得到了 e、π 和这三个无理数合作追求完整的全过程；公式没有告诉整数是多少（这个很容易算的），但是给出了更多更精确的小数位（②、③、④）；计算器比这个公式的前三项还要强大（⑤和⑥），甚至能够计算到小数点后 99 位（⑦）。用两种计算的结果（用公式计算和直接计算）符合得很好，因此我们相信手机计算的结果是正确的。

我也用电脑的计算器算过，虽然也能见证它们仨追求完整的过程，但是在小数点后几十的地方，仍然给出了不同的结果。

最后再谈谈这个公式。非常接近整数这件事，在一些讲计算数学的书里提到过，有时候用作例子说明计算的精度和稳定性。然而，关于它为什么这么接近整数，讲得就很少了。而且基本上都是先说，这个道理给普通人是讲不明白的，因为要涉及椭圆曲线、模形式以及其他一些理论。

我找到了一篇介绍性的文章 [1]，看得一头雾水，但也不是没有收获。除了找到前面用的那个公式以外，我还发现 —— 可能大家也都知道 —— 数学家的数学虽然很好，但是他们算术不一定强，经常会搞错个符号之类的。正如文中所示，作者在欢呼成功之前的最后两行公式中“196884”前面的符号不一样，肯定有一个是错的 —— 而且我肯定，是最后那行里的错了。

“大行不拘细谨，大礼不辞小让。”连数学家都不免犯错，三个努力追求完整的无理数最后没有取得特别完满的结果，也就不那么令人意外了。

的各种近似值，最后一个有 99 位小数。

参考文献

• [1] https://www.isibang.ac.in/~sury/episq163.pdf

广告声明：文内含有的对外跳转链接（包括不限于超链接、二维码、口令等形式），用于传递更多信息，节省甄选时间，结果仅供参考，IT之家所有文章均包含本声明。