分治算法之大整数乘法--big mutiply

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。

完整的代码和测试的例子：

BigIntegerMutiply

95211154和903923的乘积：

86063551957142  86063551957142
true

注意，写这个算法仅仅是看看分治是怎样实现的，我的测试用例也算不上是大整数，真正的大整数在普通的计算机上会溢出，计算不出来。 一些研究人员的实验显示，在大于600位的整数相乘时，分治策略的性能才超过普通的算法 ，显然，这么大的数字在常规的计算机上无法运算。


扩展一点，大整数加法怎么做？----------大整数乘法是为了提高效率，减少乘法运算次数而设计的，这里所说的大整数加法不是效率的问题，而是计算机无法表示的问题。例如对于100位甚至1000位的数字，计算机的表示范围显然无法达到，那么应该用数组来存放大整数的每一位，对于进位，直接向前一位加一即可。