WarShall算法求矩阵传递闭包关系

离散知识

给了你一个矩阵，你如何求他的传递闭包呢？

//求出如下矩阵的传递闭包
1 0 1 0 0
1 1 1 0 0
0 0 0 0 1
1 0 0 1 0
0 0 0 0 0

利用WarShall算法求解过程如下：

从第一列元素，开始看，哪个为1，可以看出，<1,1>,<2,1>,<4,1>为1,则分别将第1、2、4行元素与第1行元素进行逻辑加，成为新的1、2、4行的元素。

1 0 1 0 0
1 1 1 0 0
0 0 0 0 1
1 0 1 1 0
0 0 0 0 0

然后看第二列，<2,2>为1，则将第2行与第2行各个元素进行逻辑加，成为新的第2行元素

1 0 1 0 0
1 1 1 0 0
0 0 0 0 1
1 0 1 1 0
0 0 0 0 0

看第三列元素，<1,3>,<2,3>,<4,3>为1，则可得

1 0 1 0 1
1 1 1 0 1
0 0 0 0 1
1 0 1 1 1
0 0 0 0 0

看第四列，<4,4>为1，则可得

1 0 1 0 1
1 1 1 0 1
0 0 0 0 1
1 0 1 1 1
0 0 0 0 0

最后一列，<1,5>,<2,5>,< 3,5>,<4,5>，最终可得

1 0 1 0 1
1 1 1 0 1
0 0 0 0 1
1 0 1 1 1
0 0 0 0 0

代码实现

for(int i=1;i<=5;i++)
		for(int j=1;j<=5;j++)
			cin>>a[i][j];
		for(int i=1;i<=5;i++)
			for(int j=1;j<=5;j++)
			{
				if(a[j][i]==1)
					for(int k=1;k<=5;k++)
						a[j][k]=(a[j][k]||a[i][k]);
			}
			for(int i=1;i<=5;i++)
			{
				for(int j=1;j<=5;j++)
					cout<<a[i][j]<<" ";
				cout<<endl;
			}