链接:https://www.luogu.org/problemnew/show/P2865

题目描述

Bessie has moved to a small farm and sometimes enjoys returning to visit one of her best friends. She does not want to get to her old home too quickly, because she likes the scenery along the way. She has decided to take the second-shortest rather than the shortest path. She knows there must be some second-shortest path.

The countryside consists of R (1 ≤ R ≤ 100,000) bidirectional roads, each linking two of the N (1 ≤ N ≤ 5000) intersections, conveniently numbered 1..N. Bessie starts at intersection 1, and her friend (the destination) is at intersection N.

The second-shortest path may share roads with any of the shortest paths, and it may backtrack i.e., use the same road or intersection more than once. The second-shortest path is the shortest path whose length is longer than the shortest path(s) (i.e., if two or more shortest paths exist, the second-shortest path is the one whose length is longer than those but no longer than any other path).

贝茜把家搬到了一个小农场,但她常常回到FJ的农场去拜访她的朋友。贝茜很喜欢路边的风景,不想那么快地结束她的旅途,于是她每次回农场,都会选择第二短的路径,而不象我们所习惯的那样,选择最短路。 贝茜所在的乡村有R(1<=R<=100,000)条双向道路,每条路都联结了所有的N(1<=N<=5000)个农场中的某两个。贝茜居住在农场1,她的朋友们居住在农场N(即贝茜每次旅行的目的地)。 贝茜选择的第二短的路径中,可以包含任何一条在最短路中出现的道路,并且,一条路可以重复走多次。当然咯,第二短路的长度必须严格大于最短路(可能有多条)的长度,但它的长度必须不大于所有除最短路外的路径的长度。

输入输出格式

输入格式:

Line 1: Two space-separated integers: N and R

Lines 2..R+1: Each line contains three space-separated integers: A, B, and D that describe a road that connects intersections A and B and has length D (1 ≤ D ≤ 5000)

输出格式:

Line 1: The length of the second shortest path between node 1 and node N

题解:

在做本题前,一定理解最短路算法的核心思想。 当然,并不需要对所有的最短路算法有深刻了解。 所以今天我就给大家带来优先队列优化的dijsktra。(SPFA我不会)

那么问题来了,经典的dijkstra只能维护最短路,怎么去维护次短路呢?其实很简单,只要同时维护两个数组就行了。

那怎么更新呢?我们设dis1数组存最短路,dis2数组存最长路,当找到一条路,比dis1和dis2都短时,就先把目前次短路赋为目前最短路,也就是让dis2=dis1,再把目前最短路赋为找到的更短路。剩下的就是dijkstra的板子了。

当我们存边时,还是推荐大家用邻接表存边,只用写一个函数,跑的好像还比vector快,下面就为大家带来一个不怎么好看的dijkstra板:

#include<iostream>
#include<cstring>
#include<queue>
#include<algorithm>
#include<stack>
#include<vector>
#include<cstdio>
#define ll long long
#define inf 1e9+10
using namespace std;
const int maxn=1e5+;
struct mxr
{
int to,next,val;
}e[*maxn];
int n,cnt,m,vis[maxn],head[maxn],s;
int dis[maxn];
priority_queue<pair<ll,int> >q;
inline int read()
{
int s=,f=;char ch=getchar();
while(ch<''||ch>''){ch=getchar();
if(ch=='-')f=-;}
while(ch>=''&&ch<='')
s=(s<<)+(s<<)+(ch^),ch=getchar();
return s*f;
}
inline void add(int u,int v,int w)
{
e[cnt].to=v;
e[cnt].val=w;
e[cnt].next=head[u];
head[u]=cnt++;
return ;
}
inline void dij()
{
q.push(make_pair(,s));
dis[s]=;
while(!q.empty())
{
int now=q.top().second;
q.pop();
if(vis[now]) continue;
vis[now]=;
for(int i=head[now];i!=-;i=e[i].next)
{
int v=e[i].to;
if(dis[v]>dis[now]+e[i].val)
{
dis[v]=dis[now]+e[i].val;
q.push(make_pair(-dis[v],v));
}
}
}
return ;
}
int main()
{
memset(head,-,sizeof(head));
cnt=;
n=read(),m=read(),s=read();
for(int i=;i<=n;i++) dis[i]=inf;
for(int i=;i<=m;i++)
{
int a,b,w;
a=read(),b=read(),w=read();
add(a,b,w);
}
dij();
for(int i=;i<=n;i++) printf("%d ",dis[i]);
return ;
}

