如何在不使用任何库的情况下用Python实现对数运算？

当然可以不依赖任何库实现对数运算，你当前的遍历方法存在明显的效率和功能缺陷：

• 仅能处理底数的整数次幂等于输入值的特殊情况，无法计算非整数结果（比如log₁₀(2)≈0.3010，你的代码根本返回不了值）
• 遍历范围0到base完全不合理，当底数很大时，遍历次数会爆炸式增长，导致速度极慢

高效实现方案 #

我们可以利用对数的数学性质：log_b(a) = ln(a)/ln(b)，先实现自然对数（ln），再通过转换得到任意底数的对数。下面是基于级数展开的高效实现，收敛速度快且支持任意正实数输入：

实现自然对数（ln）

def ln(x):
    if x <= 0:
        raise ValueError("输入必须大于0")
    
    # 将x缩放至(0.5, 2)区间，利用ln(x) = ln(x/2^k) + k*ln2（ln2≈0.69314718056）
    k = 0
    while x >= 2:
        x /= 2
        k += 1
    while x < 0.5:
        x *= 2
        k -= 1
    
    # 令t=(x-1)/(x+1)，此时t∈(-1/3,1/3)，用级数ln(x)=2*(t + t³/3 + t⁵/5 + ...)快速收敛
    t = (x - 1) / (x + 1)
    result = 0.0
    t_power = t
    denominator = 1
    
    # 迭代直到精度满足（这里取1e-10）
    while abs(t_power / denominator) > 1e-10:
        result += t_power / denominator
        t_power *= t * t
        denominator += 2
    
    result *= 2
    result += k * 0.69314718056
    return result

实现任意底数的对数

def log(a, base):
    if a <= 0 or base <= 0 or base == 1:
        raise ValueError("a必须大于0，base必须大于0且不等于1")
    return ln(a) / ln(base)

测试示例 #

# 计算log₁₀(2)，预期结果≈0.3010
print(log(2, 10))  # 输出约0.30102999566

这个方案的优势：

• 支持所有合法的正实数输入，能返回高精度的非整数结果
• 收敛速度快，仅需几十次迭代就能达到1e-10级别的精度
• 时间复杂度与输入大小无关，仅由精度要求决定

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