warshall算法求n阶邻接方阵的传递闭包

前言

离散课上老师提到了这一个高效求关系的传递闭包的算法，对于二进制狗的我需要多掌握一些算法，于是在此做个总结。

过程

算法很简单，但是语言描述起来有点费解。简单说，用i表示行，j表示列。

循环
i=0->n（n是阶数）
如果a[j][i]=1，则将j行的每个元素与i行进行逻辑加，即a[j][k]=a[i][k] | a[j][k]（k=0->n）

优化

时间上：
1.一开始邻接矩阵全是0，直接返回。
2.某歩计算之后邻接矩阵全是1，则可以在这步之后返回。
空间上：
将二维降为一维。

核心代码

void WarShall(bool* bM) {
	int n = iLevel * iLevel;

	if (iCZ == iLevel * iLevel)return;//全是孤立点，无所谓传递闭包

	for (int i = 0; i < iLevel; i++) {
		for (int j = 0; j < iLevel; j++) {
			if (/*bM[j][i]*/bM[j*iLevel+i]) {
				for (int k = 0; k < 4; k++) {
					if (!bM[j * iLevel + k] && (bM[j * iLevel + k] = bM[i * iLevel + k] | bM[j * iLevel + k])) {
						iCZ -= 1;
						if (!iCZ) {//如果全为1，则没有必要再算下去了
							return;
						}
					}		
				}
			}
		}
	}
}

运行结果