如何不使用Math库实现log10(x)？附代码与公式运用疑问

核心思路 #

log₁₀(x)的本质是找到一个实数y，使得10ʸ = x。由于10ʸ是单调递增函数，我们可以用二分法在合理区间内逼近y的值，直到达到设定的精度要求。同时需要先实现一个支持小数指数的幂函数（原findPower仅支持整数指数），因为二分过程中需要计算10的小数次幂。

代码实现 #

class Solution(object):
    def findPower(self, number, exponent):
        # 处理整数次幂，保留原函数逻辑
        resultPower = 1
        exp = int(exponent)
        for i in range(1, exp + 1):
            resultPower *= number
        return resultPower

    def power(self, a, b):
        # 计算a^b，支持实数b
        if b == 0:
            return 1.0
        # 处理负数指数：a^(-b) = 1/(a^b)
        if b < 0:
            return 1.0 / self.power(a, -b)
        
        # 拆分指数为整数部分和小数部分
        integer_part = int(b)
        fraction_part = b - integer_part
        
        # 整数部分用原findPower计算
        integer_power = self.findPower(a, integer_part)
        
        # 小数部分用泰勒展开计算10^f（f∈[0,1)）
        # 预先定义ln10的近似值（不依赖数学函数，硬编码常数）
        ln10 = 2.302585093
        x = fraction_part * ln10
        # 泰勒展开e^x ≈ 1 + x + x²/2! + x³/3! + ... ，取前10项足够精度
        exp_x = 1.0
        term = 1.0
        for n in range(1, 10):
            term *= x / n
            exp_x += term
        fraction_power = exp_x
        
        return integer_power * fraction_power

    def log10(self, x):
        # 边界条件处理
        if x <= 0:
            raise ValueError("log10的定义域是x>0")
        if x == 1:
            return 0.0
        if x == 10:
            return 1.0
        
        epsilon = 1e-6  # 精度控制，可调整
        low, high = 0.0, 1.0
        
        # 确定二分初始区间
        if x > 1:
            while self.power(10, high) < x:
                high *= 2
        else:  # x < 1
            low = -1.0
            while self.power(10, low) > x:
                low *= 2
        
        # 二分迭代逼近
        while high - low > epsilon:
            mid = (low + high) / 2
            mid_power = self.power(10, mid)
            if mid_power < x:
                low = mid
            else:
                high = mid
        
        return (low + high) / 2

    def mySqrt(self, x):
        # 基于log10计算平方根：sqrt(x) = 10^(log10(x)/2)
        if x < 0:
            raise ValueError("负数没有实数平方根")
        if x == 0:
            return 0.0
        return self.power(10, self.log10(x) / 2)

关键细节说明 #

1. 幂函数扩展：新增power函数，将指数拆分为整数和小数部分。整数部分复用原findPower，小数部分通过泰勒展开计算e^{(f*ln10)（因为10}f = e^(f*ln10)），硬编码ln10的近似值避免依赖数学函数。
2. 二分区间确定：根据x与1的大小关系，动态扩展初始区间，确保10^low ≤ x ≤10^high。
3. 精度控制：通过epsilon控制迭代终止条件，平衡计算精度和效率。
4. 边界处理：提前处理x=1、x=10等特殊值，以及x≤0的非法输入。

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