快速幂&快速乘&矩阵快速幂

指数快速幂：

ll quickPower(ll a,ll b){
    ll ans=1,base=a;
    while(b>0){
        if(b&1)
            ans*=base;
        base*=base;
        b>>=1;
    }
    return ans;
}

矩阵快速幂：

思路：
和整数快速幂一样，唯一不同的就是存放结果的矩阵初始值为单位矩阵，通过重载运算符*后，代码可以大大简化。

#include<iostream>  
#include<cstdio>  
#include<cstring>  
using namespace std;  
const int maxn = 105;   
const int M = 1e9+7;  

struct Matrix{  
    int n,m;  
    int v[maxn][maxn];  
    Matrix(int n,int m):n(n),m(m){}  
    void init(){                    //初始化矩阵   
        memset(v,0,sizeof v);  
    }  
    Matrix operator* (const Matrix B) const {  
        Matrix C(n,B.m);            //用来存放答案   
        C.init();  
        for(int i = 0;i < n;i++)  
        for(int j = 0;j < B.m;j++)  
        for(int k = 0;k < m;k++)       
            C.v[i][j] = (C.v[i][j] + v[i][k]*B.v[k][j] % M) % M;  
        return C;  
    }  
    void print(){//输出该矩阵，用来测试   
        for(int i = 0;i < n;i++){  
        for(int j = 0;j < m;j++)  
            cout << v[i][j] << " ";  
        cout << endl;  
        }  
    }  
};  

Matrix q_pow(Matrix& A,int b,int M){  
    Matrix C(A.n,A.m);  
    C.init();  
    for(int i = 0;i < C.n;i++)   C.v[i][i] = 1;  
    while(b){  
        if(b&1) C = C*A;  
        A = A*A;  
        b >>= 1;  
    }  
    return C;  
}  

大数的快速幂（乘法快速幂）

ll mul(ll a,ll b,ll mod){//快速乘
	ll ans=0,res=a;
	while(b){
		if(b&1) ans=(ans+res)%mod;
		res=(res+res)%mod;
		b>>=1;
	}
	return ans;
}


ll quickPower(ll a,ll b,ll mod){
    ll ans=1,base=a;
    while(b>0){
        if(b&1)
            ans=mul(ans,base,mod);
        base=mul(base,base,mod);
        b>>=1;
    }
    return ans;
}