快速幂模
快速冪模
- 簡(jiǎn)述
- 師從
- 普通思路
- 缺陷一:溢出
- 缺陷二:運(yùn)算次數(shù)多
- 二分化快速冪模
簡(jiǎn)述
計(jì)算 ana^nanmod p
師從
本篇是觀Vita君算法視頻后總結(jié),他是bilibili一位小up主:小學(xué)生Vita君
正所謂“生乎吾后,其聞道也亦先乎吾,吾從而師之”,誠(chéng)然如此。
【算法小知識(shí)】如何計(jì)算快速冪(上)
普通思路
(aaaaaa……aaa) % p
缺陷一:溢出
ana^nan可能會(huì)溢出
解決方案:邊乘邊模
基于 (a*b) mod p = (a mod p)(b mod p) 成立
(aaaa……)=(a%p)(a%p)(a%p)(a%p)……
缺陷二:運(yùn)算次數(shù)多
當(dāng)n較大時(shí),運(yùn)算次數(shù)很大,速度慢
解決方案:二分化思想
ana^nan=an2a^\frac{n}{2}a2n?an2a^\frac{n}{2}a2n?
an2a^\frac{n}{2}a2n?=an4a^\frac{n}{4}a4n?an4a^\frac{n}{4}a4n?
……
(當(dāng)指數(shù)為奇數(shù)時(shí),需額外乘a)
將時(shí)間復(fù)雜度O(n) → O(logn)
二分化快速冪模
#include<iostream> using namespace std;typedef unsigned long long ull; ull binpow(ull a, ull n, ull p) {if (n == 0) return 1;a %= p;ull c = binpow(a, n / 2, p);if (n % 2 != 0) return c * c % p *( a % p );return c * c % p; } int main() {ull a, n, p;cin >> a >> n >> p;cout<< binpow(a, n, p) % p<<endl; }非遞歸:
ull binpow(ull a, ull n, ull p) {ull prod = 1;while(n > 0){if (n & 1) prod = prod * a % p;a *= a % p;n >>= 1;}return prod; }總結(jié)
- 上一篇: 两数的最大公约数算法基础及优化
- 下一篇: 颐和园什么时候开门和关门