数据结构学习之——最短路径（Dijkstra算法）

最短路径问题（Dijkstra算法）

先前的文章里介绍了两种图的遍历的方式，解决问题的能力也很有限，只能通过图向主人公反馈一个是否可以成功跳出水坑的结果，路径即使可以给出，或许也不是真正的最短路径。那要获取正真的最短路径，需要引入一种新的算法。

要解决的实例

先回顾一下先前的问题
例子是主人公站在水坑的中心点P1，水坑的长宽为5米，水坑散布着位置随机的砖块若干，主人公一次可以跳跃1米的长度来到达下一块砖，需要判断随机生成的点是否可以帮助主人公离开水坑。并且要求能给出最短的路径。

此图为随机生成的落脚点

程序设计思路

我们需要不断的确定距离中心点路径最短的点以及它距离中点的距离。
因此，对于每一个点，我们需要创建一个这样的结构：
参数1：distance，该点到中心点的最小路径
参数2：pre，该点到中点的最小路径经过的上一个点
参数3：seen，该点是否为我们已经处理过的点（处理的问题接下来讲）
初始化参数：我们需要将所有的结点以及中心点都放入到一个数组或者集合之中，方便我们对参数经行初始化。
参数初始化
1.需要将所有的distance都更新为无穷大，因为我们在后序计算中会根据路径的变化不断地跟新这个最小值，所以期望在预设时将其设定为无穷大。
2.pre可以设为空。
3.seen只是一个是否处理的标志位，可以设为-1。
定义好以上结构后我们需要创建我们用来收录这些结点的盒子，其实也可以不用，但是有可以帮助理解。
算法核心：
1.将原点放入盒子中，计算与原点直接相连的结点距离原点的距离写入distance。
2.选出所有节点中，除了原点和已经放入盒子的结点之外，与原点距离最近的结点。
3.计算该结点与其所有的未被收录的邻接点之间的距离dist1。并加上该结点的distance，若dist1+distance<邻接点的distance（这里体现为什么要将distance初始化为无穷大），则将该dist1+distance更新到该邻接点的distance。并将该邻接点的pre指向该结点。
重复经行2，3操作，直到所有结点都被收录到盒子中，则结束循环。

程序思路是这样，但是我们应当如何处理我们的实际问题。
在我们的实际问题中，因为结点的随机性，其实很多结点并不会被收录到盒子中，这些结点所在的连通分量可能并不包含中心点，因为我们限制了跳跃一步的距离只能是1m。所以有很多点可能根本无法涉足。
即使是这样也不妨碍我们利用上述思路解决问题。只是在选择邻接点时需要加上限制条件，距离需要小于1m。
最后可能还有一个问题就是我现在知道了中心点所在的连通分量中所有结点距离中心点的最短路径的距离以及这些结点的最短路径中上一个结点是哪一个。那么我如何确定哪一条是我可以跳出水坑的最短路径呢？
我们只需要分析盒子中的点中，哪些点可以跳出水坑，并且距离岸边的长度加上本身距离中心点的长度，就是该结点所在的最短路径的长度。只需要比较这些点找到其中最短的那个，便能得到最短路径以及它的长度了。

程序源码

//有权图的最短路径
void Dijkstra(void) {
	int i;
	int god;
	int Flag=-1;
	int min;
	struct Pos *pa;
	//初始化原点的参数
	P1->Stand->distance = 0;
	P1->Stand->seen = -1;
	P1->Stand->Pre = NULL;
	//寻找关于原点的邻接点
	for (i = 1; i <= FOOTHOLD; i++) {
		god = ((map->Position[i].x) - (P1->Stand->x))*((map->Position[i].x) - (P1->Stand->x)) + ((map->Position[i].y) - (P1->Stand->y))*((map->Position[i].y) - (P1->Stand->y));
		god = (int)sqrt(god);
		if (god < 10) {
			map->Position[i].Pre = P1->Stand;
			map->Position[i].distance = god;
			printf("[%d][%d] d:%d ", map->Position[i].x, map->Position[i].y, map->Position[i].distance);
		}
	}
	while (1) {
		min = 0xFFFF;
		Flag = -1;
		//寻找未收录的定点中dist最小的
		for (i = 1; i <= FOOTHOLD; i++) {
			//如果该点已被收录，则跳过该点
			if (map->Position[i].seen == 1) {
				continue;
			}
			if (map->Position[i].distance < min) {
				min = map->Position[i].distance;
				Flag = i;
			}
		}
		//若不存在未被收录的顶点，则跳出循环
		if (Flag == -1) {
			printf("Finished!");
			for (i = 1; i <= FOOTHOLD; i++) {
				min = 10;
				if (map->Position[i].seen == 1) {
					if (map->Position[i].x < min) {
						min = map->Position[i].x;
					}
					if ((50 - map->Position[i].x) < min) {
						min = 50 - map->Position[i].x;
					}
					if (map->Position[i].y < min) {
						min = map->Position[i].y;
					}
					if ((50 - map->Position[i].y) < min) {
						min = 50 - map->Position[i].y;
					}
					//若该点可以跳出水坑 打印路径
					if (min < 10) {
						map->Position[i].distance += min;
						printf("\nSucceed:[%d][%d] dist:%d\n", map->Position[i].x, map->Position[i].y, map->Position[i].distance);
						pa = &(map->Position[i]);
						while (1) {
							if (pa->x == 24&& pa->y == 24) {
								break;
							}
							printf("<-%d-[%d][%d]", pa->distance, pa->x, pa->y);
							pa = pa->Pre;
						}
					}
				}
			}
			break;
		}
		//将顶点设为已被访问（加入S[]）
		map->Position[Flag].seen = 1;
		//寻找该点所有的邻接点
		for (i = 1; i <= FOOTHOLD; i++) {
			//如果该点已被收录，则跳过该点
			if (map->Position[i].seen == 1) {
				continue;
			}
			//判断该点是否在邻接范围内
			god = ((map->Position[i].x) - (map->Position[Flag].x))*((map->Position[i].x) - (map->Position[Flag].x)) + ((map->Position[i].y) - (map->Position[Flag].y))*((map->Position[i].y) - (map->Position[Flag].y));
			god = (int)sqrt(god);
			if (god <= 10) {
				if ((god + map->Position[Flag].distance) < (map->Position[i].distance)) {
					map->Position[i].distance = god + map->Position[Flag].distance;
					map->Position[i].Pre = &(map->Position[Flag]);
				}
			}
		}
	}
}

程序有点长，思路基本上和上述以一致的。
最后挂个链接吧，这个链接我把程序更新了一些，原先广度搜索的程序可以实现搜索出需要跨步最少的路径。

链接：https://pan.baidu.com/s/1dzLrmdA0siojY55Z5S47Hw
提取码：grue