Python中如何实现超大数幂运算？常规pow方法无法得到结果

首先得说清楚：你直接用pow(a, b)得不到结果，根本原因不是Python算不了，而是这个结果大到离谱——它的位数大概是5.8×10^66位，这个量级的数字，别说你的电脑内存，就算把整个宇宙的原子都拿来当存储单元，都存不下完整结果。所以要解决问题，得先明确你到底需要这个幂运算的什么结果：

情况1：只需要幂运算的模结果（最常见场景） #

如果你是要算a^b mod m（比如密码学、数论里的需求），Python内置的pow其实已经支持三参数形式，而且是高度优化的快速幂实现，哪怕a、b、m都是超大数，速度也极快：

a = 62608558862573792084872798679396455703616395237802859621162736207631538899993
b = 93910650126758265671774994856253142403789359314618444886584691522424141933664
mod = 10**100  # 示例模数，换成你需要的
c = pow(a, b, mod)
print(c)

这个方法不需要计算完整的大整数，每一步运算都会取模，所以内存占用极小，效率拉满。

情况2：只需要结果的前几位或后几位 #

• 后n位：本质就是模10^n，用上面的三参数pow就行，比如要后20位，就用pow(a, b, 10**20)。
• 前n位：可以用对数计算，原理是log10(a^b) = b*log10(a)，这个值的整数部分是结果的位数减1，小数部分可以用来算出前几位：

  from decimal import Decimal, getcontext
  # 设置足够的精度，避免误差，这里设为比b的位数多10位
  getcontext().prec = len(str(b)) + 10
  # 计算log10(a)
  log_a = Decimal(a).ln() / Decimal(10).ln()
  # 计算log10(a^b)
  log_total = Decimal(b) * log_a
  # 取出小数部分
  fractional = log_total - log_total.to_integral_value(rounding="ROUND_FLOOR")
  # 取前5位的示例
  first_5_digits = int((Decimal(10)**fractional) * Decimal(10**4))
  print(first_5_digits)
  

  这里用decimal模块是因为普通的math.log10精度不够，处理超大数会有误差。

情况3：真的需要完整结果？ #

别想了，物理上不可能。这个结果的字节数大概是7×10^{65，而整个宇宙的原子总数才约1×10}80，就算每个原子存1字节，都得用掉宇宙7%的原子才能存下——这显然不现实。

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