矩阵快速幂

本次学习的内容是 矩阵快速幂
分为矩阵乘法和快速幂两个部分；

介绍一下单位矩阵
左上至右下对角线为1，其余部分为0；

矩阵乘法基本步骤

题目大意
n*n矩阵，求其k次方结果；

#include<iostream>
#include<cstdio>
using namespace std;

struct matrix{
    long long m[105][105];
}matri;
const int q=1e9+7;
long long n,k;
matrix ksm(matrix a,matrix b)//矩阵乘法
{
    matrix c;
    for(int i=1;i<=n;i++)   //分别行，列推移
        for(int j=1;j<=n;j++)
        {
            c.m[i][j]=0;
            for(int p=1;p<=n;p++)    
                    c.m[i][j]=(c.m[i][j]+a.m[i][p]*b.m[p][j])%q;
        }    

    return c;
}

matrix ready(matrix a,long long k)//快速幂处理
{
    k--;
    matrix ko=a,vic=a;
    for(;k;k>>=1,vic=ksm(vic,vic))
        if(k&1)
        ko=ksm(ko,vic);
    return ko;


}
int main()
{
    matrix A,C;
    scanf("%lld%lld",&n,&k);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        scanf("%lld",&A.m[i][j]);

    C=ready(A,k);
    for(int i=1;i<=n;i++)  
    {
        for(int j=1;j<=n;j++)
        printf("%lld ",C.m[i][j]);
        printf("\n");
    }
    return 0;
}