xmuC语言程序实践week 1 大作业

参考上文 

参考上文

#include<bits/stdc++.h>
#define LL long long
using namespace std;
inline int read()//快读
{
    char ch=getchar();int s=0,w=1;
    while(ch<48||ch>57){if(ch=='-')w=-1;ch=getchar();}
    while(ch>=48&&ch<=57){s=(s<<1)+(s<<3)+(ch^48);ch=getchar();}
    return s*w;
}
inline void write(int x)//快写
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+48);
}
struct Mat   //定义矩阵结构体
{
   int  m[101][101];
};
int n,b,m;   //n*n方阵 A，求矩阵A^b%m
Mat a,e; //a是输入的矩阵，e是输出的矩阵
Mat Mul(Mat x,Mat y)
{
    Mat c;
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j){
            c.m[i][j] = 0;
        }
    }
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j){
            for(int k=1;k<=n;++k){
                c.m[i][j] += (x.m[i][k]*y.m[k][j]);
                c.m[i][j]%=m;
            }
        }
    }
    return c;
}
Mat pow(Mat x,int  y)
{
    Mat ans;
        for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
           i==j?ans.m[i][j]=1:ans.m[i][j]=0;//单位矩阵*任何矩阵=任何矩阵本身 单位矩阵定义：对角线上元素为1，其他为0
        }
    }
    /* for(int i=1;i<=n;i++)      //注释掉部分可供调试
    {
        for(int j=1;j<=n;j++)
        {
           cout<<ans.m[i][j]<<" ";
        }cout<<endl;
    }*/
    if(!y)
    {
        ans=Mul(ans,x);
        return ans;
    }
    while(y)   //矩阵快速幂模板
        {
        if(y&1)
        {
            ans=Mul(ans,x);
              /*for(int i=1;i<=n;i++)      //注释掉部分可供调试
    {
        for(int j=1;j<=n;j++)
        {
           cout<<ans.m[i][j]<<" ";
        }cout<<endl;
    }*/
        }
        x = Mul(x,x);
        y>>=1;
    }
    return ans;
}

int main()
{
    n=read(),b=read(),m=read();
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
           a.m[i][j]=read();
        }
    }
   /*  for(int i=1;i<=2;i++)
    {
        for(int j=1;j<=2;j++)
        {
          cout<<a.m[i][j]<<" ";
        }
        cout<<endl;
    }*/
    Mat k=pow(a,b);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
          j==n?cout<<k.m[i][j]%m:cout<<k.m[i][j]%m<<" ";//注意格式，最后一个没有空格输出
        }
        cout<<endl;
    }
    return 0;
}