Python中如何获取无小数的精确整数立方根

问题原因 #

这是典型的二进制浮点数精度损失问题：1/3 无法用二进制浮点数精确表示，存储时会截断为近似值，做幂运算后就会出现微小的误差。8、27的计算误差刚好落在浮点数的进位阈值内，所以结果显示正常，而64的计算结果刚好差一点点到4，就会输出3.9999999999999996。

可用解决方法 #

• 方法1：使用Python标准库math的内置立方根方法（Python 3.11及以上版本支持）
  math.cbrt是专门为立方根计算优化的实现，避免了通用幂运算的精度问题：

import math
print(math.cbrt(64))
# 输出结果：4.0

• 方法2：对结果做精度舍入（适配所有Python版本）
  如果确认输入是完全立方数，可以对浮点数结果做合理位数的舍入，消掉误差位：

raw_res = 64 ** (1/3)
correct_res = round(raw_res, 9)
print(correct_res)
# 输出结果：4.0

如果需要严格处理整数完全立方根，也可以用二分法实现整数运算的立方根计算，完全避开浮点数误差：

def integer_cbrt(n: int) -> int:
    if n < 0:
        return -integer_cbrt(-n)
    low, high = 0, n
    while low <= high:
        mid = (low + high) // 2
        mid_cube = mid ** 3
        if mid_cube == n:
            return mid
        elif mid_cube < n:
            low = mid + 1
        else:
            high = mid - 1
    # 非完全立方数可根据需求返回近似值，这里示例返回最接近的整数
    return round(n ** (1/3))

print(integer_cbrt(64))
# 输出结果：4

• 方法3：使用数值计算库的立方根实现
  如果项目中已经用到numpy，可以直接调用numpy.cbrt，精度也经过优化：

import numpy as np
print(np.cbrt(64))
# 输出结果：4.0

内容的提问来源于stack exchange，提问作者Nikola。