1)大整数乘法
书上讲述的是2个位数都为n位的数相乘,我这里扩展一下,写一个更为一般的方法,任意位数的2个数相乘:
对于任意位数的2个数相乘a * b,写成:
a = a1 * 10^(n1/2) + a0 -----n1为a的位数
b = b1 * 10^(n2/2) + b0 -----n2为b的位数
分治策略就是基于以上变换,将a,b写成前一半数字和后一半数字相加的形式,例如若a = 5423678,那么a1 = 542,a0 = 3678(注意若不是偶数截取较小一半)
这样a和b相乘就可以写为:a * b = { a1 * 10^(n1/2) + a0 } * { b1 * 10^(n2/2) + b0 }
展开后整理得: a * b = a1*b1 * 10^[ (n1+n2)/2 ] + a1*b0 * 10^(n1/2) + a0*b1 * 10^(n2/2) + a0*b0 四项
这样就很容易递归的来求a * b,如果你嫌分解后的数还太大,就可以继续分解。(你可以自己规定在何时结束递归)
实现方法:我们定义一个支持方法Mul(String s1,String s2),用于在结束递归时(在本例中,我定义有一个数是1位时结束递归,直接用普通乘法)计算两个字符串的乘积(为了表示大数,用字符串来接受参数)。有了这个支持方法,分治递归实现两个大数乘法的实现如下:
public static long Mutiply(String a,String b) // 用字符串读入2个大整数
{
long result = 0 ;
if (a.length() == 1 || b.length() == 1 ) // 递归结束的条件
result = Mul(a,b);
else // 如果2个字符串的长度都 >= 2
{
String a1 = a.substring( 0 , a.length() / 2 ); // 截取前一半的字符串(较短的一半)
String a0 = a.substring(a1.length(), a.length()); // 截取后一半的字符串
// System.out.println(a1);
// System.out.println(a0);
String b1 = b.substring( 0 , b.length() / 2 );
String b0 = b.substring(b1.length(), b.length());
// 分治的思想将整数写成这样: a = a1 * 10^(n1/2) + a0, b = b1 * 10^(n2/2),相乘展开得到以下四项
// 其中n1,n2为2个整数a,b的位数
result = ( long ) (Mutiply(a1,b1) * Math.pow( 10 , a0.length() + b0.length())
+ Mutiply(a1,b0) * Math.pow( 10 , a0.length()) + Mutiply(a0,b1) * Math.pow( 10 , b0.length())
+ Mutiply(a0,b0));
}
return result;
}
可以看到如果我们2个字符串都大于1位,那么我们就继续分解a和b。
完整的代码和测试的例子:
95211154和903923的乘积:
86063551957142 86063551957142
true
注意,写这个算法仅仅是看看分治是怎样实现的,我的测试用例也算不上是大整数,真正的大整数在普通的计算机上会溢出,计算不出来。 一些研究人员的实验显示,在大于600位的整数相乘时,分治策略的性能才超过普通的算法 ,显然,这么大的数字在常规的计算机上无法运算。
扩展一点,大整数加法怎么做?----------大整数乘法是为了提高效率,减少乘法运算次数而设计的,这里所说的大整数加法不是效率的问题,而是计算机无法表示的问题。例如对于100位甚至1000位的数字,计算机的表示范围显然无法达到,那么应该用数组来存放大整数的每一位,对于进位,直接向前一位加一即可。