最小生成树Prim算法和单源最短路径Dijkstra算法

时间:2021-12-20 21:34:06

问题:

1. (最小生成树)给定一个带权的无向连通图,如何选取一棵生成树,使树上所有边上权的总和为最小,即求最小生成树。

2. (单源最短路径)给定一个权值都为正数的无向连通图和一个源点,确定它到其它点的最短距离。


之所以将这两个问题放在一起,是因为Prim算法与Dijkstra算法的思路和程序都非常相似,都是有贪心策略。


1.解法(Prim算法):

思路:设连通网络 N = { V, E },U表示已加入到生成树的顶点集合。在初始化阶段,任取一个顶点,用其关联边的权值作为初始的V-U顶点集合的数据。在每一步中,先在V-U顶点集合中选择距离U中任意点“最近”的顶点v,再把v加入到U中,最后看看在新的V-U顶点集合中,是否有哪个顶点距离v比距离U中其它点更近,若有则更新V-U顶点集合中的数据。U的元素个数为V-1,所以共要进行V-1步。


总的时间复杂度为Time = O(V)*T_EXTRACT-MIN+O(E)*T_DECREASE-KEY

若用数组作为顶点权值的数据结构,T_EXTRACT-MIN用时O(V),T_DECREASE-KEY用时O(1),总共用时O(V^2)

若用最小堆作为顶点权值的数据结构,T_EXTRACT-MIN用时O(lgV),T_DECREASE-KEY用时O(lgV),总共用时O((V+E)*lgV)

若用斐波那契堆作为顶点权值的数据结构,T_EXTRACT-MIN平均用时O(lgV),T_DECREASE-KEY平均用时O(1),总共用时O(E+V*lgV)

下面的代码使用数组作为顶点权值的数据结构:

#include <iostream>
#include <algorithm>
using namespace std;

#define MAXN 1001
#define INF 1000000

int lowcost[MAXN]; // 距离U中顶点的最小权值
bool visited[MAXN]; // 若为true表示已加入到集合U中
int mst[MAXN]; // 距离U中哪个顶点最短
int graph[MAXN][MAXN]; // 用矩阵表示图上各边权值

void Prim(int n)
{
int i, j;
memset(visited, 0, n*sizeof(bool));
visited[0] = true;
mst[0] = -1;
for (i=1; i<n; i++)
{
lowcost[i] = graph[0][i];
mst[i] = 0;
}
for (i=1; i<=n-1; i++)
{
// 取V-U中的最小权值T_EXTRACT-MIN O(V)
int min=INF, minid;
for (j=1; j<n; j++)
// 用visited数组确定V-U
if (!visited[j] && lowcost[j] < min)
{
min = lowcost[j];
minid = j;
}
visited[minid] = true;
// 减小V-U中的权值T_DECREASE-KEY O(1)
for (j=1; j<n; j++)
if (!visited[j] && lowcost[j] > graph[minid][j])
{
lowcost[j] = graph[minid][j];
mst[j] = minid;
}
}
}

int main()
{
int n, m, i, j;
cin >> n >> m;
for (i=0; i<n; i++)
for (j=0; j<n; j++)
graph[i][j] = INF;
for (i=0; i<m; i++)
{
char a[3], b[3];
int c;
scanf("%s %s %d",&a, &b, &c);
graph[a[0]-'A'][b[0]-'A'] = c;
graph[b[0]-'A'][a[0]-'A'] = c;
}
Prim(n);
int total = 0;
for (i=1; i<n; i++)
{
total += lowcost[i];
printf("%c-%c: %d\n", i+'A', mst[i]+'A', lowcost[i]);
}
printf("%d\n", total);
}


2.解法(Dijkstra算法):

思路:设连通网络 N = { V, E }和给定源点u,U表示已确定最短路径的顶点集合。在初始化阶段,用顶点u所关联边的权值作为初始的V-U顶点集合的数据。在每一步中,先在V-U顶点集合中选择距离源点u“最近”的顶点v,再把v加入到U中,最后看看在新的V-U顶点集合中,是否有哪个顶点通过顶点v到达源点u比通过U中其它点到达距离更短,若有则更新V-U顶点集合中的数据。U的元素个数为V-1,所以共要进行V-1步。


总的时间复杂度为Time = O(V)*T_EXTRACT-MIN+O(E)*T_DECREASE-KEY

若用数组作为顶点权值的数据结构,T_EXTRACT-MIN用时O(V),T_DECREASE-KEY用时O(1),总共用时O(V^2)

若用最小堆作为顶点权值的数据结构,T_EXTRACT-MIN用时O(lgV),T_DECREASE-KEY用时O(lgV),总共用时O((V+E)*lgV)

若用斐波那契堆作为顶点权值的数据结构,T_EXTRACT-MIN平均用时O(lgV),T_DECREASE-KEY平均用时O(1),总共用时O(E+V*lgV)

下面的代码使用数组作为顶点权值的数据结构:

#include <iostream>
#include <algorithm>

using namespace std;

#define MAXN 1001
#define INF 1000000

int lowcost[MAXN]; // 距离源点u的最短距离
bool visited[MAXN]; // 若为true表示已加入到集合U中
int minroad[MAXN]; // 在最短路径上顶点所连接的前一个顶点
int graph[MAXN][MAXN]; // 用矩阵表示图上各边权值

void Dijkstra(int n)
{
int i, j;
memset(visited, 0, n*sizeof(bool));
visited[0] = true;
minroad[0] = -1;
for (i=1; i<n; i++)
{
lowcost[i] = graph[0][i];
minroad[i] = 0;
}
for (i=1; i<=n-1; i++)
{
// 取V-U中的最小权值T_EXTRACT-MIN O(V)
int min=INF, minid;
for (j=1; j<n; j++)
// 用visited数组确定V-U
if (!visited[j] && lowcost[j] < min)
{
min = lowcost[j];
minid = j;
}
visited[minid] = true;
// 减小V-U中的权值T_DECREASE-KEY O(1)
for (j=1; j<n; j++)
if (!visited[j] && lowcost[j] > lowcost[minid]+graph[minid][j])
{
lowcost[j] = lowcost[minid]+graph[minid][j];
minroad[j] = minid;
}
}
}

int main()
{
int n, m, i, j;
cin >> n >> m;
for (i=0; i<n; i++)
for (j=0; j<n; j++)
graph[i][j] = INF;
for (i=0; i<m; i++)
{
char a[3], b[3];
int c;
scanf("%s %s %d",&a, &b, &c);
graph[a[0]-'A'][b[0]-'A'] = c;
graph[b[0]-'A'][a[0]-'A'] = c;
}
Dijkstra(n);
int total = 0;
for (i=1; i<n; i++)
{
total += lowcost[i];
printf("%c-%c: %d\n", i+'A', minroad[i]+'A', lowcost[i]);
}
printf("%d\n", total);
}

大家有没有觉得两个问题很像啊?哈哈~~