次短路问题(一次最短路 or 两次相反方向最短路)

次短路

模板题poj3255

题意：

给一张n个点的路径图，求点1到n的次短路

解法：

1.两次最短路
以1为源点求出d1[],以n为源点求出d2[]
设次短路长度为ans，遍历所有边，设当前边为w[u,v]
d1[u]+w[u,v]+d2[v]等于1到n最短路，那么说明w[u,v]在最短路上，跳过.
否则更新ans=min(ans,d1[u]+w[u,v]+d2[v])
2.一次dijkstra+同时更新最短/次短路长度
d1[],d2[]依次为最短距离和次短距离
跑dj的时候同时更新d1[],d2[]
但是个人觉得代码细节太多了，不够无脑，不太喜欢
代码自带细节注释

两次最短路：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#define set0(a) memset(a,0,sizeof(a))
#define set1(a) memset(a,-1,sizeof(a))
#define setinf(a) memset(a,inf,sizeof(a))
#define ll long long
const int inf=0x3f3f3f3f;
const int inn=0x80808080;
using namespace std;
const int maxm=2e5+5;
int head[maxm],to[maxm],nt[maxm],w[maxm];
int mark[maxm];
int d1[maxm],d2[maxm];
int cnt;
void add(int x,int y,int z){
    cnt++;nt[cnt]=head[x];head[x]=cnt;to[cnt]=y,w[cnt]=z;
    cnt++;nt[cnt]=head[y];head[y]=cnt;to[cnt]=x,w[cnt]=z;
}
void init(){
    set1(head);
    setinf(d1);
    setinf(d2);
    cnt=0;
}
void spfa(int st,int *d){
    queue<int>q;
    q.push(st);
    set0(mark);
    mark[st]=1;
    d[st]=0;
    while(!q.empty()){
        int x=q.front();;
        q.pop();
        mark[x]=0;
        for(int i=head[x];i!=-1;i=nt[i]){
            int v=to[i];
            if(d[v]>d[x]+w[i]){
                d[v]=d[x]+w[i];
                if(!mark[v]){
                    mark[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
int main(){
    init();
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    spfa(1,d1);
    spfa(n,d2);
    int ans=inf;
    for(int x=1;x<=n;x++){
        for(int i=head[x];i!=-1;i=nt[i]){
            int v=to[i];
            int temp=d1[x]+d2[v]+w[i];
            if(temp>d1[n]&&temp<ans){
                ans=temp;
            }
        }
    }
    printf("%d\n",ans);
    return 0;
}

dijkstra：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#define set0(a) memset(a,0,sizeof(a))
#define set1(a) memset(a,-1,sizeof(a))
#define setinf(a) memset(a,inf,sizeof(a))
#define ll long long
const int inf=0x3f3f3f3f;
const int inn=0x80808080;
using namespace std;
const int maxm=2e5+5;
int head[maxm],to[maxm],nt[maxm],w[maxm];
int d1[maxm],d2[maxm];
int cnt;
void add(int x,int y,int z){
    cnt++;nt[cnt]=head[x];head[x]=cnt;to[cnt]=y,w[cnt]=z;
    cnt++;nt[cnt]=head[y];head[y]=cnt;to[cnt]=x,w[cnt]=z;
}
void init(){
    set1(head);
    setinf(d1);
    setinf(d2);
    cnt=0;
}
void dj(int st){
    priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > >q;
    d1[st]=0;//注意:这个地方只能初始化d1[st]不能初始化d2[st]
    q.push(make_pair(0,st));
    while(!q.empty()){
        int d=q.top().first;
        int x=q.top().second;
        q.pop();
        if(d>d2[x])continue;//如果距离比次短路还长那么肯定无法更新
        for(int i=head[x];i!=-1;i=nt[i]){
            int v=to[i];
            int temp=d+w[i];
            if(temp<d1[v]){
                swap(d1[v],temp);//这个地方一定是交换，用旧的最短路更新次短路
                q.push(make_pair(d1[v],v));
            }
            if(d1[v]<temp&&d2[v]>temp){
                d2[v]=temp;
                q.push(make_pair(d2[v],v));
            }
        }
    }
}
int main(){
    init();
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    dj(1);
    printf("%d\n",d2[n]);
    return 0;
}