快速幂&矩阵快速幂

快速幂：我的理解就是可以快速的计算a^b; 而这里所用的方法涉及到了二进制；

如：2^11: 11的二进制值数是1011， 那么11可以写成1*2^3+0*2^2+1*2^1+1*2^0, 即：8+0+2+1=11；因此， 2^11=2^(8+2+1)=2^8*2^2*2^1; 很明显， 前者需要连续乘11个2， 后者只需运算3次；

int power(int a, int b){
    int ans=1, res=a;
    while(b){
        if(b&1) ans*=res;//（当二进制位是1时， 要乘）
        b>>=1;
        res*=res;//(这里一定要注意， 不要乘错， 为何要乘res， 而不是a呢？——x^(2^i)=(x^(2^(i-1)))^2）
    }
    return ans;
}

然而， a^b是个很大的数， 我们通常要对某数取模：

int power(int a, int b, int c){//（对c取模）
    int ans=1, res=a%c;
    while(b){
        if(b&1) ans*=res, ans%=c;
        b>>=1;
        res*=res;
        res%=c;
        cout << res << endl;
    }
    return ans;
}

现在来看看矩阵快速幂；

node mul(node a, node b){
    node c;
    c.x=a.x;
    c.y=b.y;
    memset(c.z, 0, sizeof(c.z));
    for(int i=0; i<c.x; i++){
        for(int j=0; j<c.y; j++){
            for(int k=0; k<b.y; k++)
                c.z[i][j]+=a.z[i][k]*b.z[k][j];
        }
    }
    return c;
}
node power(node a, int n){
    node b;
    b.x=a.x;
    b.y=a.y;
    for(int i=0; i<a.x; i++){
        b.z[i][i]=1;
    }
    b.z[0][1]=b.z[1][0]=0;
    while(n>0){
        if(n&1) b=mul(b, a);
        a=mul(a, a);
        n>>=1;
    }
    return b;
}