煜虽无意逐鹿

filein=open(".in","r")
fileout=open(".out","w")

int 2147483647

矩阵快速幂

#include<bits/stdc++.h>
using namespace std;
long long mod;
long long n,m;
struct jz{
  long long f[101][101];
};
jz e,p,ans;
jz cheng(jz a,jz b){
	jz hh;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++){
		hh.f[i][j]=0;
		for(int k=1;k<=n;k++){
			hh.f[i][j]+=a.f[i][k]*b.f[k][j];
			hh.f[i][j]%=mod;
		}
	}
	return hh;
}
jz qpow(jz a,long long b){
	jz hh=p;
	while(b>0){
		if(b&1) hh=cheng(hh,a);
		a=cheng(a,a);
		b>>=1;
	}
	return hh;
}
int main(){
	scanf("%lld%lld%lld",&n,&m,&mod);//n*n矩阵,m次幂
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	scanf("%d",&e.f[i][j]);
	for(int i=1;i<=n;i++)
	for(int i=1;i<=n;i++)
	p.f[i][i]=1;
	if(mod==1) printf("0");
	if(n==0) printf("1");
	ans=qpow(e,m);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)
		printf("%lld ",ans.f[i][j]);
		printf("\n");
	}
	return 0;
}

快速幂

long long mypow(long long a,long long b){
	if(!b) return 1;
	long long f=mypow(a,b/2);
	if(b%2) return f*f%k*a%k;
	else return f*f%k;
}

stable_sort()//稳定快排

while(cin>>ws,!cin.eof())//读入直到文件末尾

ios::sync_with_stdio(false);//优化输出时间，使cin cout时间与scanf printf相差无几

map<string,int>a;
map<string,int>::iterator it;

freopen(".in","r",stdin);
freopen(".out","w",stdout);

int mkbk(string s){//将字符串转为整数
	int ls=s.length();
	int m[ls],tot=0;
	for(int i=0;i<ls;i++){
		tot*=10;
		m[ls-i-1]=s[i]-'0';
		tot+=m[ls-i-1];
	}
	return tot;
}

string add(string s1,string s2){
    long long a[10010]={},b[10010]={},c[10010]={};
    long long la=s1.length(),lb=s2.length(),lc=max(la,lb),x=0;
    string ans="";
    for(int i=0;i<la;i++)a[i]=s1[la-i-1]-'0';
    for(int i=0;i<lb;i++)b[i]=s2[lb-i-1]-'0';
    for(int i=0;i<lc;i++){
        c[i]=a[i]+b[i]+x;
        x=c[i]/10;
        c[i]=c[i]%10;
    }
    if(x)c[lc++]=1;
    for(int i=lc-1;i>=0;i--)ans+=(c[i]+'0');
    return ans;
}

string sub(string s1,string s2){
    long long a[10010]={},b[10010]={},c[10010]={};
    long long la=s1.length(),lb=s2.length(),lc=max(la,lb),x=0,f=0;
    string ans="";
    if((s1<s2&&s1.length()==s2.length())||s1.length()<s2.length()){
        f=1;
        string t=s1;
        s1=s2;
        s2=t;
        swap(la,lb);
        cout<<'-';
    }
    for(int i=0;i<la;i++)a[i]=s1[la-i-1]-'0';
    for(int i=0;i<lb;i++)b[i]=s2[lb-i-1]-'0';
    for(int i=0;i<lc;i++){
        c[i]=a[i]-b[i]-x;
        if(c[i]<0){
            x=1;
            c[i]+=10;
        }else x=0;
    }
    while(c[lc-1]==0&&lc>1)lc--;
    for(int i=lc-1;i>=0;i--)ans+=(c[i]+'0');
    return ans;
}

string mul(string s1,string s2){
    long long a[10010]={},b[10010]={},c[20010]={};
    long long la=s1.length();
    long long lb=s2.length();
    long long lc=la+lb;
    string ans="";
    for(int i=0;i<la;i++) a[i]=s1[la-1-i]-'0';
    for(int i=0;i<lb;i++) b[i]=s2[lb-1-i]-'0';
    for(int i=0;i<la;i++){ 
        for(int j=0;j<lb;j++){
            c[i+j]+=a[i]*b[j];
            c[i+j+1]+=c[i+j]/10;
            c[i+j]%=10;
        }
    }
    while(c[lc-1]==0&&lc>1)
    lc--;
    for(int i=lc-1;i>=0;i--) ans+=(c[i]+'0');
    return ans;
}

string div(string s,int b){
    int n,a[10100],c[10100],la,lc;
    la=s.length(),lc=0;
    for(int i=0;i<la;i++)a[i]=s[i]-'0';
    long long x=0;
    string ans="";
    for(int i=0;i<la;i++){
        c[i]=(x*10+a[i])/b;
        x=(x*10+a[i])%b;
    }
    while(c[lc]==0&&lc+1<la) lc++;
    for(int i=lc;i<la;i++)ans+=(c[i]+'0');
    return ans;
}

骗分过样例，暴力出奇迹，暴搜挂着机，打表出省一（>.<）