题目背景
矩阵快速幂
题目描述
给定n*n的矩阵A,求A^k
输入输出格式
输入格式:
第一行,n,k
第2至n+1行,每行n个数,第i+1行第j个数表示矩阵第i行第j列的元素
输出格式:
输出A^k
共n行,每行n个数,第i行第j个数表示矩阵第i行第j列的元素,每个元素模10^9+7
输入输出样例
输入样例#1:
2 1
1 1
1 1
输出样例#1:
1 1
1 1
n<=100, k<=10^12, |矩阵元素|<=1000 算法:矩阵快速幂
题解:
代码:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#define LL long long
#define MOD 1000000007
using namespace std;
LL a[109][109],ans[109][109],bak[109][109],K;
int n;
void Fast_Pow()
{
while(K)
{
if(K%2)
{
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
bak[i][j]=ans[i][j],ans[i][j]=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++)
ans[i][j]=(ans[i][j]+bak[i][k]*a[k][j]%MOD)%MOD;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
bak[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]+bak[i][k]*bak[k][j]%MOD)%MOD;
K/=2;
}
}
int main()
{
scanf("%d%lld",&n,&K);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
scanf("%lld",&a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
ans[i][j]=a[i][j];
K--;Fast_Pow();
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++) printf("%lld ",ans[i][j]);
puts("");
}
return 0;
}