008-大整数乘法-分治法-《算法设计技巧与分析》M.H.A学习笔记

设u和v是两个n位二进制整数，传统乘法算法需要Θ(n²)的时间来计算u和v的乘积。

但我们还有一个Θ(n^1.59)的改进算法。

 

算法分析：

对于n位二进制u和v的乘法，我们有如下的规则：

 

当然这个规则对于其他进制的乘法也是成立的，为了便于理解，代码部分我们计算十进制的乘法，但算法分析部分依旧使用二进制的乘法。

分析过程书里写得很明白，就直接截图了：

 

 

C++代码：

/************************************************************************/
//函数功能：分治法求两个N为的整数的乘积
//输入参数：X,Y分别为两个N为整数
//算法思想：
//时间复杂度为：T(n)=O(nlog3)=O(n1.59)
/************************************************************************/
#define SIGN(A) ((A > 0) ? 1 : -1)
int IntegerMultiply(int X, int Y, int N)
{
int sign = SIGN(X) * SIGN(Y);
int x = abs(X);
int y = abs(Y);
if((0 == x) || (0 == y))
return 0;
if (1 == N)
return x*y;
else
{
int XL = x / (int)pow(10., (int)N/2);
int XR = x - XL * (int)pow(10., N/2);
int YL = y / (int)pow(10., (int)N/2);
int YR = y - YL * (int)pow(10., N/2);
int XLYL = IntegerMultiply(XL, YL, N/2);
int XRYR = IntegerMultiply(XR, YR, N/2);
int XLYRXRYL = IntegerMultiply(XL - XR, YR - YL, N/2) + XLYL + XRYR;
return sign * (XLYL * (int)pow(10., N) + XLYRXRYL * (int)pow(10., N/2) + XRYR);
}
}
int _tmain(int argc, _TCHAR* argv[])
{
int x = 1234;
int y = 4321;
cout<<"x * y = "<<IntegerMultiply(x, y, 4)<<endl;
cout<<"x * y = "<<x*y<<endl;
return 0;
}