模板之矩阵快速幂（luogu P3390【模板】矩阵快速幂）

#include<cstdio>
const int c=1000000007;
long long n,k;
void read(long long &y)
{
    y=0;char x=getchar();
    while(x<‘0‘||x>‘9‘) x=getchar();
    while(x>=‘0‘&&x<=‘9‘)
    {
        y=y*10+x-‘0‘;
        x=getchar();
    }
}
struct mat
{
    long long a[105][105];
    mat operator*(const mat &b)const
    {
        mat x;
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                x.a[i][j]=0;
                for(int k=1;k<=n;k++)
                {
                    x.a[i][j]+=a[i][k]*b.a[k][j];
                    x.a[i][j]%=c;
                }
            }
        }   
        return x;
    }
}m;
mat ksm(mat x,long long y)
{
    mat ans=x,s=x;
    for(;y;y>>=1)
    {
        if(y&1) ans=ans*s;
        s=s*s;
    }
    return ans;
}
int main()
{
    read(n);read(k);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++) read(m.a[i][j]);
    }
    m=ksm(m,k-1);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++) printf("%d ",m.a[i][j]);
        printf("\n");
    }
    return 0;
}