最短路的几种算法及其优化(模板)

一.Dijkstra 算法

    dijkstra算法适用于边权为正的情况,求单源最短路,适用于有向图和无向图

    模板伪代码:

                  清除所有点的标号

                  设d[0]=0,其余d[i]=INF;

                  循环n次{

                             在所有未标记的节点中,寻找d[i]最小的点x

                             给x做标记

                            对于从x出发的所有边(x,y)更新d[y]=min(d[y],d[x]+w[x,y]);

                             }

                  

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 j=0;j<n;++j)if(!visit[j]&&d[j]<m){m=d[j];x=j;}visit[x]=1;for(int j=0;j<n;++j)d[j]=min(d[j],d[x]+w[x][j]);} 

     简单说一下dijkstra的优化:

              1.储存结构上:邻接矩阵是很占空间的(众所周知),所以我们一般采用邻接表或者边表

              2.堆优化:因为在dijkstra中有循环n次寻找最小dict的过程,我们可以维护一个小根堆来实现,也就把时间复杂度从n^2降到了n*(logn+E)。

     优化后的dijkstra,自己写的:

 

/*
建图用的邻接表,复杂度O(E*logE)
*/struct pnode {int num;int len;pnode() {}pnode(int a, int b) : num(a), len(b) {}//初始化结构体用的,把a复制给num,把b复制给len; bool operator < (const pnode tmp) const {return len > tmp.len;}
};int dis[N];
bool vis[N];
int n;void dijkstra(int s) {priority_queue<pnode> q;q.push(pnode(s, 0));pnode u;int v, i, res = inf;for(i = 0; i <= n; ++i) dis[i] = inf, vis[i] = false;dis[s] = 0;while(!q.empty()) {u = q.top(); q.pop();if(u.len != dis[u.num]) continue;/*这是应对优先队列中的重复入队的点,只要最新的那个点就可以了*/if(vis[u.num])  continue;vis[u.num] = true;for(i = head[u.num]; i != -1; i = g[i].next) {v = g[i].to;if(dis[v] > u.len + g[i].val) {dis[v] = u.len + g[i].val;q.push(pnode(v, dis[v]));}}}
}

二.Bellman-Ford的优化(也就是SPFA,直接看三吧)

三.SPFA模板及SPFA的优化

   1.普通SPFA模板(队列化的Bellman-Ford算法):

int visit[N],dis[N];
bool SPFA(int s)
{queue<int>q;memset(dis,127,sizeof(dis));memset(visit,false,sizeof(visit));
memset(cnt,0s,sizeof(cnt));dis[s]
=0;visit[s]=true;q.push(s);while(!q.empty()){int k=q.front();q.pop();visit[k]=false;for(int i=head[k];i;i=edge[i].last)/*边表*/{if(dis[k]+edge[i].w<dis[edge[i].v]){dis[edge[i].v]=dis[k]+edge[i].w;if(!visit[edge[i].v]){q.push(edge[i].v);visit[edge[i].v]=true;
if(++cnt[edge[i].v]>n) /*如果某一个点的入队次数超过了n次,说明存在负环,返回false*/
return false;}}}}return true;/*安全跑完了,不存在环*/ }

    2.SPFA的优化

      SPFA算法有两个优化策略SLF和LLL——SLF:Small Label First 策略,设要加入的节点是j,队首元素为i,若dist(j)<dist(i),则将j插入队首,否则插入队尾; LLL:Large Label Last 策略,设队首元素为i,队列中所有dist值的平均值为x,若dist(i)>x则将i插入到队尾,查找下一元素,直到找到某一i使得dist(i)<=x,则将i出队进行松弛操作。 SLF 可使速度提高 15 ~ 20%;SLF + LLL 可提高约 50%。 在实际的应用中SPFA的算法时间效率不是很稳定,为了避免最坏情况的出现,通常使用效率更加稳定的Dijkstra算法。实际上dijkstra算法+heap优化后是一定快于一般SPFA的,而且更加稳定。

