Warshall算法简要陈述如上;那么,我们应该怎么去完成呢?
其实,我认为最好的理解方法就是跟着算法的思维去手动解决这个问题。
首先,我们知道i=1,这是随循环增加的一个数值,它是算法进行的一个限制条件,当然,它也可以作为运算时利用的条件。
当j行i列为1,就能触发Warshall进行运算。这时进行的是关系运算,设k,将其依次循环,把矩阵j行k列和i行k列关系运算后重新赋值给j行k列,循环结束后i+1。这个关系运算是怎么进行的呢?这里我们不做特别详细的说明,因为,最开始这个问题也困绕了笔者很久,但是当自己慢慢的、一个一个的去将每一个矩阵的每一个步骤列出来的时候,一切就清晰了。如下:
由上图,关系运算是能够看出来的。
当我们看懂关系运算进行后进行总结:warshall运算时,就是把第j行+第i行,送回第j行
#Warshall
print("请输入矩阵的数据") #获取矩阵数据
row=int(input("请输入矩阵行数"))
column=int(input("请输入矩阵的列数"))
print("\n")
n=row
matrix=[]
for i in range(1,n+1):
matrix.append(input("row{0}".format(i)).split(","))
print("原矩阵为:",(matrix))
a = [[0 for i in range(n)] for j in range(n)] #置新矩阵
for j in range(row): #进入循环
for i in range(n): #限制条件,i<=n
if int(matrix[j][i])==1: #进行运算
for k in range(n):
a[j][k]=int(matrix[j][k]) or int(matrix[i][k])
#第j行+第i行,送回第j行
print("M t(R)=",a)
下面给出测试集