时间复杂度取决于大数乘法的大整数除法.docVIP

时间复杂度取决于大整数乘法的大整数除法 技术领域 本文主要提供一种特殊的分治递归运算的大整数除法（包括整除的商值及整除后取余数）。这种算法主要应用于计算机领域，尤其是加密技术中。 论文背景 当前公开或使用的大整数除法即使采用Knuth的经典名著“The Art of Computer Programming”的第4.3.1节中的猜商值公式来计算，也需要多次试运算，所有这些大整数除法最大的缺点是时间复杂度T(n) = O(n2)这与大整数乘法相比差远了。因此减少试运算的次数，降低时间复杂度时间复杂度时间复杂度时间复杂度T(n) = O(nlog23)，如果采用其他更快的大整数乘法，本算法会更快，这就使数值很长的大整数除法的计算速度大幅提高，例如当被除数的长度为4194304个四字节，除数长度为2097152个四字节时，本算法中采用Karatsuba乘法加速时的速度为现有大整数除法速度的30倍左右。 论文内容 本文在这部分将推出一种全新的大整数除法计算理论，并利用该理论构造分治的大整数除法，在分治中充分利用大整数乘法，以达到加速目的。请看下面的推理： 设大整数运算过程中采用h（大于1的整数）进制， 正整数集合：Z={x| x＞0并且 [x]=x } 集合1：H= {x| 0≤x＜h 并且x为整数} 集合2：H2={x| 1≤x＜h 并且x为整数} 任何有限长度的大整数除法都可转化为非负大整数的除法，因此除了最后一段外，下文只涉及非负大整数除法，非负大整数除法可表示如下： C=A/B (A≥0, B＞0, C≥0并且A和B都是整数) （1） 当A=0时，显然C=0， （2） 当A＜B时，显然C的整数部分为0， （3） 当A≥B时，可通过将被除数和除数都乘以hm从而使A和B都转化为不小于hm的整数，并使结果C保持不变，其中m≥2且m ∈Z，因此（1）式在排除（2）、（3）后，可转化为下式： C=A/B (A≥B，B≥hm, m≥2, C＞0且A、B和m都是正整数) （4） 在（4）中A=(A1A2A3…Ai-1Ai)h可用多项式表示如下： A=A1hi-1＋ A2hi-2＋ A3hi-3＋…＋ Ai-1h1＋ Ai （5） 其中i∈Z，A1∈H2， Ar∈H ( 1＜r≤i,且r为整数), h为进制基数， 用AH表示A的前m项，即 AH=A1hi-1＋ A2hi-2＋ A3hi-3＋…＋ Am-1hi-(m-1)＋ Amhi-m （6） 设A=AH＋AE （7） 则有 AE=Am＋1hi-(m＋1)＋Am＋2hi-(m＋2)＋Am＋3hi-(m＋3)＋…＋Ai-1h＋Ai （8） 在（4）中B=(B1B2B3…Bj-1Bj)h可用多项式表示如下： B=B1hj-1＋ B2hj-2＋ B3hj-3＋…＋ Bj-1h1＋ Bj （9） 其中j∈Z，B1∈H2， Br∈H ( 1＜r≤j,且r为整数), h为进制基数 用BH=(B1B2B3…Bm-1Bm)h×hj-m表示B的前m项，即 BH=B1hj-1＋B2hj-2＋ B3hj-3＋…＋Bm-1hj-(m-1)＋Bmhj-m （10） 其中m＜j 设B=BH＋BE 则有 BE=Bm＋1hj-(m＋1)＋Bm＋2hj-(m＋2)＋Bm＋3hj-(m＋3)＋…＋Bj-1h＋Bj （11） 在（4）中C=(C1C2C3…Ck-1Ck)h＋fC可用多项式表示如下： C=C1hk-1＋C2hk-2＋C3hk-3＋…＋Ck-1h1＋Ck＋fC （12） 其中k∈Z，C1∈H2 ，Cr∈H ( 1＜r≤k,且r为整数), 0≤fC＜1， h为进制基数 设 W=A/BH 则W=( W1W2W3…Wn-1Wn)h＋fW可用多项式表示为 W=W1hn-1＋W2hn-2＋W3hn-3＋…＋Wn-1h1＋Wn＋fW （12_H） 其中n∈Z，W1∈H2， Wr∈H ( 1＜r≤n,且r为整数), 0≤fW＜1，其中h为进制基数 根据前面的陈述可知: A/BH≥A/B 即W≥C 下面进行理论推导的第一大步骤,讨论i、j和k的关系: C=A/B=(A1hi-1＋A