1)SPFA的SLF优化,(很简单的,只要使用双端队列就可以实现了)。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <deque>
using namespace std;
const int N=501;
const int NN=100001;
const int inf=0x7fffffff;
int n,nu;
typedef struct node
{int adj,val;struct node *next;
};
node node[NN],*p[N];
int SPFA()
{deque<int> qu;int x,i,a,b;int vis[N],dis[N],num[N];struct node *head[N];for(i=1;i<=n;i++){vis[i]=0;num[i]=0;dis[i]=inf;head[i]=p[i];}dis[1]=0;vis[1]=1;num[1]++;qu.push_back(1);while(!qu.empty()){x=qu.front();qu.pop_front();vis[x]=0;head[x]=p[x];while(head[x]){a=head[x]->adj;b=head[x]->val;if(dis[a]>dis[x]+b){dis[a]=dis[x]+b;if(!vis[a]){vis[a]=1;num[a]++;if(num[a]>=n)return 1;if(!qu.empty()){if(dis[a]>dis[qu.front()])qu.push_back(a);elsequ.push_front(a);}elsequ.push_back(a);}}head[x]=head[x]->next;}}return 0;
}
int main()
{int t,i,m,w,a,b,c;scanf("%d",&t);while(t--){memset(node,0,sizeof(node));memset(p,0,sizeof(p));nu=0;scanf("%d%d%d",&n,&m,&w);for(i=0;i<m;i++){scanf("%d%d%d",&a,&b,&c);node[nu].adj=b;node[nu].val=c;node[nu].next=p[a];p[a]=&node[nu];nu++;node[nu].adj=a;node[nu].val=c;node[nu].next=p[b];p[b]=&node[nu];nu++;}for(i=0;i<w;i++){scanf("%d%d%d",&a,&b,&c);node[nu].adj=b;node[nu].val=-c;node[nu].next=p[a];p[a]=&node[nu];nu++;}if(SPFA())puts("YES");elseputs("NO");}return 0;
}
网上找的
#include<iostream>
using namespace std;
#include<deque>
#include<cstdio>
#include<cstring>
int n,m;
#define N 1001
struct Edge{int u,v,w,last;
}edge[N];
int head[N];
int dict[N];
bool visit[N];
void input()
{scanf("%d%d",&n,&m);/*n个点,m条边*/for(int i=1;i<=m;++i){scanf("%d%d%d",&edge[i].u,&edge[i].v,&edge[i].w);edge[i].last=head[edge[i].u];head[edge[i].u]=i;}    
}
bool SPFA()
{deque<int>q;memset(dict,127,sizeof(dict));memset(visit,false,sizeof(visit));dict[1]=0;visit[1]=true;;int cnt[N];memset(cnt,0,sizeof(cnt));/*判断有无环的标志*/++cnt[1];while(!q.empty()){int k=q.front();q.pop_front();visit[k]=false;/*取出后,不要忘记标志位*/for(int l=head[k];l;l=edge[l].last)/*边表*/{int p=edge[l].v;if(dict[p]>dict[k]+edge[l].w){dict[p]=dict[k]+edge[l].w;++cnt[p];if(cnt[p]>n) return true;/*如果某个点的入队次数超过了n,那么一定存在环*/if(!visit[p]){visit[p]=true;if(!q.empty())/*这就是SLF Small Label First 策略.的核心,把将要入队的元素的dict与队首元素相比较,如果将要入队的元素的dict大的话,就放在队尾,否则就放在队首 */{if(dict[p]>dict[q.front()])/*这样可以保证始终用dict小的更新,也节约了时间*/q.push_back(p);else q.push_front(p);/*所以必须使用双端队列*/}else q.push_back(p);/*不要忘记考虑队列为空的情况*/}}}}return false;
}
int main()
{input();if(SPFA())printf("circle");else{for(int i=1;i<=n;++i)printf("%d ",dict[i]);}return 0;
}
自己打了一遍,还可以理解

2.SPFA的LLL优化

   在网上实在没找到可靠的。

3.SPFA的DFS优化:

   在很多的题目中,SPFA都是用BFS来实现的,对于随机图来说,BFS的速度会远大于DFS,但是对于某些特殊的图结构来说,DFS也是一个很好的选择

    例如 1):题目中要求在最短的时间内,判断有无环,DFS明显会比BFS快(例题是POj上一个判断单词接龙的题目)

          2):对于网格型的图,DFS的速度比BFS快

          模板:  

