快速幂

快速幂

#include<iostream>
using namespace std;
const int _mod=1e7;
long long mul_mod(long long a,long long b,long long m){
    long long a1=a/_mod,a2=a%_mod,b1=b/_mod,b2=b%_mod,res=0;
    res=a1*b1%m*_mod%m*_mod%m;
    res=(res+a1*b2%m*_mod%m);
    res=(res+a2*b1%m*_mod%m);
    res=(res+a2*b2%m);
    return res%m;
}
long long pow_mod(long long a,long long b,long long m){
       long long t=a%m,res=1;
       while(b>=1){
           if(b&1) res=mul_mod(res,t,m);
           t=mul_mod(t,t,m);
           b>>=1;
       }
       return res;
}
int main(){
}