为何Python中(1+(1/x))**x与math.e比较会出现真假振荡？

先看你给出的代码：

import math
e = lambda x: (1+(1/x))**x
for x in range(1, 20):
    print(math.e > e(10**x))

数学上的结论没错：$(1+\frac{1}{x})^x$确实是单调递增且极限为自然常数e的序列，但这是实数域精确计算的结果。代码里出现振荡，完全是计算机浮点数运算的问题，核心原因有两个：

• 浮点数精度有限：Python用的是64位双精度浮点数，只有约15-17位有效数字。当x是10的高次幂时，1/x的数值极小，小到和1相加后，浮点数无法区分出和1的差异——比如x=10^18时，1/x=1e-18，这个值远小于1附近浮点数的最小精度单位（ULP，约2.2e-16），此时1 + 1/x会被直接舍入成1，计算结果就是1^x=1，肯定小于math.e，输出True。而当x的规模刚好让1+1/x能被浮点数区分开时，计算出的结果又可能因为误差偏大，输出False。

• 幂运算的舍入误差：就算1+1/x能被浮点数正确表示，计算它的x次方时，浮点数幂运算的过程也会引入舍入误差。不同大小的x，误差的方向可能随机变化——有时候误差让结果偏小，有时候偏大，这就导致和math.e比较时出现True和False交替的情况。

简单说，就是数学上的完美序列，到了计算机里因为浮点数的近似计算，变成了带有随机误差的波动序列，自然会出现和e比较时的振荡现象。

内容的提问来源于stack exchange，提问作者Thiago Martins