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 大整数乘法的分治策略其实目的都是减少乘法,通过加减法代替乘法,在数量级面前,加减法远快于乘法。