深夜敲模板_1——快速幂 && 矩阵的快速幂

快速幂：

///a^n%m
long long quickpow(long long a,int n,int m){
    long long ans=1;
    while(n){
        if(n&1) ans = (ans*a)%m;
        a = (a*a)%m;
        n>>=1;
    }
    return ans;
}

矩阵的快速幂：

///用结构体保存一个矩阵
struct Matrix{
    LL mat[N][N];
    ///两个矩阵乘法
    Matrix operator *(const Matrix& a) const{
        Matrix c;
        memset(c.mat,0,sizeof(c.mat));
        for(int i=0;i<N;i++){
            for(int k=0;k<N;k++) if(mat[i][k]){
                for(int j=0;j<N;j++){
                    c.mat[i][j]+=mat[i][k]*a.mat[k][j];
                }
            }
        }
        return c;
    }
    ///矩阵的幂，用二分的思想写，就是快速幂算法
    Matrix operator ^(int t) const{
        Matrix c;
        ///初始化为单位矩阵
        for(int i=0;i<N;i++){
            for(int j=0;j<N;j++){
                c.mat[i][j]=(i==j);
            }
        }
        ///就是快速幂
        Matrix a = *this;///毕竟直接敲的模板，没运行就写上去了，发现这里少了一个*号
        while(t)
        {
            if(t%2) c=c*a;
            a=a*a;
            t=t>>1;
        }
        return c;
    }
};