最短路径--------Floyd算法剖析

微博:http://weibo.com/375061590

QQ :375061590

 

用到两个重要矩阵:

 

         1.d[numVex][numVex]  (numVex图的顶点数):最开始该矩阵就是图的邻接矩阵,经过Floyd算法处理开后,d[numVex][numVex]中的d[i][j],表示着从顶点i到j的最短路径的权重。

 

         2.p[numVex][numVex]:p[i][j]表示从i到j的最短路径上 i的后继,例如1到5最短路劲为1-2-4-5  那么p[1][5]==2 ,最开始构建的p矩阵中p[i][j]= j

 

算法核心思想:  三圈for循环

 

for (int k = 0; k < graph.getNumVex(); k++) {

 

                     for (int v = 0; v < graph.getNumVex(); v++) {

 

                            for (int w = 0; w < graph.getNumVex(); w++) {

 

                                   if (d[v][w] > d[v][k] + d[k][w]) {

 

                                          d[v][w] = d[v][k] + d[k][w];

 

                                          p[v][w] = p[v][k];// p[v][w]是v--w最短路径上 v的下一顶点

 

                                   }

 

                            }

 

                     }

 

              }

 

第一层 k是作为中间顶点

 

第二层 v是作为起始顶点

 

第三层 w是作为终点顶点

 

 

 

内层核心代码

 

以v为起点,w为终点,再以k作为v和w之间的中间点,去判断d[v][ w]和d[v][k] + d[k][w]的大小关系,如果d[v][w] > d[v][k] + d[k][w],说明找到从v→w的更短路径了,此时更改d[v][w]的值为d[v][k] + d[k][w]。

 

 

 

p[v][w]的值也要相应改成p[v][k]的值,因为 p[v][k]的值是v→k最短路径上v的后继顶点,而v→w这段最短路径是连接在v→k这段路径后面的,所以令所当然p[v][w]也要指向p[v][k]。

 

注意:最外层的k循环,前面的n此循环的结果跟后面n+1次循环的错做过程是息息相关,

 

       三次循环完成后,各个顶点之间的最短路径权重会存储在d矩阵中:d[i][j]表示i→j的最短路径权重。

 

p中则存储着路径上的顶点,如果把i→j最短路径上的顶点列出来呢?

代码:

 

//求start→end 最短路径上的顶点

 

StringBuilder path = new StringBuilder();

 

int index = start;//起始点

 

path.append(start + " → ");

 

 

 

              while (index != end) {

 

              //循环取出路径上的各个顶点

 

                     index = p[index][end];

 

                     if(index != end){

 

path.append(index + " →");

 

}

 

                    

 

              }

 

用一个while循环循环 index = p[index][end];直到达到终点

 

假设该最短路径为 start→A→B→C→end

 

则执行过程是:

 

index = p[start][end]; →A

 

path== start→A →

 

 

 

index = p[A][end]; →B

 

path== start→A →B →

 

 

 

index = p[B][end]; →C

 

path== start→A →B →C→

 

 

 

index = p[C][end]; →end

 

path== start→A →B →C→end

 

 

 

这个就是p矩阵的妙处

 

 

 

 

 

 

 

测试:

 

图:



 

 

请输入定点的数目:5
顶点数为:5
请输入边数:7
边数为:7
请输入(Vi,Vj)上下标i 和  j,以及权重,用逗号隔开
0,1,5
0,4,7
1,2,4
4,2,8
1,3,2
2,3,6
4,3,1
初始的d矩阵
 
0 5 9999 9999 7
 
5 0 4 2 9999
 
9999 4 0 6 8
 
9999 2 6 0 1
 
7 9999 8 1 0
 
初始的p矩阵
 
0 1 2 3 4
 
0 1 2 3 4
 
0 1 2 3 4
 
0 1 2 3 4
 
0 1 2 3 4
 
处理后的d矩阵
 
0 5 9 7 7
 
5 0 4 2 3
 
9 4 0 6 7
 
7 2 6 0 1
 
7 3 7 1 0
 
处理后的p矩阵
 
0 1 1 1 4
 
0 1 2 3 3
 
1 1 2 3 3
 
1 1 2 3 4
 
0 3 3 3 4
 
求最短路径
请输入起点:
0
请输入终点:
2
从0到2的最短路径为9
该路劲为:0 → 1 →2
是否继续计算其他最短路径 Y/N?
y
求最短路径
请输入起点:
0
请输入终点:
3
从0到3的最短路径为7
该路劲为:0 → 1 →3
是否继续计算其他最短路径 Y/N?
y
求最短路径
请输入起点:
4
请输入终点:
1
从4到1的最短路径为3
该路劲为:4 → 3 →1
是否继续计算其他最短路径 Y/N?
y
求最短路径
请输入起点:
2
 
请输入终点:
4
从2到4的最短路径为7
该路劲为:2 → 3 →4
是否继续计算其他最短路径 Y/N?

 

源代码:

 

package DataStructure;

import java.util.Scanner;

public class Floyd {
	private Graph graph;
	private int[][] d;// 用来存储顶点到顶点之间最短路径的权重
	private int[][] p;// p[1][5]表示1到5的最短路径上 1的后继,例如1到5最短路劲为1-2-4-5 那么p[1][5]==2

