Given a n × n matrix A and a positive integer k, find the sum S = A + A² + A³ + … + A^k.

The input contains exactly one test case. The first line of input contains three positive integers n (n ≤ 30), k (k ≤ 10⁹) and m (m < 10⁴). Then follow n lines each containing n nonnegative integers below 32,768, giving A’s elements in row-major order.

Output the elements of S modulo m in the same way as A is given.

#include<iostream>
#include<string.h>
#include<map>
#include<cstdio>
#include<cstring>
#include<stdio.h>
#include<cmath>
#include<ctype.h>
#include<math.h>
#include<algorithm>
#include<set>
#include<queue>
typedef long long ll;
using namespace std;
const ll mod=1000;
const int maxn=30+10;
const int maxk=5e3+10;
const int maxx=1e4+10;
const ll maxe=1000+10;
#define INF 0x3f3f3f3f3f3f
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
struct matrix
{int a[maxn<<1][maxn<<1];
}ans,res,c;
int n,k,m;
int b[maxn<<1][maxn<<1];
void init()//这里全都是给矩阵赋值的过程
{memset(ans.a,0,sizeof(ans.a));memset(res.a,0,sizeof(res.a));memset(c.a,0,sizeof(c.a));for(int i=0;i<n;i++){for(int j=0;j<n;j++)cin>>b[i][j];//输入的矩阵
    }for(int i=n;i<n*2;i++){for(int j=0;j<n;j++){c.a[i][j]=b[i-n][j];//代表第二个矩阵，也就是不变的那个
        }}for(int i=0;i<n;i++){res.a[i][i]=res.a[i][i+n]=1;//单位矩阵的赋值
    }for(int i=n;i<n*2;i++){for(int j=n;j<n*2;j++){res.a[i][j]=b[i-n][j-n];//转移矩阵
        }}for(int i=0;i<n*2;i++)ans.a[i][i]=1;//单位矩阵
}
matrix multiply(matrix x,matrix y)
{matrix temp;memset(temp.a,0,sizeof(temp));for(int i=0;i<n*2;i++){for(int j=0;j<n*2;j++){for(int l=0;l<n*2;l++){temp.a[i][j]=(temp.a[i][j]+x.a[i][l]*y.a[l][j]%m)%m;}}}return temp;
}
void Quickpow()//这里每个矩阵快速幂都是一样的
{while(k){if(k&1)ans=multiply(ans,res);res=multiply(res,res);k>>=1;}
}
void solve()//这就是最后两个矩阵相乘的过程了
{ans=multiply(ans,c);for(int i=0;i<n;i++){for(int j=0;j<n;j++)cout<<ans.a[i][j]<<" ";//最后只要输出左上角那个矩阵就好了cout<<endl;}
}
int main()
{cin>>n>>k>>m;init();Quickpow();solve();return 0;
}