【矩阵快速幂】

1.21矩阵快速幂模板

typedef long long LL;

const int mod =  ;

struct Matrix{
    LL v[4][4];
};

Matrix multiply(Matrix x, Matrix y, int Matrix_len);/*矩阵乘法*/
Matrix matrix_pow(Matrix a, int k, int Matrix_len);/*矩阵快速幂*/

Matrix multiply(Matrix x, Matrix y, int Matrix_len){
    Matrix z;
    memset(z.v, 0, sizeof(z.v));
    for(int i = 0; i < Matrix_len; i++){
        for(int j = 0; j < Matrix_len; j++){
            for(int k = 0; k < Matrix_len; k++){
                z.v[i][j] += x.v[i][k] * y.v[k][j];
                z.v[i][j] %= mod;
            }
        }
    }
    return z;
}
Matrix matrix_pow(Matrix a, int k, int Matrix_len){
    Matrix b;
    for(int i = 0; i < Matrix_len; i++){
        for(int j = 0; j < Matrix_len; j++){
            i == j? b.v[i][j] = 1: b.v[i][j] = 0;
        }
    }
    while(k){
        if(k & 1)
            b = multiply(b, a, 2);
        a = multiply(a, a, 2);
        k >>= 1;
    }
    return b;
}