最后一个点过不了
查看原帖
最后一个点过不了
1389794
Lei_ZT楼主2024/7/16 09:43
#include<bits/stdc++.h>
using namespace std;

long long a[999][999],sum[999][999],c[999][999];
const int m=1e9+7;
void jzcf2(long long n){
	for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      c[i][j]=sum[i][j],sum[i][j]=0;
    
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      for(int k=1;k<=n;k++)
      	sum[i][j]=(sum[i][j]+(c[i][k]*a[k][j])%m)%m;
}
void jzcf1(long long n){
	for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      c[i][j]=a[i][j],a[i][j]=0;
    
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      for(int k=1;k<=n;k++)
       a[i][j]=(a[i][j]+(c[i][k]*c[k][j])%m)%m;
}
void fastpow(long long n,long long k){
	while(k){
		if(k%2==1){
			jzcf2(n);
		}
		k/=2;
		jzcf1(n);
	}
}
int main(){long long n,k;
	cin>>n>>k;
	if(k==0){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				cout<<1<<" ";
			}
			cout<<endl;
		}
		return 0;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cin>>a[i][j];
			sum[i][j]=a[i][j];
		}
	}
	fastpow(n,k-1);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cout<<sum[i][j]<<" ";
		}
		cout<<endl;
	}
	return 0;
}
2024/7/16 09:43
加载中...