分治法妙用-1

Divide and Conquer

求斐波那契数

转化为矩阵乘法,用分治法计算,复杂度为O(logn)O(logn)

分治法妙用-1

位乘问题

设X,Y为两个n位二进制数,n=2kn = 2^k,求XY

传统算法

W(n)=O(n2)W(n) = O(n^2)

分治法

把两个数都对半截,分解为4个子问题
X=An/2+B,Y=Cn/2+DXY=AC2n+(AD+BC)2n2+BDW(n)=4W(n/2)+cnW(n)=O(nlog4)=O(n2) X = A^{n/2}+B , Y = C^{n/2}+D\\ XY = AC2^n + (AD+BC)2^{n-2}+BD\\ W(n) = 4W(n/2) + cn \\ 由主定理 W(n) = O(n^{log4})= O(n^2)

,没有本质变化,继续改进:
ACADBCBD3AD+BC=(AB)DC)+AC+BDW(n)=3W(n/2)+cnW(n)=O(nlog3)=O(n1.59) 原本有四个子问题,AC,AD,BC,BD,通过代数变换可以减少为3个 \\ AD+BC = (A-B)(D-C)+AC+BD\\ 则W(n) =3W(n/2)+cn\\ W(n) = O(n^{log3})= O(n^{1.59})

矩阵乘法

A,B为两个n阶矩阵,计算C = AB

暴力算法

不表,O(n3)O(n^3)

分治法

分治法妙用-1

分解为8个子问题,仍然是O(n3)O(n^3)

分治法改进:Strassen算法(1969)

通过简单的线性变换,将八个子问题变为七个子问题,时间复杂度:
W(n)=7W(n/2)+18(n/2)2W(n)=O(nlog27)=O(n2.8075) W(n) = 7W(n/2)+18(n/2)^2\\ 由主定理 W(n) = O(n^{log_2^7})= O(n^{2.8075})
线性变换如下:

分治法妙用-1

现代进展

分治法妙用-1

Coppersmith–Winograd 算法把 N* N大小的矩阵乘法的时间复杂度降低到了:分治法妙用-1,而2010年,Andrew Stothers再度把复杂度降低到了分治法妙用-1,一年后的2011年,Virginia Williams把复杂度最终定格为:分治法妙用-1

但这些算法常数太大,实际应用较少。

GItHub主页上的原文链接