void SPFA(int k)
{flag[k]=true;for(int l=head[k];l;l=edge[l].last){int v=edge[l].v;/*找出第一个可以被更新的点*/if(dis[v]>dis[k]+edge[l].w){dis[v]=dis[k]+edge[l].w;if(!flag[v]){SPFA(v);/*接着深搜下去*/}else /*这表明从某个点开始DFS,结果又搜到了这个点,说明存在负权回路*/{printf("cycle");return;}}}flag[k]=false;/*不要忘记再把k设为false,为的是让他能够重复入队*/
}

 

四,Floyd算法模板(没有任何优化,就是邻接矩阵和n^3,一般情况单源最短路是绝对不会用的)

for(int k=1;k<=n;++k)/*注意这个k必须在最外层循环才行*/for(int i=1;i<=n;++i)for(int j=1;j<=n;++j)dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]); 

 

  

转载于:https://www.cnblogs.com/c1299401227/p/5401240.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/357308.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

极客青年说,北京沙龙

大家好&#xff0c;我是孙叫兽&#xff0c;本期内容给大家分享infoq写作平台在北京站的沙龙主题。主要流程如下&#xff1a;早上起来&#xff0c;吃了点饭&#xff0c;然后就去北京望京悠乐汇的A1202,这次是个轰趴馆&#xff0c;比较适合年轻的程序员放松。去的稍微早了点&…

html css 知识回顾2

在学习代码期间&#xff0c;练习代码是必不可少的内容&#xff0c;同样的对知识点的回顾也是更加重要的&#xff0c;在练代码的同时回顾知识点效果会更好。&#xff08;知识会有重复的&#xff0c;我遵循的是不论是什么时候的知识都会从头来一遍&#xff0c;对于新掌握的我会再…

vue点击弹窗自动触发点击事件的解决办法

业务场景&#xff1a;使用vue element ui 的el-dialog&#xff0c;点击弹窗之后&#xff0c;默认加载第一个按钮的数据进行初始化。 div 指令&#xff1a; // 自动触发点击事件directives:{trigger:{inserted(el,binging){// console.log("自动触发事件")el.click()}…

谈谈yii2-gii如何自定义模板

作者&#xff1a;白狼 出处&#xff1a;http://www.manks.top/article/yii2_gii_custom_template本文版权归作者&#xff0c;欢迎转载&#xff0c;但未经作者同意必须保留此段声明&#xff0c;且在文章页面明显位置给出原文连接&#xff0c;否则保留追究法律责任的权利。 yii2中…

vue自定义指令directives实现自动点击事件及自动点击第一个按钮

业务场景:点击弹窗默认加载第一个按钮的数据。vue自定义指令directives实现这个需求 目录 自动点击所有的按钮。 自动点击第一个按钮, 自动点击所有的按钮。 <ul class="areaList"><li:class="{active: index

人员雇佣 网络流_雇用Java EE开发人员的一些面试问题

人员雇佣 网络流互联网上充斥着Java开发人员的面试问题。 这些问题的主要问题是&#xff0c;它们仅证明候选人具有良好的记忆力&#xff0c;并记住所有语法&#xff0c;结构&#xff0c;常量等。对他/她的逻辑推理没有真正的评估。 我在下面列举一些面试问题的例子&#xff0c;…

vue根据表格字段不同的状态显示不同的颜色。

业务需求:根据后台返回的数据,对表格中的严重等级和问题状态做一下颜色区分。数据很大,大概一年左右的数据,在二级弹窗中,数据滚动的形式。 大家好,我是孙叫兽 不加状态前: 我这个使用vue+div循环的实现,很便捷,使用element ui也是可以的。 然后再computed:{}中添加…

[IT业界] 网盘纷纷停止服务,下一个是谁?

今天新浪微盘发公告称即日起调整服务功能&#xff0c;停止普通用户的存储服务&#xff0c;也就是将关闭新浪微盘的搜索、分享功能。2016年以来&#xff0c;已经有115网盘、UC网盘等多家网络数据在线存储产品宣布关闭文件存储、分享功能&#xff0c;这对这些产品本身而言是致命的…

Apache JMeter教程

要负载测试您的Web应用程序吗&#xff1f; 然后&#xff0c;您应该了解Apache JMeter &#xff01; 该JMeter教程介绍了基本概念以及如何创建测试计划以及如何执行它以对应用程序进行基准测试。 这是我在Disy的Tech-Blog上的文章的交叉张贴 -谢谢您让我在公司时间写文章&…

