快速幂&&快速计算（a*b）mod m

  概述：起初在dalao的博客里看到这个（a*b）mod m就觉得很不解，为什么（a*b）对m求余还需要用一个函数来实现，非常的不解。想想之后才发现自己真的蠢，要是m很大的话，虽然有a,b<m，但是a*b或许就可以超出long long可以表示的范围，所以需要用特殊的方法计算。至于说为什么把快速幂和a*b的问题扯到一起，简而言之，两者之间唯一的差别就是一个是b个a相乘，一个是b个a相加，都是用位运算处理的。

  简要思路：。由于二者思路类似，所以这里只说快速幂，记得以前知道这样一个概念，用1，2，4，8...2^n的线性组合可以表示2^n以内任意的数，快速幂正是利用了这个原理。这是基于二进制的一个原理，比如11的二进制是1011，也就是11=2^0+2^1+2^3,那么a^11就可以表示成为a^(2^0+2^1+2^3)=a^2^0*a^2^1*a^2^3。实现起来就是每次用&计算b的二进制的最后一位，若为1，就说明最后的乘积项中含有a^2^x这一项，将其乘进ans。不管二进制最后一位是奇数还是偶数，都对a进行平方操作，为的是得到a^2^x。由于是乘法，可以直接对每一步进行求余数操作。

  算法的时间复杂度：log2(n)。

  快速幂代码模板：

  long long qpow(long long a,long long b){
    long long ans=1;
    long long k=a;
    while(b){
        if(b&1)ans=ans*k%1000000009;
        k=k*k%1000000009;
        b>>=1;
    }
    return ans;
 }

 这里m=1000000009，如果需要对其他数求余数，函数里面加一个参数即可。

  （a*b）mod m的模板：

long long qmul(long long a,long long b,long long m){
    a%=m;
    b%=m;
    long long ans=0;
    while (b){
        if (b&1){
            ans=(ans+a)%m;
        }
        a=(a<<=1)%m;
        b>>=1;
    }
    return ans;
}

这里，除了一开始加的mod优化，基本就是把快速幂的乘法改成了加法，当然，如果想对快速幂优化也有类似的，b%euler(m)，这个定理及其相关将在其他博客给出。