hdu 1575 矩阵快速幂

最新推荐文章于 2022-07-29 11:41:21 发布

原创最新推荐文章于 2022-07-29 11:41:21 发布 · 485 阅读

0 ·

CC 4.0 BY-SA版权

编程同时被 2 个专栏收录

457 篇文章

订阅专栏

数学题-技巧-思维

24 篇文章

订阅专栏

本文介绍了一种使用矩阵快速幂算法高效计算矩阵的幂，并通过该算法求解矩阵的迹（对角元素之和），特别关注了在计算过程中取模运算以避免整数溢出的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<string.h>
#define N 100
int map[N][N];
int power(int n,int c) {
int sum[N][N],i,j,k,mapan[N][N];
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
sum[i][j]=map[i][j];
c--;
while(c) {
if(c&1) {
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
mapan[i][j]=sum[i][j];
memset(sum,0,sizeof(sum));
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
for(k=1;k<=n;k++)
sum[i][j]=(sum[i][j]+mapan[i][k]*map[k][j])%9973;
}
memset(mapan,0,sizeof(mapan));
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
for(k=1;k<=n;k++)
mapan[i][j]=(mapan[i][j]+map[i][k]*map[k][j])%9973;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
map[i][j]=mapan[i][j];
c/=2;
}

int su=0;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
if(i==j)
su=(su+sum[i][j])%9973;
return su;
}
int main() {
int n,t,k,i,j;
scanf("%d",&t);
while(t--) {
scanf("%d%d",&n,&k);
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
scanf("%d",&map[i][j]);
printf("%d\n",power(n,k));
}
return 0;
}