C++算法——龟速乘

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll mod;
ll mul(ll a,ll b)
{
    ll ans=0;
    while(b)
    {
        if(b&1)
            ans=(ans+a)%mod;
        a=(a+a)%mod;
        b>>=1;
    }
    return ans;
}