	public Floyd() {
		this.graph = new Graph();
		d = graph.getArc();
		p = new int[graph.getNumVex()][graph.getNumVex()];
		initP();// 初始化矩阵p
		System.out.println("初始的d矩阵\n");
		for (int i = 0; i < graph.getNumVex(); i++) {
			for (int j = 0; j < graph.getNumVex(); j++) {
				System.out.print(d[i][j] + " ");
			}
			System.out.println("\n");
		}
		System.out.println("初始的p矩阵\n");
		for (int i = 0; i < graph.getNumVex(); i++) {
			for (int j = 0; j < graph.getNumVex(); j++) {
				System.out.print(p[i][j] + " ");
			}
			System.out.println("\n");

		}
		work();
		
		System.out.println("处理后的d矩阵\n");
		for (int i = 0; i < graph.getNumVex(); i++) {
			for (int j = 0; j < graph.getNumVex(); j++) {
				System.out.print(d[i][j] + " ");
			}
			System.out.println("\n");
		}
		
		System.out.println("处理后的p矩阵\n");
		for (int i = 0; i < graph.getNumVex(); i++) {
			for (int j = 0; j < graph.getNumVex(); j++) {
				System.out.print(p[i][j] + " ");
			}
			System.out.println("\n");
		}
	}

	/**
	 * 初始化p矩阵
	 * 
	 */
	private void initP() {
		for (int i = 0; i < graph.getNumVex(); i++) {
			for (int j = 0; j < graph.getNumVex(); j++) {
				p[i][j] = j;
			}
		}
	}

	/**
	 * 对d和p进行变化
	 * 
	 */
	private void work() {
		for (int k = 0; k < graph.getNumVex(); k++) {
			for (int v = 0; v < graph.getNumVex(); v++) {
				for (int w = 0; w < graph.getNumVex(); w++) {
					if (d[v][w] > d[v][k] + d[k][w]) {
						d[v][w] = d[v][k] + d[k][w];
						p[v][w] = p[v][k];// p[v][w]是v--w最短路径上 v的下一顶点
					}
				}
			}
		}
	}

	/**
	 * 获取最短路劲
	 * 
	 */
	public void getShortestPath(int start, int end) {
		StringBuilder path = new StringBuilder();
		int index = start;// 起始点
		path.append(start + " → ");

		while (index != end) {
			// 循环取出路径上的各个顶点
			index = p[index][end];
			if (index != end) {
				path.append(index + " →");
			}else {
				path.append(index);
			}

		}

		System.out.println("从" + (start) + "到" + (end) + "的最短路径为"
				+ d[start][end] + "\n该路劲为:" + path.toString());
	}

	public static void getShortestPath(Floyd floyd) {
		Scanner scanner = new Scanner(System.in);
		System.out.println("求最短路径\n请输入起点:");
		int start = scanner.nextInt();
		System.out.println("请输入终点:");
		int end = scanner.nextInt();
		floyd.getShortestPath(start, end);
		System.out.println("是否继续计算其他最短路径 Y/N? ");
		String tag = scanner.next();
		if (tag.toLowerCase().equals("y")) {
			getShortestPath(floyd);
		}

	}

	/**
	 * 图内部类
	 * 
	 * @author ccf
	 * 
	 */
	class Graph {
		/**
		 * 定点数
		 * 
		 */
		private int numVex = 0;
		private int arc[][] = null;
		private int numEdge = 0;
		private final int INFINITY = 9999;

		public Graph() {
			System.out.print("请输入定点的数目:");
			Scanner scanner = new Scanner(System.in);
			this.numVex = scanner.nextInt();
			arc = new int[numVex][numVex];
			for (int i = 0; i < numVex; i++) {
				for (int j = 0; j < numVex; j++) {
					arc[i][j] = INFINITY;
				}
			}
			for (int i = 0; i < numVex; i++) {
				arc[i][i] = 0;

			}
			System.out.println("顶点数为:" + this.numVex);
			System.out.print("请输入边数:");
			scanner = new Scanner(System.in);
			this.numEdge = scanner.nextInt();
			System.out.println("边数为:" + this.numEdge);

			System.out.println("请输入(Vi,Vj)上下标i 和  j,以及权重,用逗号隔开");
			for (int i = 1; i <= numEdge; i++) {
				scanner = new Scanner(System.in);
				String a = scanner.nextLine();
				String[] b = a.split(",");
				// System.out
				// .println("输入了:" + Integer.parseInt(b[0]) + " "
				// + Integer.parseInt(b[1]) + " "
				// + Integer.parseInt(b[2]));
				arc[Integer.parseInt(b[0])][Integer.parseInt(b[1])] = Integer
						.parseInt(b[2]);
				arc[Integer.parseInt(b[1])][Integer.parseInt(b[0])] = Integer
						.parseInt(b[2]);

			}
			
		}

		public int[][] getArc() {
			return arc;
		}

		public int getNumVex() {
			return numVex;
		}

	}

	public static void main(String[] args) {
		Floyd floyd = new Floyd();
		getShortestPath(floyd);
	}

}

 

 

 

你可能感兴趣的:(算法之道,算法,图,最短路径,Floyd)