快速幂模板

Posted 2021-01-23 jcrl

math中的pow是底数不停地相乘，这样数据大了（一般都会让你模）不仅可能超空间还可能超时。

所以，快速幂诞生了；

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
template <typename T>inline void read(T &a)
{
    char ch=getchar();bool f=0;a=0;
    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){a=a*10+ch-‘0‘;ch=getchar();}
    if(f==1)a=-a;
}
 long long mod;
 long long mi( long long a, long long b)
{
     long long tot=1,base=a;
    while(b)
    {
        if(b&1)
        {
            tot*=base;
            tot%=mod;
        }
        base*=base;
        base%=mod;
        b>>=1;
    }
    return tot%mod;//一般网上都不会mod，但我强烈建议这里要模一下，如果b直接为0，岂不是就不模了？？？
}
int main()
{
 long long a,b;
read(a);
read(b);
read(mod);
cout<<a<<"^"<<b<<" mod "<<mod<<"="<<mi(a,b);
}

把long long换成自己自己定义的类型就行了