由此可推知:
證明了以上的公式以後,我們可以先讓a關於c取餘,這樣可以大大減少a的大小
於是我們做出了初次的改進:
long long ans = 1;//power(a,b)
long long a = a % c; //加上這一句
for(int i = 1;i<=b;i++)
ans = ans % c;
快速冪演算法依賴於以下明顯的公式
1.如果b是偶數,我們可以記k = (a^2)mod c,那麼求((k)^b/2) mod c就可以了。
2.如果b是奇數,我們也可以記k = (a^2)mod c,那麼求((k)^b/2 mod c * a) mod c就可以了。
long long ans=1;
a=a%c;
if(b%2==1)
k=(a*a)%c; //我們取a^2而不是a
for(int i=1;i<=b/2;i++)
ans = ans % c;
當然,這樣子治標不治本。但我們可以看到,當我們令k = (a * a) mod c時,狀態已經發生了變化。
我們所要求的最終結果即為:((k)^b/2) mod c而不是按原來的方法爆破(a^b)mod c,所以我們發現這個過程是可以迭代下去的。當然,對於奇數的情形會多出一項a mod c,所以為了完成迭代。
當b是奇數時,我們通過ans = (ans * a) % c;來彌補多出來的這一項,此時剩餘的部分就可以進行迭代了。
#include #include using namespace std;
/*樸素演算法*/
/*表示a的b次冪然後對c取餘的結果*/
int power1(int a, int b, int c)
/*快速冪演算法*/
int power2(int a, int b, int c)
return res;
}int main()
return 0;
}
copy的howe_young的blog
copy的howe_young的blog
快速冪 快速冪取模
快速冪的思想在於快速求解高冪指數的冪運算 複雜度為o log2n 與樸素運算相比有很大的改進 接下來給出 其中有詳解 include include using namespace std typedef long long ll ll pow1 int a,int b 最常規的方法 將冪指數轉化為...
快速冪 快速冪取模
求x m 一般方法是 xm x xm 1,這樣需要做m次乘法,未免過慢。加速方法有兩種。1.基於當m為偶數時,xm x2 m 2 當m為奇數時,xm x xm 1。顯然當m為偶數時m會減半,當m為奇數時,下次就是偶數。m可以很快收斂到0.表示冪 2.將m看成二進位制串mkmk 1 m1m0,那麼xm...
快速冪 快速冪取模
原文 快速冪這個東西比較好理解,但實現起來到不老好辦,記了幾次老是忘,今天把它系統的總結一下防止忘記。首先,快速冪的目的就是做到快速求冪,假設我們要求a b,按照樸素演算法就是把a連乘b次,這樣一來時間複雜度是o b 也即是o n 級別,快速冪能做到o logn 快了好多好多。它的原理如下 假設我們...