Karatsuba算法主要应用于两个大数的相乘,原理是将大数分成两段后变成较小的数位,然后做3次乘法,并附带少量的加法操作和移位操作。
现有两个大数,x,y。
首先将x,y分别拆开成为两部分,可得a,b,c,d。他们的关系如下:
x = a * 10m + b;
y = c * 10m + d。其中m为正整数,m < n,且a,c 小于 10m。
那么
xy = (a * 10^m + b)(c * 10^m + d)
xy = (10^m * a + b)(10^m * c + d)
xy = 10^2m * ac + 10^m * ad + 10^m * bc + bd
xy = 10^2m * ac + 10^m * (ad + bc) + bd
设:
z0 = ac
z1 = ad + bc
z2 = bd
利用(a + b)(c + d) = ac + ad + dc + bd 可以求得 z1
z1 = ad + bc = (a + b)(c + d) - ac - bd
最终
xy = 10^2m * ac + 10^m * ((a + b)(c + d) - ac - bd) + bd
实例展示
设x = 12345,y=6789,令m=3。那么有:
12345 = 12 * 1000 + 345;
6789 = 6 * 1000 + 789。
下面计算:
z2 = 12 * 6 = 72;
z0 = 345 * 789 = 272205;
z1 = (12 + 345) * (6 + 789) - z2 - z0 = 11538。
然后我们按照移位公式(xy = z2 * 10^(2m) + z1 * 10^(m) + z0)可得:
xy = 72 * 1000000 + 11538 * 1000 + 272205 = 83810205。