【题目链接】
ybt 1498：Roadblocks
洛谷 P2865 [USACO06NOV] Roadblocks G
【题目考点】
1. 图论：严格次短路径
严格次短路的路径长度必须大于最短路的路径长度。
非严格次短路的路径长度大于等于最短路的路径长度。
【解题思路】
每个交叉路口是一个顶点，每条路是无向边，求从顶点1到顶点n的严格次短路。
使用Dijkstra堆优化算法。
设Path类，包含属性u和d，表示存在一条从源点出发到达顶点u，长度为d的路径。
设优先队列pq，优先队列中保存的元素为Path类型对象，路径长度d更小的Path对象更优先。
设dis1，dis2数组，dis1[i]表示从源点到顶点i的最短路径长度，dis2[i]表示从源点到顶点i的次短路径长度。
首先将dis1和dis2数组每个元素都设为无穷大。
已知源点1到顶点1自己的最短路径长度为0，设dis1[1]=0，那么存在一条从源点到顶点1，长为0的路径，将对象Path{1, 0}入队到优先队列pq。
每次循环从优先队列出队长度最短的路径，取出该路径为从源点到达顶点u，路径长度为d。
访问顶点u的每个邻接点，顶点u到其邻接点v的边权为w。那么就存在一条从源点到顶点v的长为d+w的路径。
- 如果该到达顶点v的长为d+w的路径比从源点到顶点v的最短路径长度
dis1[v]更小，那么原来的最短路径长度变为次短路径长度，即dis2[v] = dis1[v]，当前的最短路径长度是d+w，设dis1[v] = d+w。
现在存在一条新的到达顶点v长为dis1[v]的路径，将对象Path{v, dis1[v]}入队。 - 如果该到达顶点v的长为d+w的路径比从源点到顶点v的最短路径长度
dis1[v]更大（注意不能等于dis1[v]），d+w比次短路径长度dis2[v]更小，那么当前的次短路径长度应该为d+w，设dis2[v] = d+w。
现在存在一条新的到达顶点v长为dis2[v]的路径，将对象Path{v, dis2[v]}入队。
最后结果为源点1到顶点n的严格次短路长度，即dis2[n]。
关于使用Dijkstra堆优化算法求次短路时，不能设vis数组进行优化
如果使用Dijkstra堆优化算法求最短路径，每个顶点只出队1次即可，第2次出队时没有必要继续扩展。可以设vis数组记录顶点是否已出队。而在求次短路过程中不能使用该方法进行优化。
已知从优先队列中出队的各个Path对象的路径长度d属性的单调递增的。
如果出队的Path对象为到达顶点u路径长度为d1。根据优先队列的比较规则，此时d1小于等于优先队列中所有Path对象的d属性。
通过顶点u扩展出的新的路径的长度为d1+w（w为顶点u到其某个邻接点的边权），因为Dijkstra的前提是图中没有负权边，所以d1+w一定大于d1，因此新入队的Path对象的d属性也大于d1。
因此接下来出队的Path对象的d属性一定大于等于d1，即按出队顺序看，Path对象的d属性是单调递增的。
顶点u第一次出队时，出队的Path对象为到达顶点u有长为d1的路径。顶点u第二次出队，出队的Path对象为到达顶点u有长为d2的路径。那么根据上述原理，一定有 d 2 ≥ d 1 d2\ge d1 d2≥d1。认为到顶点u有长为d1的路径，接下来扩展得到到其它顶点的最短路径长度，一定小于等于认为到顶点u有长为d2的路径，接下来扩展得到到其它顶点的最短路径长度。因此如果一个顶点第二次出队，就没有必要再继续进行扩展了。
vis[u]表示顶点u是否已经出队。如果顶点u已经出队过了，则不再访问更新其邻接点。否则设vis[u]为真，标记顶点u已出队，接下来访问更新其邻接点。
而求次短路时不应设vis数组记录顶点是否出队，因为当顶点u第二次出队时，如果是到顶点u有长为d2的路径，基于该存在的路径虽然不可能再更新各顶点的最短路径dis1的长度，但可能更新各顶点的次短路径dis2的长度。因此求次短路时不能进行该优化过程。
【题解代码】
解法1：Dijkstra堆优化算法
cpp
#include<bits/stdc++.h>
using namespace std;
#define N 5005
struct Edge
{
int v, w;
};
struct Pair
{
int u, d;//u：顶点 d：sv到u有一条长为d的路径
bool operator < (const Pair &b) const
{
return b.d < d;
}
};
vector<Edge> edge[N];
int n, m, dis1[N], dis2[N];//dis1[i]:源点到i的最短路径长度 dis2[i]：源点到i的严格次短路长度
void dijkstra(int sv)
{
priority_queue<Pair> pq;
memset(dis1, 0x3f, sizeof(dis1));
memset(dis2, 0x3f, sizeof(dis2));
dis1[sv] = 0;
pq.push(Pair{sv, dis1[sv]});
while(!pq.empty())
{
int u = pq.top().u, d = pq.top().d;
pq.pop();
for(Edge e : edge[u])
{
int v = e.v, w = e.w;
if(dis1[v] > d+w)
{
dis2[v] = dis1[v];
dis1[v] = d+w;
pq.push(Pair{v, dis1[v]});
}
else if(dis1[v] < d+w && d+w < dis2[v])
{
dis2[v] = d+w;
pq.push(Pair{v, dis2[v]});
}
}
}
}
int main()
{
int f, t, w;
cin >> n >> m;
for(int i = 1; i <= m; ++i)
{
cin >> f >> t >> w;
edge[f].push_back(Edge{t, w});
edge[t].push_back(Edge{f, w});
}
dijkstra(1);
cout << dis2[n];
return 0;
}