原以为次短路算法会很麻烦,结果是随便乱搞(雾)。
次短路算法可以拓宽图论时我们的思维和乱搞能力(雾),看一道引题。
【引题】用一张图和一些边,求一列点之间的最短通路,要求初始点和终点必须都是集合中的点。
如果不能从原图中有效地抽离出那些点的话,直接用floyd算法,复杂度就是O(n3)。但是如果建立一个超级源,从源点出发,向每一个点连一条边长为0的边,再求到每点之间的次短路就可以很好地解决这个问题。(当然集合中的元素必须远小于n,不然从原题中抽离出的思维过程就是无意义的)。
【算法】
从源点和汇点之间各跑一遍最短路,再枚举每条边,枚举出"路径长是大于最短路且最小的路径"就是次短路。
限制:源点和汇点都是必须确定的。
时间复杂度:O(2最短路算法+2m)=O(最短路算法)。
【code】
最短路算法推荐还是dijkstra,这里为了编者方便 还是展示最短路为spfa的算法
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define int long long
using namespace std;
const int maxm=2e5+100;
const int maxn=5100;
struct node{
int v,nxt,c;
}edge[maxm];
int head[maxn],tot;
int n,m;
int dist[maxn],dist1[maxn];
bool inq[maxn];
inline void read(int &x){
x=0;char tmp=getchar();
while(tmp<'0'||tmp>'9') tmp=getchar();
while(tmp>='0'&&tmp<='9') x=(x<<1)+(x<<3)+tmp-'0',tmp=getchar();
}
inline void add_edge(int x,int y,int c){
edge[tot].v=y;
edge[tot].c=c;
edge[tot].nxt=head[x];
head[x]=tot++;
}
signed main(){
memset(head,-1,sizeof(head));
cin>>n>>m;
for(int i=1;i<=m;i++){
int x,y,c;
read(x),read(y),read(c);
add_edge(x,y,c),add_edge(y,x,c);
}
memset(dist,0x3f,sizeof(dist));
dist[1]=0;
queue<int>q;q.push(1);
while(!q.empty()){
int u=q.front();q.pop();inq[u]=0;
for(int i=head[u];i!=-1;i=edge[i].nxt){
int v=edge[i].v,c=edge[i].c;
if(dist[v]>dist[u]+c){
dist[v]=dist[u]+c;
if(!inq[v]) inq[v]=1,q.push(v);
}
}
}
memset(dist1,0x3f,sizeof(dist1));
dist1[n]=0; q.push(n),inq[n]=1;
while(!q.empty()){
int u=q.front();q.pop();inq[u]=0;
for(int i=head[u];i!=-1;i=edge[i].nxt){
int v=edge[i].v,c=edge[i].c;
if(dist1[v]>dist1[u]+c){
dist1[v]=dist1[u]+c;
if(!inq[v]) inq[v]=1,q.push(v);
}
}
}
int ans=1LL<<62;
for(int i=1;i<=n;i++)
for(int j=head[i];j!=-1;j=edge[j].nxt){
if(dist1[edge[j].v]+edge[j].c+dist[i]>dist[n])
ans=min(ans,dist1[edge[j].v]+edge[j].c+dist[i]);
}
cout<<ans<<endl;
return 0;
}