a的b次方怎么求
pow（a, b）是數學頭文件math.h里面有的函數
可是它返回值是double類型，數據有精度誤差
那就自己寫for循環咯
LL pow(LL a, LL b){//a的b次方
LL ret = 1;
for(LL i = 1; i <= b; i ++){
ret *= a;
}
return ret;
}
完美
可是題目是b的范圍是1 <= b <= 1e9（＃°Д°）
超時，妥妥的。。。
看個例子
比如計算
2*2*2*2*2*2*2*2*2*2*2
可以這樣算
原式=4*4*4*4*4*2
=8*8*4*2
=16*4*2
你看，相同的可以先合並，減少計算步驟
如果題目說數據很大，還需要求余，那么代碼就可以這么寫
1 LL pow_mod(LL a, LL b){//a的b次方 2 if(b == 0) return 1; 3 LL ret = pow_mod(a, b/2); 4 ret = ret * ret % MOD; 5 if(b % 2 == 1) ret = ret * a % MOD; 6 return ret; 7 }
這是遞歸寫法
然后還有遞推寫法
1 LL pow_mod(LL a, LL b){//a的b次方 2 LL ret = 1; 3 while(b != 0){ 4 if(b % 2 == 1){ 5 ret = (ret * a) % MOD ; 6 } 7 a = (a * a ) % MOD ; 8 b /= 2; 9 } 10 return ret; 11 }
對於位運算熟的小盆友，還可以寫成位運算形式，速度又快，又好理解，在加一個求余p，代碼如下
1 LL pow_mod(LL a, LL b, LL p){//a的b次方求余p 2 LL ret = 1; 3 while(b){ 4 if(b & 1) ret = (ret * a) % p; 5 a = (a * a) % p; 6 b >>= 1; 7 } 8 return ret; 9 }
有了快速冪，於是，快速乘誕生了
1 LL mul(LL a, LL b, LL p){//快速乘，計算a*b%p 2 LL ret = 0; 3 while(b){ 4 if(b & 1) ret = (ret + a) % p; 5 a = (a + a) % p; 6 b >>= 1; 7 } 8 return ret; 9 }
(*´Д｀*)快速乘應該不怎么會用，無意義的東西，說不定哪天用的上
這些知識到底算不算數論呢？？？不管了（´∀｀*)