直接给出Dijkstra算法的伪代码,他可以算出从单个源点出发,到所有结点的最短路。该算法同时适用于有向图和无向图
清楚所有点的坐标
设d[0]=0,其他d[i]=INF
循环n次
{
在所有未标记结点中选出d值最小的结点x
给结点x标记
对于从x出发的所有边(x,y),更新d[y]=min(d[y], d[x]+w[x,y])
}
下面是伪代码对应的程序。假设起点是结点0,他到结点i的路径长度为d[i]。未标号结点的v[i]=0,已标号结点的v[i]=1。w[x][y]=-1表示边(x,y)不存在。
除了求出最短路的长度以外,使用Dijkstra算法也能方便的地打印出结点0到所有结点的最短路本身,方法就是从终点出发,不断顺着d[i]+w[i][j]=d[j]的边(i,j)从结点j退回到结点i,直到回到起点。另外也可以用空间换时间,在更新d数组时维护父亲结点。即
memset(v,0,sizeof(v));
for(int i = 0; i < n; i++) d[i]=(i==0 ? 0 : INF);
for(int i = 0; i < n; i++)
{
int x, m=INF;
for(int y = 0; y < n; y++) if(!v[y] && d[y] <= m) m=d[x=y];
v[x] = 1;
for(int y = 0; y < n; y++) if(w[x][y]!=-1 && d[y] > d[x]+w[x][y])
{
d[y] = d[x]+w[x][y];
fa[y] = x; //记录自己的父亲结点
}
}
稀疏图的邻接表:
不难看出,上面的程序的时间复杂度为O(n*n)。由于最短路算法实在太重要,我们有必要把他优化到O(mlongn),其中m是边的个数。在学习稀疏图的时候,首先要掌握图的新表示法。既然m远小于n*n,那么邻接矩阵中会有大量表示“此边不粗拿在”的元素,不仅浪费了空间,而且也降低了时间效率。为了解决这个问题,我们可以是用vector<int> g[maxn]存储稀疏图,但这里想介绍的是另外一种流行了表示法,邻接表。在这种表示方法中,每个结点i都有一个链表,里面存储着从i出发的所有边。对于无向图来说,每条边会在邻接表中出现两次。和前面一样,这里继续是用数组实现链表。首先对于每条边编号,然后用first[u]保存结点u的第一条边的编号,next[e]表示编号e的边的下一条边的编号。下面是构建图的代码。
int n, m;
int first[maxn];
int u[maxm],v[maxm],w[maxm],next[maxm]
void read(){
scanf("%d%d",&n,&m);
for(int i = 0; i < n; i++) first[i] = -1;
for(int i = 0; i < m; i++){
scanf("%d%d%d",&u[i],&v[i],&w[i]);
next[i] = first[u[i]];
first[u[i]] = i;
}
}
接下来是用优先队列加速Dijkstra算法。主要是在循环中获取最小的d[i]时,用优先队列来实现。