Washall算法可由一个图的邻接矩阵求它的可达矩阵,形式上很简单,就是一个三层循环,
问题引出:在一个图结构中,常常需要找两个节点之间有没有一条通路(通路长任意),也就是任意两点之间的可达情况,我们很自然会联想到用邻接矩阵来求可达矩阵,从而得出两点之间的可达情况。
如:
在这个图中,可以写出其邻接矩阵为:
0 1 1 1
1 0 0 0
1 0 0 1
1 0 1 0
记为矩阵A。那么很明显在图中,V2与V1存在直接通路(即长是1),V1与V4之间也有直接通路,但是V2与V4之间没有长度是1的通路,但是存在边长是2的通路,然而这在矩阵A中没有相应的体现,因为矩阵元素M(2,4) == 0,于是我们找任意图中一点Vk,只要满足V(2k)乘以V(k4)!= 0[注意:括号中指的是下标],则证明V2可以借助Vk(此图中指的是V1)与V4连通。
以此类推,当想得到任意一点与其他各点通路是2的连接的时候(这里指的是Vi,Vj两点)就可以逐个计算满足:V(ik)乘以V(kj)!= 0的,而V(ik)乘以V(kj)正好是求矩阵A乘以矩阵A的过程,所以矩阵A的平方就可以反应出通过边长为2图中各点之件的可达情况,以此类推,A的3次方、A的4次方……,分别表示了边长是3、边长是4……的可达情况,现在把所有可能长度的通路可达情况加起来(即A1方 + A2方 + A3方 +……),就可以得到近似可达矩阵。
为什么这里说是近似的可达矩阵?因为一个复杂的图中(如互联网模型),通路的长度成千上百都可能,那到底应该加到A的几次方才算有保证呢?很明显,这个算法存在局限性,对于简单的图还凑合,复杂的图就没有全面的保障了。Warshall算法就是为了求一个可靠的可达矩阵。
算法的自己的语言描述:依次遍历邻接矩阵中的所有元素(M[i,j]),比如按照先列后行进行,如果M[i,j] != 0,那么就把i行加到第j行上。