分治法求兩個大整數的乘積

分治法求兩個大整數的乘積<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" />

 

分治法求兩個大整數的乘積 ，把兩個n位的大整數u.v分割成長度爲n/3的三段，可以用五次n/3位整數的乘法求得uv的值。按此思想設計一個兩個大整數乘積的分治算法

 

把U分成U1、U2、U3三段，即U=U1*2^2n/3+U2*2^n/3+U3

把V分成V1、V2、V3三段，即V=V1*2^2n/3+V2*2^n/3+V3

則：

UV=U1V1*2^4n/3+(U2V1+U1V2)*2ⁿ+(U1V3+U2V2+U3V1)*2^2n/3+(U3V2+U2V3)*2^n/3+U3V3

  由上邊式子知：U1V1和U3V3是避免不了的要進行乘法運算。

  另外在其餘三個式子中要通過拆分消除缺少的分量。

  例如在U1V3+U2V2+U3V1中要在所有分量中消除U1V1 U1V2 U2V1 U2V3 U3V2 U3V3。

  則可以構造(U1+U2+U3)*(V1+V2+V3)和(U1-U2+U3)*(V1-V2+V3)兩個式子通過正負抵消在所有分量中消除多餘分量。

  在U2V1+U1V2和U3V2+U2V3中還要構造彌補抵消U2V2的式子，通過演算給分量配係數可得構造式子爲(U1-2U2+4U3)*(V1-2V2+4V3)。

綜上：可令

M1=U1V1

M2=U3V3

M3=(U1+U2+U3)*(V1+V2+V3)

M4=(U1-U2+U3)*(V1-V2+V3)

M5=(U1-2U2+4U3)*(V1-2V2+4V3)

  則有： U1V1=M1

         U2V1+U1V2=(<?xml:namespace prefix = st1 ns = "urn:schemas-microsoft-com:office:smarttags" />3M1-12M2+2M3-6M4+M5)/6

         U1V3+U2V2+U3V1=(M3+M4)/2-M1-M2

         U3V2+U2V3=(-3M1+12M2+M3+3M4-M5)/6

         U3V3=M2

M1,M2,M3,M4,M5即爲五次n/3位整數的乘法。