下面是本题代码:

#include<iostream>
#include<cstring>
#include<queue>
#include<algorithm>
#include<stack>
#include<vector>
#include<cstdio>
#define ll long long
#define inf 1e9+10
using namespace std;
const int maxn=1e5+;
struct mxr
{
int to,next,val;
}e[*maxn];
int n,cnt,m,head[maxn];
int dis1[maxn],dis2[maxn];
priority_queue<pair<ll,int> >q;
inline char nc()//fread,跑的快
{
static char buf[],*L=buf,*R=buf;
return L==R&&(R=(L=buf)+fread(buf,,,stdin),L==R)?EOF:*L++;
}
inline int read()
{
int ret=,f=;char ch=nc();
while (ch<''||ch>''){if (ch=='-') f=-;ch=nc();}
while (ch>=''&&ch<='') ret=ret*+ch-'',ch=nc();
return ret*f;
}
inline void write(int x)
{
if(x<) putchar('-'),x=-x;
if(x>) write(x/);
putchar(x%+'');
return ;
}
inline void add(int u,int v,int w)//邻接表存边,貌似比vector快
{
e[cnt].to=v;
e[cnt].val=w;
e[cnt].next=head[u];
head[u]=cnt++;
return ;
}
inline void dij()//优先队列优化dijkstra
{
for(register int i=;i<=n;++i) dis1[i]=inf,dis2[i]=inf;
dis1[]=;
q.push(make_pair(,));
while(!q.empty())
{
int now=q.top().second,dis=-q.top().first;
q.pop();
for(register int i=head[now];i!=-;i=e[i].next)
{
int v=e[i].to;
if(dis1[v]>dis+e[i].val)
dis2[v]=dis1[v],dis1[v]=dis+e[i].val,q.push(make_pair(-dis1[v],v));//找到更优,更新dis1与dis2
if(dis2[v]>dis+e[i].val&&dis1[v]<dis+e[i].val)
dis2[v]=dis+e[i].val,q.push(make_pair(-dis2[v],v));//找到比dis2优,比dis1劣,更新dis2
}
}
return ;
}
int main()
{
memset(head,-,sizeof(head));
cnt=;
n=read(),m=read();
for(register int i=;i<=m;++i)
{
int a,b,w;
a=read(),b=read(),w=read();
add(a,b,w),add(b,a,w);
}
dij();
write(dis2[n]),puts(" ");
return ;
}

谢谢大家!

最新文章

  1. 七天学会ASP.NET MVC (一)——深入理解ASP.NET MVC
  2. jQuery-1.9.1源码分析系列(六) 延时对象应用——jQuery.ready
  3. Java 集合类 TreeSet、TreeMap
  4. Java绘图
  5. 0930MySQL中实现高性能高并发计数器方案(例如文章点击数)
  6. @RequestMapping详解
  7. Hibernate工作原理及为什么要用?
  8. JavaWeb学习记录(十九)——jsp标签库
  9. UIPickerView常见属性、常见方法(包括代理方法和数据源方法)的一些说明
  10. Bluetooth LE(低功耗蓝牙) - 第六部分(完)
  11. ASCII是指128个字符(不是256个)和ASCII Extended Characters(就是那些奇怪的外文字符)
  12. python2与python3
  13. 路由知识 静态路由 rip eigrp ospf
  14. 利用Azure嵌套虚拟化,解决公有云上机器不能启动的问题
  15. js万年历,麻雀虽小五脏俱全,由原生js编写
  16. Win7系统用户文件夹多出一个Administrator.xxx开头的文件怎么解决
  17. 【题解】魔板—洛谷P1275。
  18. java 异步查询转同步多种实现方式:循环等待,CountDownLatch,Spring EventListener,超时处理和空循环性能优化
  19. 潭州课堂25班:Ph201805201 爬虫高级 第四课 sclapy 框架 crawispider类 (课堂笔记)
  20. day9-复习学习python实例

热门文章

  1. 网络监控之一:netstat命令
  2. 1125 Chain the Ropes
  3. django-上传文件 fromdata(头像实例)
  4. Docker - 避免启动container后运行shell脚本执行完成后docker退出container
  5. Nginx启动提示找不到libpcre.so.1解决方法
  6. leetcode594
  7. ffmpeg截取一段视频中一段视频
  8. 【291】Python 中字符串添加到剪贴板
  9. json_encode和json_decode和isset和array_key_exists
  10. Log4php使用指南