一个我自己写的矩阵快速幂模板

/* n阶方形矩阵快速幂模板，如果m*n的矩阵可以考虑将m，n存入结构体matrix中
   函数分为3个，分别是矩阵相乘，转化单位矩阵，快速幂。
   测试结果  AC  
   BY SHU_ONISAC */
#include
#include
using namespace std;

struct matrix{
	long long a[15][15];
};

matrix matrix_muti(matrix a,matrix b,int n)//矩阵相乘（这个只要注意i，j，k的顺序就行了） 
{
	matrix c;
	memset(c.a,0,sizeof(c.a));
	for(int i=0;i>1（相当于k=k/2），
	  k最终一定会先变为1再变为0，所以最终结果可以存在r中。*/
	  a=matrix_muti(a,a,n);
	  k=k>>1;
	}
	return r;
}

int main()
{   int n,k;
    matrix ori,ans;
    while(cin>>n>>k)
    {
    	for(int i=0;i>ori.a[i][j];
   
    	ans=fast_power(ori,n,k);//运行快速幂 

    	for(int i=0;i