Day18——图的连通性检测

谈“图”色变……
今天看到公式的时候呆了半天，一直不明白为什么能这样做，后来慢慢的才理出点头目来，看来很有必要来复习一下离散数学了。
今天判断图的连通性主要用的公式：
$$\begin{array}{rl}& M_a=M^0+M^1+M^2+\cdots +M^{n-1}\\ & 其中，M为图的邻接矩阵，n为顶点数\end{array}$$
图的连通性：
无向图：

• 如果图中任意两点是连通的，则称图是连通图。

有向图：

• 如果图中任意两个顶点v到顶点u有路径，且顶点u到顶点v有路径，则图是强连通图。

无向图可以看作是有向图的特例，在这里我们直接讨论有向图。
现在来解释上面的公式：
首先要明白图的邻接矩阵的含义，如：
$$M=\left[\begin{array}{ccc}0& 1& 0\\ 1& 0& 1\\ 0& 1& 0\end{array}\right]$$
$M_{ij}$的值表示顶点 i 是否能到达顶点 j ，若$M_{ij}=1$,表示能；否则不能。
$M^0$是一个单位阵，主对角线的值为1，表示每个顶点都能到达其自身。
$$\begin{array}{rl}& M^2=M\times M\\ & =\left[\begin{array}{ccc}0& 1& 0\\ 1& 0& 1\\ 0& 1& 0\end{array}\right]\times \left[\begin{array}{ccc}0& 1& 0\\ 1& 0& 1\\ 0& 1& 0\end{array}\right]\end{array}$$
把矩阵乘法拆开：
$$M_{ij}^2=M_{i1}\times M_{1j}+M_{i2}\times M_{2j}+M_{i3}\times M_{3j}$$
再结合矩阵的含义：
$M_{i}1表示i\to 1是否有路径，M_{1}j表示1\to j是否有路径$
两种情况：

• i 可到 1 ，1 可到 j，则i 可到 j，值大于0，表示i可以经过结点1到达j;
• i 到 j 没有路径

其它同理，求和后的$M_{ij}^2$的值表示 i 经过任意一个中转结点可到达 j 的条数。
可得，$M^l中的M_{ij}^l表示顶点i到j的长度为l的通路数$

所以$M_a=M^0+M^1+M^2+\cdots +M^{n-1}$是概括了所有情况，如果$M_a$中还存在值为0，则该图一定是非连通图。

package day18;

import day17.IntMatrix;

public class Graph {

	/**
	 * The connectivity matrix.
	 */
	IntMatrix connectivityMatrix;

	/**
	 * 
	 *********************
	 * The first constructor.
	 * 
	 * @param paraNumNodes The number of nodes in the graph.
	 *********************
	 *
	 */
	public Graph(int paraNumNodes) {
		connectivityMatrix = new IntMatrix(paraNumNodes, paraNumNodes);
	}// Of first constructor

	/**
	 * 
	 *********************
	 * The second constructor.
	 * 
	 * @param paraMatrix The data matrix.
	 *********************
	 *
	 */
	public Graph(int[][] paraMatrix) {
		connectivityMatrix = new IntMatrix(paraMatrix);
	}// Of second constructor

	/**
	 * Overrides the method claimed in Object,the superclass of any class.
	 */
	public String toString() {
		String resultString = "This is the connectivity matrix if the graph.\r\n" + connectivityMatrix;
		return resultString;
	}// Of toString

	/**
	 * 
	 *********************
	 * @Title: getConnectivity
	 * @Description: TODO(Get the connectivity of the graph.)
	 * 
	 * @throws Exception for internal error.
	 *********************
	 *
	 */
	public boolean getConnectivity() throws Exception {
		// Step 1.Initialize accumulated matrix:M_a=I.
		IntMatrix tempConnectivityMatrix = IntMatrix.getIdentityMatrix(connectivityMatrix.getData().length);

		// Step 2.Determine the actual connectivity
		IntMatrix tempMultipliedMatrix = new IntMatrix(connectivityMatrix);

		// Step 3.Determine the actual connectivity.
		for (int i = 0; i < connectivityMatrix.getData().length - 1; i++) {
			// M_a=M_a+M^k
			tempConnectivityMatrix.add(tempMultipliedMatrix);

			// M^k
			tempConnectivityMatrix = IntMatrix.multiply(tempConnectivityMatrix, connectivityMatrix);
		} // Of for i

		// Step 4.Check the connectivity.
		System.out.println("The connectivity matrix is: " + tempConnectivityMatrix);
		int[][] tempData = tempConnectivityMatrix.getData();
		for (int i = 0; i < tempData.length; i++) {
			for (int j = 0; j < tempData.length; j++) {
				if (tempData[i][j] == 0) {
					System.out.println("Node " + i + " cannot reach " + j);
					return false;
				} // Of if
			} // Of for j
		} // Of for i

		return true;
	}// Of getConnectivity

	/**
	 * 
	 *********************
	 * @Title: getConnectivityTest
	 * @Description: TODO(Unit test for getConnectivityTest)
	 *
	 *********************
	 *
	 */
	public static void getConnectivityTest() {
		// Test an undirected graph.
		int[][] tempMatrix = { { 0, 1, 0 }, { 1, 0, 1 }, { 0, 1, 0 } };
		Graph tempGraph2 = new Graph(tempMatrix);
		System.out.println(tempGraph2);

		boolean tempConnected = false;
		try {
			tempConnected = tempGraph2.getConnectivity();
		} catch (Exception e) {
			System.out.println(e);
		} // Of try.

		System.out.println("Is the graph connected? " + tempConnected);

		// Test a directed graph
		// Remove one arc to form a directed graph.
		tempGraph2.connectivityMatrix.setValue(1, 0, 0);

		tempConnected = false;
		try {
			tempConnected = tempGraph2.getConnectivity();
		} catch (Exception e) {
			System.out.println(e);
		} // Of try.
		System.out.println("Is the graph connected? " + tempConnected);
	}// Of getConnectivityTest

	/**
	 * 
	 *********************
	 * @Title: main
	 * @Description: TODO(The entrance of the program)
	 *
	 * @param args Not used now.
	 *********************
	 *
	 */
	public static void main(String args[]) {
		System.out.println("Hello!");
		Graph tempGraph = new Graph(3);
		System.out.println(tempGraph);

		// Unit test.
		getConnectivityTest();
	}// Of main
}// Graph

运行结果：

总结：
反推这个结论推了快一个小时，对于这块内容还是不敏感，推得我真的难受，去找个同学学习学习。