方阵乘法、快速幂

#include<cstdio>
#include<cstring>
const int maxn = 100;
typedef struct{
   long long Mat[maxn][maxn];
}MAT;
long long MOD,d; // d是方阵的行数或者列数
// 方阵的乘法，模MOD
MAT mul(MAT a,MAT b){
   MAT c;
   memset(c.Mat,0,sizeof(c.Mat));
   for(int i = 1 ; i <= d ; i++)
   for(int j = 1 ; j <= d ; j++)
   for(int k = 1 ; k <= d ; k++)
   c.Mat[i][j] = (c.Mat[i][j] + a.Mat[i][k] * b.Mat[k][j]) % MOD;
   return c;
}
// 方阵快速幂
MAT Quick_pow(MAT a,long long b){
   MAT c; memset(c.Mat,0,sizeof(c.Mat));
   for(int i = 1 ; i <= d ; i++) c.Mat[i][i] = 1;
   while(b){
      if(b & 1) c = mul(c,a);
      a = mul(a,a);
      b >>= 1;
   }
   return c;
}