快速幂（log2 p）模板

#include <cstdio>
int main()
{
    long long b, p, m;//虽然看似是int范围，但如果b = 2^32 - 2 = m - 1，那么就会越界
    scanf("%lld%lld%lld", &b, &p, &m);
    long long ans = 1;
    printf("%lld^%lld mod %lld=", b, p, m);//使用这个输入输出的时候，是比较接近底层的。所以要自己考虑类型的问题
    //b %= m;//由于已经改成ll了，所以这个保护不需要（也本没有实效
    while (p > 0)
    {
        if (p & 1)//查看对应二进制位上是否为1，确定乘不乘这个b^{p^n}
            ans = ans * b % m;//少用 *= %=，它们会把数据放回再取出。 应当利用好缓存
        b = b * b % m;//维护{p^n}
        p >>= 1;
    }
    printf("%lld", ans % m);//%m 要考虑没有进入过循环的情况
    return 0;
}