Karatsuba’s 大整数乘法的三种方法(python)

1.问题阐述

Karatsuba’s 大整数乘法原理是将大数分成两段后变成较小的数位，然后做3次乘法，并附带少量的加法操作和移位操作。

分段工作是将大整数X、Y（分别为nx,ny位）分解为a、b、c、d。值得注意的是如果位数nx或ny为奇数，则a为前n/2+1或m/2+1位，n/2或m/2向下取整。

方法一：

为了降低复杂度可以利用如下公式计算ad+bc：

由此可以推出以下两种方法：

方法二：

方法三：

2.代码

import time

def Multiply1(X,Y,nx,ny):
    if nx==1 or ny==1:
        return X*Y
    else:
        m=int(max(nx,ny)/2)
        a,b,c,d=int(X/(10**m)),int(X%(10**m)),int(Y/(10**m)),int(Y%(10**m))
        e,f=Multiply1(a,c,len(str(a)),len(str(c))),Multiply1(b,d,len(str(b)),len(str(d)))
        g,h=Multiply1(b,c,len(str(b)),len(str(c))),Multiply1(a,d,len(str(a)),len(str(d)))
        return 10**(2*m)*e+(g+h)*(10**m)+f

def Multiply2(X,Y,nx,ny):
    if nx==1 or ny==1:
        return X*Y
    else:
        m=int(max(nx,ny)/2)
        a,b,c,d=int(X/(10**m)),int(X%(10**m)),int(Y/(10**m)),int(Y%(10**m))
        e,f=Multiply2(a,c,len(str(a)),len(str(c))),Multiply2(b,d,len(str(b)),len(str(d)))
        return 10**(2*m)*e+((a+b)*(c+d)-e-f)*(10**m)+f

def Multiply3(X,Y,nx,ny):
    if nx==1 or ny==1:
        return X*Y
    else:
        m=int(max(nx,ny)/2)
        a,b,c,d=int(X/(10**m)),int(X%(10**m)),int(Y/(10**m)),int(Y%(10**m))
        e,f=Multiply3(a,c,len(str(a)),len(str(c))),Multiply3(b,d,len(str(b)),len(str(d)))
        return 10**(2*m)*e+(e+f-(a-b)*(c-d))*(10**m)+f

X=int(input('请输入任意整数X：'))
Y=int(input('请输入任意整数Y：'))
start_time = time.perf_counter()
result=Multiply1(X,Y,len(str(X)),len(str(Y)))
tt= time.perf_counter() - start_time
print(result)
print('使用方法一耗时',tt*1000,'ms')

start_time = time.perf_counter()
result=Multiply2(X,Y,len(str(X)),len(str(Y)))
tt= time.perf_counter() - start_time
print('使用方法二耗时',tt*1000,'ms')

start_time = time.perf_counter()
result=Multiply3(X,Y,len(str(X)),len(str(Y)))
tt= time.perf_counter() - start_time
print('使用方法三耗时',tt*1000,'ms')

测试数据1：

测试数据2：

3.总结

Karatsuba’s 大整数乘法的分治策略其实目的都是减少乘法，通过加减法代替乘法，在数量级面前，加减法远快于乘法。