50. Pow(x, n) - x * x 和 x **2 的区别

求问大家 x * x 和 x ** 2 的区别。
这个题如果我在第20行（代码如下）改成 cur_prod ** 2, 会给我报numerical out of range的错。不知道这是为什么？
代码：
class Solution(object):
    def myPow(self, x, n):
        """
        :type x: float
        :type n: int
        :rtype: float
        """   
        if n == 0:
            return 1
        if n < 0:
            return 1/self.myPow(x, -n)

        total_pow = 0
        total_prod = 1
        while (total_pow < n):
            cur_pow = 1
            cur_prod = x
            while (2 * cur_pow < n - total_pow):
                cur_pow += cur_pow
                cur_prod = cur_prod * cur_prod
            total_pow +=  cur_pow
            total_prod *= cur_prod

        return total_prod

楼主问的问题可以简化为下面两行代码：

1. # 输出：inf
   
2. print 2e+154 * 2e+154
   
3. 
   
4. # 34号（溢出）错误
   
5. print 2e+154 ** 2
   

复制代码

下面是原理部分：
--
(1) ** 的Python实现在这里（不同Python版本会有变动，但总体上是差不多的）：https://github.com/python/cpytho ... /mathmodule.c#L2561
用 2e+154 ** 2 作为例子，那么上面链接的函数里会调用C语言的pow()函数先计算 2e+154 ** 2 得到一个double类型的结果，其值为infinity —— 这是IEEE标准下一个合法的浮点值：https://en.wikipedia.org/wiki/IE ... d_negative_infinity

然后在第2616行有一个Py_IS_INFINITY检查，发现double类型结果为infinity时就会设置错误代码为 ERANGE，也就是34号错误。

--
(2) 对于乘法操作，代码在这里：https://github.com/python/cpytho ... /floatobject.c#L568
这里并没有检查infinity而是直接返回乘法结果了

大哥爱挖挖

盲猜是因为double精度炸了……

Jarvet 发表于 2019-6-6 07:19
盲猜是因为double精度炸了……

为何俩 x相乘就木有关系 = =

magicsets 发表于 2019-6-6 12:06
楼主问的问题可以简化为下面两行代码：
[mw_shl_code=python,true]# 输出：inf
print 2e+154 * 2e+154

哇！满分！ 非常感谢！

狂暴CNM地 发表于 2020-12-9 14:22
大哥爱挖挖

流口水大哥爱挖挖