21计科程一帆


私信TA

用户名:uq_88617846948

访问量:5255

签 名:

搞哥毛哥在上,俺寻思俺是一个最大最强的技术小子

等  级
排  名 958
经  验 3423
参赛次数 2
文章发表 52
年  龄 19
在职情况 学生
学  校 石河子大学
专  业 计算机科学与技术

  自我简介:

憨憨一个,欢迎大佬指正

解题思路:这一题的思路相当简单,只需要接收题给矩阵,按照要求求出对应的矩阵幂就行,我们需要做的就是处理矩阵幂的函数,这里我们细分为两部分,首先做出矩阵乘积的计算函数,然后套用该函数做出矩阵幂函数,具体的实现细节比较简单,详细解释部分见于代码的对应注释

注意事项:刚开始不习惯可以将矩阵的行列初始都加上一个无用的0,这样后面理坐标的时候下标不用从0开始,不容易搞混

参考代码:

a,b=map(int,input().split())
c=[list(map(int,input().split())) for _ in range(a)]#接收题给矩阵数据
def juzhenji(A,B):#矩阵乘积计算
   n=len(A)#这里本需要处理出矩阵A的行和矩阵B的列,还有矩阵A的列或者矩阵B的行,但是题给的都是方阵,行列相等
   d=[[0]*n for _ in range(n)]#与矩阵A同规模的空矩阵,用于接收计算完成的新矩阵
   for i in range(a):#矩阵A的行
       for j in range(a):#矩阵B的列数
           for k in range(a):#矩阵A的列或者矩阵B的行

               d[i][j]+=A[i][k]*B[k][j]
   return d
def juzhenmi(A,b):#矩阵幂,只需考虑方阵就行,非方阵计算不了幂
   a=len(A)#取行/列
   result=[[1 if i==j else 0 for j in range(a)]for i in range(a)]#做一个单位矩阵
   while b>0:
       if b%2==1:
           result=juzhenji(result,A)#剩余要处理的幂数为奇数则乘以本矩阵
       A=juzhenji(A,A)#偶数则将要处理的矩阵平方
       b//=2
   return result
juzhen=juzhenmi(c,b)
for i in juzhen:
   print(" ".join(map(str,i)))

 

0.0分

1 人评分

  评论区

  • «
  • »