SPFA算法(SLF优化)2022.7.8更新
SPFA可能会被卡掉,能用dijkstra就别用SPFA,代码较长,但我已尽力做到解释,请耐心看下去,存储为邻接表存储。
#include<bits/stdc++.h>
#define inf 0x3f3f3f3f//(宏定义一个很大的值,例如0x3f3f3f3f等)
using namespace std;
int n,m,cnt;//cnt 计数器(有cnt条边)
struct edge//结构体定义
{
int v,w,nxt;//v 目标点 w 边权 nxt 这条边的上一条边(遍历)
};
edge e[3000010];//存边(边表)
int dis[3001000];//记录起点到第x个点的距离
int h[1001000];//记录第x个点所发出的最后一条边
bool f[1001000];//判断是否在队列内
deque<int> q;//双端队列
void add(int,int,int);
void spfa()//SLF优化
{
for(int i=1;i<=n;++i)
{
dis[i]=inf;
}//初始化,也可以用memset()
dis[1]=0;//起点到自己的距离为0 (该题起点为1)
q.push_back(1);//起点入队
f[1]=1;//标志起点已入队
while(!q.empty())
{
int top=q.front();//取队首元素
q.pop_front();//踢出队列
int w=dis[top];//w 起点的值
f[top]=0;//代表该元素已出队
for(int i=h[top];i;i=e[i].nxt)//邻接表遍历每一条边
{
int v=e[i].v;//目标点
int di=e[i].w;//边权
if(dis[v]>w+di)//松弛操作
{
dis[v]=w+di;//更新起点到点v的值
if(!f[v])//判断是否入队,没入队便入队
{
if(!q.empty()&&dis[v]<dis[q.front()])//如果比队首元素小,从队首入队
{//这里一定要把!q.empty()放在前面,编译时它会从前往后读,如果把它放后面而队列又为空,dis[v]<dis[q.front()]调用时会RE
q.push_front(v);//入队操作
}
else q.push_back(v);//否则从队尾入队
f[v]=1;//标志已入队
}
}
}
}
}
int main()
{
memset(h,0,sizeof h);//数组初始化
memset(f,false,sizeof f);//数组初始化
scanf("%d%d",&n,&m);
for(int a,b,c,i=1;i<=m;++i)
{
scanf("%d%d%d",&a,&b,&c);
add(a,b,c); //加边
add(b,a,c);//加边,如果是有向图,删掉这一步操作
}
spfa();//调用函数
printf("%d",dis[n]);//输出,根据题目要求,这里是输出1到n的最短距离
return 0;
}
void add(int u,int v,int w)
{
e[++cnt].v=v;//cnt 边的编号
e[cnt].w=w;
e[cnt].nxt=h[u];//指向上一条边
h[u]=cnt;//更新最后一条边
}
2022.7.8更新:
新增判断负环的写法和dfs写法
dfs写法判断负环更快,但求答案很慢
bfs队列写法求答案快,但判断负环很慢
根据题目情况来
bfs队列写法:
bool spfa()
{
for(rint i=1;i<=n;++i)
dis[i]=0x7ffffff;
dis[0]=0;
q.push_back(0);
vis[0]=true;
++in[0];//差记录入队次数
while(!q.empty())
{
int u=q.front();
q.pop_front();
vis[u]=false;
for(rint i=h[u];i;i=e[i].nxt)
{
int v=e[i].v;ll w=e[i].w;
if(dis[v]>dis[u]+w)
{
dis[v]=dis[u]+w;
if(!vis[v])
{
++in[v];
if(in[v]>=n+1) return false;//判断负环
if(!q.empty()&&dis[v]<dis[q.front()]) q.push_front(v);
else q.push_back(v);
vis[v]=true;
}
}
}
}
return true;
}
dfs写法:
void dfs_spfa(int u)
{
if(fg) return;
vis[u]=true;
for(rint i=h[u];i;i=e[i].nxt)
{
int v=e[i].v;
ll w=e[i].w;
if(dis[v]>dis[u]+w)
{
dis[v]=dis[u]+w;
if(vis[v]==true)//如果这个点被访问过,就说明这是负环
{
fg=true;//打标记
return;
}
else dfs_spfa(v);
}
}
vis[u]=false;
}
SPFA算法(SLF优化)2022.7.8更新的更多相关文章
- 关于SPFA算法的优化方式
关于SPFA算法的优化方式 这篇随笔讲解信息学奥林匹克竞赛中图论部分的求最短路算法SPFA的两种优化方式.学习这两种优化算法需要有SPFA朴素算法的学习经验.在本随笔中SPFA朴素算法的相关知识将不予 ...
- spfa的SLF优化
spfa的SLF优化就是small label first 优化,当加入一个新点v的时候如果此时的dis[v]比队首dis[q.front()]还要小的话,就把v点加入到队首,否则把他加入到队尾,因为 ...
- HDU4725(KB4-P SPFA+LLL+SLF优化)
The Shortest Path in Nya Graph Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K ...
- 并不对劲的图论专题(三):SPFA算法的优化
1.bzoj1489-> 这是个新套路. 我们希望找到最小的x,那么可以二分x,然后判断是否存在圈的边权的平均值小于等于x. 设圈的边权依次为w1,w2,w3,…,wk,平均值为p, 则有p= ...
- 《SPFA算法的优化及应用》——姜碧野(学习笔记)
一.核心性质:三角不等式.最短路满足d[v]<=d[u]+w(u,v) 二.SPFA两种实现: 常见的是基于bfs的,这是直接对bellman-ford用队列维护.根据最短路的长度最长为(n-1 ...
- 蓝书3.3 SPFA算法的优化
T1 最小圈 bzoj 1486 题目大意: 一个环的权值平均值为定义为一个这个环上所有边的权值和除以边数 求最小的环的权值平均值 思路: 二分一个值 把所有边减去这个值 判断是否有负环 #inclu ...
- 队列优化dijsktra(SPFA)的玄学优化
转载:大佬博客 最近想到了许多优化spfa的方法,这里想写个日报与大家探讨下 前置知识:spfa(不带任何优化) 由于使用较多 STLSTL ,本文中所有代码的评测均开启 O_2O2 优化 对一些数 ...
- SPFA算法的SLF优化 ——loj#10081. 「一本通 3.2 练习 7」道路和航线
今天做到一道最短路的题,原题https://loj.ac/problem/10081 题目大意为给一张有n个顶点的图,点与点之间有m1条道路,m2条航线,道路是双向的,且权值非负,而航线是单向的,权值 ...
- SPFA求最短路——Bellman-Ford算法的优化
SPFA 算法是 Bellman-Ford算法 的队列优化算法的别称,通常用于求含负权边的单源最短路径,以及判负权环.SPFA 最坏情况下复杂度和朴素 Bellman-Ford 相同,为 O(VE), ...
随机推荐
- 【Pandas vs SQL】数据分析代码逐行比对,孰优孰劣?
在数据分析领域,pandas是python数据分析基础工具,SQL是数据库最常用分析语言.二者有相通的地方,也有很大的语法不同,做起数据分析来,谁将更胜一筹呢? 做过业务开发.跟数据库打交道比较多的小 ...
- 在SpringBoot中使用logback优化异常堆栈的输出
一.背景 在我们在编写程序的过程中,无法保证自己的代码不抛出异常.当我们抛出异常的时候,通常会将整个异常堆栈的信息使用日志记录下来.通常一整个异常堆栈的信息是比较多的,而且存在一些没用的信息.那么我们 ...
- hooks 与 animejs
hooks 与 animejs 本文写于 2020 年 1 月 13 日 animejs 是现如今非常不错的一个 js 动画库.我们将其与 React Hooks 融合,使它更方便的在 React 中 ...
- 关于win10安装wsl子系统Ubuntu图形界面的错误解决
解决了https://blog.csdn.net/weixin_30834783/article/details/102144314Xserver个人使用的是VcXsrv. 在WSL中配置环境变量DI ...
- mysqldump还原备份数据时遇到一个问题
问题描述 ERROR 1839 (HY000) at line 24: @@GLOBAL.GTID_PURGED can only be set when @@GLOBAL.GTID_MODE = O ...
- 个人冲刺(四阶段)——体温上报app(一阶段)
任务:完成了后台数据库的类模块 MyDBHelper.java package com.example.helloworld; import android.content.Context; impo ...
- declare-声明限定类型变量
用于声明变量并设置变量的属性. 语法 declare [+/-][rxi][变量名称=设置值] declare -f 特殊符号 +/- "-"可用来指定变量的属性,"+& ...
- 如何让 Windows 把 TypeScript 文件当作文本文件
TL;DR 修改注册表项 HKEY_CLASSES_ROOT\.ts 为 HKEY_CLASSES_ROOT\.txt 的值 起因 Windows10 总把 TypeScript 文件自动当成视频,放 ...
- Java_Scanner的使用
目录 Scanner对象 scanner.next()和scanner.nextln()的区别 scanner.hasNext()和scanner.hasNextln() Scanner拓展 视频课程 ...
- 阿里巴巴开源限流组件Sentinel初探
1 Sentinel主页 https://github.com/alibaba/Sentinel/wiki/主页 1.1 Sentinel介绍 随着微服务的流行,服务和服务之间的稳定性变得越来越重要. ...