矩阵快速幂模板

int mod;
const int NUM;//定义矩阵能表示的最大维数
int N;//N表示矩阵的维数，以下的矩阵加法、乘法、快速幂都是按N维矩阵运算的
struct Mat{
  //矩阵的类
    LL a[NUM][NUM];
    void init()//将其初始化为单位矩阵
    {
        memset(a,0,sizeof(a));
        for(int i=0;i<NUM;i++)
        {
            a[i][i]=1;
        }
    }
};
Mat add(Mat a,Mat b)//(a+b)%mod  矩阵加法
{
    Mat ans;
    for(