目录
一.单源最短路
1.dijkstra算法及实现
2.spfa算法及实现
(1)spafa负环判断及实现
二.多源最短路
1.floyd算法及实现
一.单源最短路
1.dijkstra算法及实现
求源点到图中其余各顶点的最短路径
dfs效率慢,解决规模小,bfs只能边权为1的图
Dijkstra算法——迪杰斯塔拉算法(非负全图)
基本思想:
首先假定源点为u,顶点集合V被划分为两部分:集合 S 和 V-S.
初始时S中仅含有源点u,其中S中的顶点到源点的最短路径已经确定。
集合S 和V - S中所包含的顶点到源点的最短路径的长度待定,
称从源点出发只经过S中的点到达V - S中的点的路径为特殊路径,
并用dist[]记录当前每个顶点对应的最短特殊路径长度。
实现从1到n的最短路径
输入:
3 3
1 2 5
2 3 5
3 1 2
输出:
2
#include<iostream>
#include<cstring>
using namespace std;
const int N = 1001;
const int M = 10001;
int head[N];
int size1;
int n, m;
struct edge {
int v, w, next;
edge() {};
edge(int _v, int _w, int _next) {
v = _v;
w = _w;
next = _next;
}
}e[M*2];
void init() {
memset(head, -1, sizeof(head));
size1 = 0;;
}
void insert(int u,int v,int w) {
e[size1] = edge(v, w, head[u]);
head[u] = size1++;
}
void insert2(int u, int v, int w) {
insert(u, v, w);
insert(v, u, w);
}
int dis[N];
bool vis[N];
void dijkstra(int u)
{
memset(vis, false, sizeof(vis));
memset(dis, 0x3f, sizeof(dis)); //0x3f用以表示正无穷
dis[u] = 0;
int mind = 1000000000;
for (int i = 0; i < n; i++) {
int minj = -1;
for (int j = 1; j <= n; j++) {
if (!vis[j] && dis[j] < mind) {
minj = j;
mind = dis[j];
}
}
if (minj == -1)
return;
vis[minj] = true;
for (int j = head[minj];~j ; j = e[j].next) { ///~j等价于j!=-1;
int v = e[j].v;
int w = e[j].w;
if (!vis[v] && dis[v] > dis[minj] + w)
dis[v] = dis[minj] + w;
}
}
}
int main()
{
init();
int u, v, w;
cin >> n >> m;
while (m--) {
cin >> u >> v >> w;
insert2(u, v, w);
}
dijkstra(1);
cout << dis[n] << endl;
return 0;
}
2.SPFA算法——Shortest Path Faster Algorithm
思路:
用数组dis记录每个结点的最短路径估计值,用邻接表或邻接矩阵来存储图G。
我们采取的方法是动态逼近法:设立一个先进先出的队列用来保存待优化的结点,
优化时每次取出队首结点u,
并且用u点当前的最短路径估计值对离开u点所指向的结点v进行松弛操作,
如果v点的最短路径估计值有所调整,且v点不在当前的队列中,就将v点放入队尾。
这样不断从队列中取出结点来进行松弛操作,直至队列空为止
优点:
(可以解决带负权的有向图并且在稀疏图优于dijsktra)
反例:
1.出题特殊数据使SPFA算法慢,2.有负环不能解决
推荐后续优化的dijkstra算法
实现:
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int N = 1001;
const int M = 10001;
int head[N];
int size1;
int n, m;
struct edge {
int v, w, next;
edge() {};
edge(int _v, int _w, int _next) {
v = _v;
w = _w;
next = _next;
}
}e[M*2];
void init() {
memset(head, -1, sizeof(head));
size1 = 0;;
}
void insert(int u,int v,int w) {
e[size1] = edge(v, w, head[u]);
head[u] = size1++;
}
void insert2(int u, int v, int w) {
insert(u, v, w);
insert(v, u, w);
}
int dis[N];
bool vis[N];
void spfa(int u) {
memset(vis, false, sizeof(vis));
vis[u] = true;
memset(dis, 0x3f, sizeof(dis));
dis[u] = 0;
queue<int>q;
q.push(u);
while (!q.empty()) {
u = q.front();
q.pop();
vis[u] = false;
for (int j = head[u]; ~j; j = e[j].next) {
int v = e[j].v;
int w = e[j].w;
if (dis[v] > dis[u] + w) {
dis[v] = dis[u] + w;
if (!vis[v]) {
q.push(v);
vis[v] = true;
}
}
}
}
}
int main()
{
init();
int u, v, w;
cin >> n >> m;
while (m--) {
cin >> u >> v >> w;
insert2(u, v, w);
}
spfa(1);
cout << dis[n] << endl;
return 0;
}
(1)spfa判断负环
在进行spfa用一个数组cnt标记每个顶点的入队次数,如果一个顶点的入队次数大于顶点总数
则表示该图含有负环
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int N = 1001;
const int M = 10001;
int head[N];
int size1;
int n, m;
struct edge {
int v, w, next;
edge() {};
edge(int _v, int _w, int _next) {
v = _v;
w = _w;
next = _next;
}
}e[M*2];
void init() {
memset(head, -1, sizeof(head));
size1 = 0;;
}
void insert(int u,int v,int w) {
e[size1] = edge(v, w, head[u]);
head[u] = size1++;
}
void insert2(int u, int v, int w) {
insert(u, v, w);
insert(v, u, w);
}
int dis[N], in[N];
bool vis[N];
bool spfa(int u) {
memset(vis, false, sizeof(vis));
vis[u] = true;
memset(dis, 0x3f, sizeof(dis));
dis[u] = 0;
memset(in, 0, sizeof(in));
in[u] = u;
queue<int>q;
q.push(u);
while (!q.empty()) {
u = q.front();
q.pop();
vis[u] = false;
for (int j = head[u]; ~j; j = e[j].next) {
int v = e[j].v;
int w = e[j].w;
if (dis[v] > dis[u] + w)
dis[v] = dis[u] + w;
if (!vis[v]) {
q.push(v);
vis[v] = false;
++in[v];
if (in[v] > n)
return true;
}
}
}
return false;
}
int main()
{
init();
int u, v, w;
cin >> n >> m;
while (m--) {
cin >> u >> v >> w;
insert2(u, v, w);
}
if (spfa(1))
cout << "Yes" << endl;
else
cout << "No" << endl;
return 0;
}
二.多源最短路
1.floyd算法——(Floyd-Warshall algorithm),中文亦称弗洛伊德算法
利用动态规划的思想解决带权图中任意两个点之间最短路径算法
优势:代码简短,高效,可以解决带负权文章来源:https://www.toymoban.com/news/detail-619878.html
反例:不能解决带负环
思路:
dp[0][i][j]=图
dp[k][i][j]=从i到j可能经过1...k最短路径,
min(dp[k][i][j],dp[k-1][i][k]+dp[k-1][k][j])
优化为二维:
min(dp[i][j],dp[i][k]+dp[k][j])
实现:
文章来源地址https://www.toymoban.com/news/detail-619878.html
#include<iostream>
#include<cstring>
using namespace std;
const int N = 101;
int g[N][N];
void floyd(int n) {
for (int k = 1; k <= n; k++) {
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= n; j++) {
g[i][j] = min(g[i][j], g[i][k] + g[j][k]);
}
}
}
}
int main()
{
memset(g, 0x3f, sizeof(g));
for (int i = 0; i < N; i++) {
g[i][i] = 0;
}
int n, m;
int u, v, w;
cin >> n >> m;
while (m--) {
cin >> u >> v >> w;
g[u][v] = g[v][u] = w;
}
floyd(n);
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= n; j++) {
cout << g[i][j] << ' ';
}
cout << endl;
}
return 0;
}
到了这里,关于求最短路径的三种算法的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!