1. 程式人生 > >快速冪取模演算法

快速冪取模演算法

演算法2:另一種演算法利用了二分的思想,可以達到O(logn)。
可以把b按二進位制展開為:b = p(n)*2^n  +  p(n-1)*2^(n-1)  +…+   p(1)*2  +  p(0)
其中p(i) (0<=i<=n)為 0 或 1

這樣 a^b =  a^ (p(n)*2^n  +  p(n-1)*2^(n-1)  +...+  p(1)*2  +  p(0))
               =  a^(p(n)*2^n)  *  a^(p(n-1)*2^(n-1))  *...*  a^(p(1)*2)  *  a^p(0)
對於p(i)=0的情況, a^(p(i) * 2^(i-1) ) =  a^0  =  1,不用處理
我們要考慮的僅僅是p(i)=1的情況
化簡:a^(2^i)  = a^(2^(i-1)  * 2) = (  a^(  p(i)  *  2^(i-1)  )  )^2
(這裡很重要!!具體請參閱秦九韶演算法:http://baike.baidu.com/view/1431260.htm