最短路算法总结(Dijkstra)

如何求任意两点之间的最短路径呢?在之前的学习里,知道可以通过深搜或者广搜求出两点之间的最短路径。但学习了Dijkstra这个新的算法以后,会更方便。

Dijkstra算法的基本思想:
每次找到离源点最近的一个顶点,然后以该顶点为中心进行扩展,最终得到源点到其余所有点的最短路径。

基本步骤:
1. 将所有顶点分为两部分:已知最短路程的顶点集合P和未知最短路径的顶点集合Q。最开始,已知最短路径的顶点集合P中只有源点一个顶点。我们这里用一个book数组来记录哪些顶点再集合P中。例如对于某个顶点i,如果book[i]为1则表示这个顶点再集合P中,如果book[i]为0则表示这个顶点再集合Q中。
2. 设置源点x到自己的最短路径为0即 dis[x]=0。若存在有源点能直接到达的顶点i,则把dis[i]设为g[x][i]。同时把所有其他(源点不能直接到达的)顶点的最短路径设为∞。
3. 在集合Q的所有顶点中选择一个离源点x最近的顶点u(即dis[u]最小)加入到集合P。并考察所有以点u为起点的边,对每一条边进行松弛操作。例如存在一条从u到v的边,那么可以通过将u->v添加到尾部来拓展一条从s到v的路径,这条路径的长度时dis[u]+g[u][v]。如果这个值比目前已知的dis[v]的值要小,我们可以用新值来替代当前dis[v]中的值。
4. 重复第3步,如果集合Q为空,算法结束。最终dis数组中的值就是源点到所有顶点的最短路径。

其实它的过程是很简单的,下面是代码的实现:

Floyd-Warshall:

#include<stdio.h>
int main()
{
	int i,j,k,m,n,t1,t2,t3,e[10][10];
	int inf = 99999999;
	scanf("%d%d",&n,&m);
	for(i = 1; i <= n; i ++)
		for(j = 1; j <= n; j ++)
			if(i == j)
				e[i][j] = 0;
			else
				e[i][j] = inf;
	for(i = 1; i <= m; i ++)
	{
		scanf("%d%d%d",&t1,&t2,&t3);
		e[t1][t2] = t3;//双向e[t1][t2] = e[t2][t1] = t3;
	}
	//Floyd-Warshall核心语句 
	for(k = 1; k <= n; k ++)
		for(i = 1; i <= n; i ++)
			for(j = 1; j <= n; j ++)
				if(e[i][j] > e[i][k] + e[k][j])
					e[i][j] = e[i][k] + e[k][j];
	for(i = 1; i <= n; i ++)
	{
		for(j = 1; j <= n; j ++)
			printf("%d",e[i][j]);
		printf("\n");
	}
	return 0;
}

Dijkstra:

#include<stdio.h>
int m,n,e[10][10],dis[10],book[10];
int inf = 99999999;
void Dijkstra(int x)
{
	int i,j,u,min;
	for(i = 1; i <= n; i ++)
	{
		dis[i] = e[x][i];
		book[i] = 0;
	}
	dis[x] = 0;
	book[x] = 1;
	for(i = 1; i < n; i ++)
	{
		min = inf;
		for(j = 1; j <= n; j ++)
			if(book[j] == 0 && dis[j] < min)
			{
				min = dis[j];
				u = j;
			}	
		book[u] = 1;
		for(j = 1; j <= n; j ++)
			if(book[j] == 0 && dis[j] > dis[u] + e[u][j])
				dis[j] = dis[u] + e[u][j];
	}
}
int main()
{
	int i,j,t1,t2,t3;
	scanf("%d%d",&n,&m);
	for(i = 1; i <= n; i ++)
		for(j = 1; j <= n; j ++)
			if(i == j)
				e[i][j] = 0;
			else
				e[i][j] = inf;
	for(i = 1; i <= m; i ++)
	{
		scanf("%d%d%d",&t1,&t2,&t3);
		e[t1][t2] = t3;//双向e[t1][t2] = e[t2][t1] = t3
	}
	Dijkstra(1);
	return 0; 
}

另外,还有一个负权边的情况,就是在给出的距离上会出现负的,此时就不能再用上面的代码来实现了,因为扩展到负权边的时候会产生更短的路程,有可能就破坏了已经更新的点路程不会改变的性质。

此时,用Bellman-Ford算法就能很好的解决了。

代码实现:

#include<stdio.h>
int main()
{
	int i,j,k,m,n,count,flag;
	int dis[10],u[10],v[10],w[10];
	int inf = 99999999;
	scanf("%d%d",&n,&m);
	for(i = 1; i <= m; i ++)
		scanf("%d%d%d",&u[i],&v[i],&w[i]);
	for(i = 1; i <= n; i ++)
		dis[i] = inf;
	dis[1] = 0;
	//Bellman-Ford核心语句 
	for(k = 1; k < n; k ++)
	{
		count = 0;
		for(i = 1; i <= m; i ++)
			if(dis[v[i]] > dis[u[i]] + w[i])
			{
				dis[v[i]] = dis[u[i]] + w[i];
				count = 1;
			}
		if(count == 1)
			break;
	}
	flag = 0;
	for(i  =1; i <= m; i ++)
		if(dis[v[i]] > dis[u[i]] + w[i])
			flag = 1;
	if(flag == 1)
		printf("此图含有负权回路");
	return 0; 
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值