经典分治算法——Karatsuba乘法算法

勤劳的进取家

已于 2025-03-31 13:32:11 修改

阅读量525

点赞数 5

CC 4.0 BY-SA版权

分类专栏：算法设计文章标签：算法数据结构概率论线性代数线性回归数学建模

于 2025-03-31 13:31:44 首次发布

28 篇文章

订阅专栏

基本思想：

给定两个 n 位整数 x 和 y，将它们各自拆分成两部分：

$x = a \times 10^{m} + b,\quad y = c \times 10^{m} + d,$

其中 $m = \lceil n/2 \rceil$ ，a, b, c, d分别表示整数 x 和 y 的高位和低位部分。

按照普通乘法的公式：

$x \times y = a \times c \times 10^{2m} + (a \times d + b \times c) \times 10^{m} + b \times d.$

传统方法需要计算 4 个大数乘法：a×c 、a×d 、b×c 和 b×d 。

Karatsuba的优化思路在于减少乘法次数。注意到：

$(a+b) \times (c+d) = a \times c + a \times d + b \times c + b \times d,$

因此可以推出：

$a \times d + b \times c = (a+b) \times (c+d) - a \times c - b \times d$

这样只需计算三个乘法：

然后利用上述关系，计算中间项，从而减少乘法运算的次数。

分治递归关系的推导：

设 T(n) 表示两个 n 位数相乘的时间复杂度。按照 Karatsuba 算法的步骤：

因此递归关系为：

T(n)=3T(n2)+O(n)

利用主定理分析：

对于递归关系 $T(n) = a\, T\left(\frac{n}{b}\right) + O(n^d)$ ，这里 a=3 ，b=2 ，d=1 。

计算临界指数：

$\log_b a = \log_2 3 \approx 1.585.$

因为 $\log_2 3 > 1$ （即 d），所以根据主定理：

$T(n) = O\left(n^{\log_2 3}\right) \approx O\left(n^{1.585}\right)$

这比传统的 O(n²) 算法更高效，尤其在大整数乘法中优势明显。

基本情况：
- 如果 n 足够小（例如 n = 1 或 n低于某个阈值），直接使用常规乘法返回结果。
分解数字：
- 设 m = ⌈n/2⌉，将整数 x 和 y 分解为：
  $x = a \times 10^{m} + b,\quad y = c \times 10^{m} + d.$
递归计算三个乘积：
- 计算 $ac = \text{Karatsuba}(a, c)$ 。
- 计算 $bd = \text{Karatsuba}(b, d)$ 。
- 计算 $(a+b)(c+d) = \text{Karatsuba}(a+b, c+d)$ 。
合并中间结果：
- 计算中间项：
  $ad+bc = (a+b)(c+d) - ac - bd.$
- 计算最终结果：
  $x \times y = ac \times 10^{2m} + (ad+bc) \times 10^{m} + bd.$
返回结果：
- 返回合并后的结果。

算法思路： 利用分治思想，将大整数分解为高低位两部分，通过三个子问题（而不是四个）递归计算乘积，利用加减法合并结果。
理论推导： 递归关系为 T(n)=3T(n/2)+O(n) ，由主定理得出时间复杂度为 $O(n^{\log_2 3}) \approx O(n^{1.585})$ 。
算法步骤： 包括基本情况处理、分解数字、递归计算、合并结果以及返回最终乘积。