C++ 矩阵快速幂

在学这个之前学了好一阵子矩阵，仍然不太懂。。后来发现解斐波那契数列好像只需要方形矩阵就行了，就学了一下方形矩阵的乘法。

思路：首先，假定我们会快速幂。
那么如何写矩阵快速幂呢？
我们需要明确一件事，我们写的快速幂实际上是：十进制数字快速幂
十进制数字的单位是1。而矩阵的单位是单位矩阵
那么，是不是我们只要把十进制数字的单位"1"换成"单位矩阵"就可以写矩阵快速幂了？
答案是是的。不过需要注意一个细节：由于十进制和矩阵的乘法规则不一样，因此矩阵的乘法需要额外写一个方法。

上代码吧（未AC）

#include<cstdio>
#include<iostream>
#define ll long long
#define mod 1000000007
using namespace std;
struct Mat{
	ll m[101][101];
};
Mat initMat(int n);
Mat Mul(Mat a,Mat b,int n);
Mat poww(Mat a,ll b);
Mat a,e;//a:输入 e:单位矩阵 
int n;//矩阵长度 (每行n个数 ) 
int main(){
	Mat a;
	int k;//k是Mat的指数 
	cin>>n>>k;//input 矩阵长度以及指数
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			cin>>a.m[i][j];
		}
	}
	
	//初始化e
	f