快速幂

---

1 概述

  快速幂可以从二进制的角度理解。比如要计算$a^b$，其中指数$b={\sum }_{i=0}{b}_i\times 2^i$，而$a^{2^0},a^{2^1},a^{2^2},...$在$\mathcal{O}(\log (b))$内就能全部算出来，那么计算$a^b$的时间复杂度也就变成了（$\log b$）。

---

2 整数快速幂

2.1 C++实现

int qpow(int a,int b){
    int res=1,base=a;
    while(b){
        if(b&1) res*=base;
        base*=base;
        b>>=1;
    }
    return res;
}

---

3 矩阵快速幂

  将整数替换成矩阵即可，但是矩阵要实现乘法运算。具体参考方阵四则运算的代码：link

3.1 C++实现

const int MAXN=15;
const int MOD=7;
int A,B,n,sz;
struct Matrix{
    int dt[MAXN][MAXN];
    Matrix(int tp=0){
        for(int i=1;i<=sz;i++){
            for(int j=1;j<=sz;j++){
                if(i==j) dt[i][j]=tp;
                else dt[i][j]=0;
            }
        }
    }
    Matrix operator+(const Matrix& a){
        Matrix res;
        for(int i=1;i<=sz;i++){
            for(int j=1;j<=sz;j++){
                res.dt[i][j]=(dt[i][j]+a.dt[i][j])%MOD;
            }
        }
        return res;
    }
    Matrix operator-(const Matrix& a){
        Matrix res;
        for(int i=1;i<=sz;i++){
            for(int j=1;j<=sz;j++){
                res.dt[i][j]=(dt[i][j]-a.dt[i][j])%MOD;
            }
        }
        return res;
    }
    Matrix operator*(const Matrix& a){
        Matrix res;
        for(int i=1;i<=sz;i++){
            for(int j=1;j<=sz;j++){
                for(int k=1;k<=sz;k++){
                    res.dt[i][j]=(res.dt[i][j]+dt[i][k]*a.dt[k][j])%MOD;
                }
            }
        }
        return res;
    }
};
Matrix qpow(const Matrix& a,int b){
    Matrix res(1),rem=a;
    while(b){
        if(b&1) res=res*rem;
        rem=rem*rem;
        b>>=1;
    }
    return res;
}