what's 快速幂?
相信大家都看得懂这是什么:
a^b mod p=?
注:a^b表示a的b次方,mod即求余
时间复杂度为O(b)的方法相信大家都能写出:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a,b,p;
int main(){
cin>>a>>b>>p;
for(ll i=1;i<=b;i++)
a=(a%p*a%p)%p;
cout<<a;
return 0;
}
很明显的暴力相乘,当然这种方法还是有点弱,除非数据很水qwq。
但是,接下来要介绍的就是时间复杂度为O(log n(对数底数为2))的方法:
上个小学一年级的题目:
7 1 * 2 9 ———————————————————— 6 3 9 1 4 2 ———————————————————— 2 0 5 9
如果我们依照竖式的思路的话,先算71×9,这样用加法就算9次,再算71×20,71×2算2次。还剩个10,很简单,往左移一位即可,时间复杂度O(log n(对数底数为10)),当然这个时间复杂度常数会非常大,假设计算99999×99999,那么常数就非常大,因为要加9次,所以我们不妨用二进制来算,将两个因数的最小的一个变为二进制数,因为99999=99999,所以我们就把99999化为二进制数,则二进制数为: 11000011010011111
从最右边的第一位开始,如果是1,就加上另外一个因数,反之则不加,然后再将这个因数加上自己变为2倍的自己。
这样的话时间复杂度就会降到O(log n(对数底数为2)),而且常数也很小。
但是一个是O(log n(对数底数为2))常数小,一个是O(log n(对数底数为10))常数大,到底哪个效率高。
要考虑效率谁高,我们就得从他们的最坏的时间效率算起。
O(log n(对数底数为2))常数最大是1,因为最多加1次。O(log n(对数底数为10)常数最大是9,因为最多加9次
好了,现在我们可以得出来两个式子 log n(对数底数为2)和9 log n(对数底数为10)我现在得比较这两个数的大小
如果我将这两个数同时除以9,就可以得到:
(log n (对数底数为2))/9 log n(对数底数为10)
其实就是
log n (对数底数为512) log n(对数底数为10)
这个就很好想了,因为512>10,所以还是O(log n(对数底数为2))效率高
但是有没有可能三进制,四进制……效率比二进制还高。
其实就是证明二进制效率是不是最高的。
其实就是证明2^(n-1)>n 即2^n>2n (n>2)
画个函数图像就行了,当n=2时函数f(x)=2^n 和 f(x)=2n 相交,等到n>2 时 2^n就开始指数爆炸单调递增了。而2n还是按照原来速度单调递增
(百度百科直接抄的)因为本蒟蒻是一个数学蒟蒻,所以暂时无法直接证明。望大家谅解。(
说实话函数什么的也是我上网查的)
上代码:
#include<bits/stdc++.h>
using namespace std;
long long b,p,k,ans=1;
int main()
{
cin>>b>>p>>k;
while(p>0)
{
if(p&1)
ans=ans%k*b%k;
b=b%k*b%k;
p=p>>1;
}
cout<<ans%k;
return 0;
}
阿这,不会有人看不懂位运算吧,感谢提供者TXW。
那就——拜拜喽qwq