浏览器在线PDF预览取消下载按钮

业务场景:需要在线PDF下载的功能,刚开始用的是第三方的PDF.JS,按照网上的方法操作,排查,发现没有作用,经过和后端小伙伴沟通,原来这个地址是浏览器PDF预览。 目录 浏览器PDF取消办法: 在需要URL地址后边加上代码 PDF.JS取消下载的办法:</

Android深度探索(卷1)HAL与驱动开发 读书笔记(第四章)

第四章 源代码的下载和编译 本章主要介绍使用Git下载两套源代码。一套是Android 源代码&#xff0c;另一套是Linux 内核源代码。主要介绍如何下载和编译Android源代码和Linux内核源代码。 4.1下载Android源代码的环境 第一步&#xff1a;创建一个存放下载脚本文件repo的目录 #…

【Java从入门到天黑|04】JavaSE入门之数组

目录 数组概述 数组的四个基本特点: 数组声明创建 1、声明数组 2、创建数组

stage3d 骨骼优化

用过Away3D的朋友估计都会发现&#xff0c;&#xff0c;在Away3D里面使用超过一定骨骼数量的角色&#xff0c;当场景里面角色的数量稍微多一点&#xff0c;整个场景就会很卡。 对于这个现象&#xff0c;我之前得出的结论是。Stage3D的VC缓存器数量的限制&#xff0c;造成了对需…

morphia_Morphia和MongoDB:不断发展的文档结构

morphia在我先前关于Morphia的帖子中 &#xff0c;我介绍了一些典型用法&#xff0c;并提到了一些已知问题的警告和解决方法。 我展示了使用Morphia的工作有多么容易&#xff0c;以及它与Java世界的交互方式有多么干净。 为了跟进该帖子的后续内容&#xff0c;我将讨论如何处理…

PC端网页vue项目,页面滚动点击悬浮按钮最快的方法

业务场景,可视化大屏右下角添加一个悬浮按钮,点击按钮页面滚动到顶部。 经过实验,使用锚点是最简单的办法。 方法一:锚点 给顶部的合适位置添加一个 name="top"的锚点,maodianTop用于控制样式,可以不设置。 <a class="maodianTop" name="…

Echats给柱状图及提示文字添加百分号(%)的解决办法

业务需求&#xff1a;给柱状图添加百分号。 目录 柱状图上添加%效果图 悬浮标签添加%效果图&#xff1a; 刚开始从后台取的数据带%&#xff0c;我这边取这个值的时候显示undifined&#xff0c; 后来就让这哥们把这类的数据从数据库把%去掉。这样我这边就取到了数据&#xff0…

前端使用linux命令更新项目生产包与测试包命令

业务需求&#xff1a;把vue开发的项目打成dist.zip文件&#xff0c;丢到服务器去&#xff0c;通过域名进行访问。 首先登录云管平台的账号和密码&#xff0c;找到对应的服务器&#xff0c;然后连接到堡垒机。 在命令行进入到域名下的地址及目录。 切换到测试文件夹 点击左上角…

servle 3.0 新特性之一 对上传表单的支持

1. 上传 * 上传对表单的要求&#xff1a; > method"post" > enctype"multipart/form-data"&#xff0c;它的默认值是&#xff1a;application/x-www-form-urlencoded > <input type"file" name"必须给"/> * 上传Servle…

manifest.mf._MANIFEST.MF和feature.xml版本控制规则

manifest.mf.我永远都忘记了OSIF插件和功能的 MANIFEST.MF和feature.xml中的依赖项声明的规则是什么。 谷歌搜索经常导致沮丧而不是答案。 因此&#xff0c;因为今天我实际上找到了这些规则的简要列表&#xff0c;所以我想在这里重新发布它们&#xff0c;并进行一些较小的修改以…

CSDN学习神器——CSDN浏览器助手测评体验

导读&#xff1a;大家好&#xff0c;我叫孙叫兽&#xff0c;本期内容给大家分享一下CSDN浏览器助手这款插件&#xff0c;看看最新版的插件已经比较完善了&#xff0c;还有一些可以优化的地方。下面简单进行体验一把最新版&#xff01; 体验时间&#xff1a;2021年5月31日。 目录…