最短路相关模型

严格次短路

简洁明了，就是求图$G$中$s$到$t$的次短路，不过它在长度上要严格大于最短路。

我们维护两个数组，$dis{t}_i$表示$s$到$i$的最短距离，$dis{t}_i^{\prime }$表示$s$到$i$的严格次短距离，在图$G$上跑一遍$s$到$t$的$Dijksgtra$，每次更新距离时优先更新$dis{t}_i$，其次再更新$dis{t}_i^{\prime }$。

其实次短距离（不一定比最短距离长）做法十分相似，稍微改一下就行，详见代码。

代码

int n, m;
LL dist[N], dist2[N];
vector<Edge> e[N];

void Dijkstra(){
    rep(i, n + 1) dist[i] = dist2[i] = inf;
    priority_queue<pair<LL, int> > pq;
    dist[1] = 0;
    pq.push(MP(-dist[1], 1));
    while (!pq.empty()){
        int cur = pq.top().second;
        int curdis = -pq.top().first;
        pq.pop();
        if (curdis > dist2[cur]) continue;
        for (Edge u : e[cur]){
            if (curdis + u.len < dist[u.to]){
                dist2[u.to] = dist[u.to];
                dist[u.to] = curdis + u.len;
                pq.push(MP(-dist[u.to], u.to));
            } else if (curdis + u.len < dist2[u.to] && curdis + u.len > dist[u.to]){
                dist2[u.to] = curdis + u.len; // curdis + u.len > dist[u.to]表示严格次短
                pq.push(MP(-dist2[u.to], u.to));
            }
        }
    }
}

最短路中的差分约束

在有向图$G$中，我们记$s$到$i$的最短距离为$d_i$。若存在一条从$u$到$v$权值为$w$的边，则有$d_u+w\ge d_v$成立。

在满足所有这些不等式的$d$中，$d_t-d_s$的最大值就是$s$到$t$的最短距离。

考虑粗略地证明一下。首先，若$d_i$就是最短距离，则一定满足$d_u+w\ge d_v$，转化形式如下。
d v ≤ min ⁡ e = ( u , v ) ∈ E { d u + w } d_v \le \min_{e = (u, v) \in E}\{d_u + w\} dv​≤e=(u,v)∈Emin​{du​+w}
由于求最短路时$d$都初始化为$\inf$，如果一些$d_u$未在$d_v$之前求出，由于约束是取$min$，所以对$d_v$的最大值是没有影响的。因此我们可以按一定的顺序来逐个求解$d_i$的最大值。

若$s$到$i$的最短距离不是$d_i-d_s$的最大值呢？那么$\forall e=(u,v)\in E,d_u+w>d_v$，也就不存在路径了。

来看一道很好的例题。

P4878 [USACO05DEC]Layout G

题目描述

在$FarmerJohn$的农场上，有$N$头奶牛在排队等候喂食，编号为$1,2,\cdots ,N$，编号小的奶牛一定排在编号大的奶牛之前。

有$M_L$对奶牛希望彼此之间不要超过一定的距离，还有另外$M_D$对奶牛则希望相互之间不要小于一定的距离。

$FarmerJohn$想请你来求出$1$号奶牛和$N$号奶牛之间的最大距离。

输入格式

第一行三个整数$N,M_L,M_D$。

接下来$M_L$行，每行三个整数$A,B,D$，分别表示$A$号奶牛和$B$号奶牛之间的最大距离。

再接下来$M_D$行，每行三个整数$A,B,D$，分别表示$A$号奶牛和$B$号奶牛之间的最小距离。

输出格式

若无法排成一队，输出$-1$；若$1$号奶牛和$N$号奶牛之间的距离可以无穷大，输出$-2$；否则输出最大可能距离。

数据范围与提示

$2\le N\le 1000$

$1\le M_L,M_D\le 10^4$

$1\le A,B\le N$

$1\le D\le 10^6$

题解

记$1$号奶牛到$i$号奶牛的最大距离为$d_i$，我们可以上述约束转化为三类不等式。

1. 编号小的排在编号大的前，等价于$d_i<d_{i+1}\to d_{i+1}\ge d_i$（后面会解释为什么可以与等号）。
2. $M_L$对约束，$d_A+D\ge d_B$。
3. $M_D$对约束，$d_A+D\le d_B\to d_B-D\ge d_A$。

于是建图。

1. 从$i+1$向$i$连一条权值为$0$的边。
2. 从$A$向$B$连一条权值为$D$的边。
3. 从$B$向$A$连一条权值为$-D$的边。

从$1$到$N$跑一边$Bellman\_Ford$就行了，因为图中有负权边，所以不能跑$Dijkstra$。

坑

关于$-1$和$-2$的判断问题。$-1$代表图中存在负环，$-2$代表原图不连通。

来谈一谈“从$i+1$向$i$连一条权值为$0$的边”的精妙之处。

首先如果我们对$i$有了约束，那么对$1\sim i-1$就都有了约束，所以这并不影响图的连通性的判定。

另外，由于$D$的取值范围缘故，只要$A$号奶牛和$B$号奶牛之间有距离，$i(A<i<B)$号奶牛就可以找到自己的位置，所以对答案的贡献为$0$。

代码

#include <bits/stdc++.h>
#define PB push_back
#define MP make_pair
using namespace std;

#define rep(i, n) for (int i = 0; i < (int)(n); ++ i)
#define rep1(i, n) for (int i = 1; i < (int)(n); ++ i)
#define foreach(itr, c) for (__typeof((c).begin()) itr = (c).begin(); itr != (c).end(); ++ itr)

typedef long long LL;
typedef pair<int, int> pii;

const int N = 1005;
const LL inf = 1e16 + 7;

int n, ml, md;
LL dist[N];

struct Edge {
    int to;
    LL len;
};
vector<Edge> e[N];

LL Bellman_Ford(){
    rep(i, n + 1) dist[i] = inf;
    dist[1] = 0;
    rep(cnt, n){
        bool update = false;
        rep1(i, n + 1) for (Edge u : e[i]) if (dist[i] + u.len < dist[u.to]){
            dist[u.to] = dist[i] + u.len;
            update = true;
            if (cnt == n - 1) return -1;
        }
        if (!update) break;
    }
    if (dist[n] == inf) return -2;
    return dist[n];
}

int main(){
    scanf("%d%d%d", &n, &ml, &md);
    rep(i, ml){
        int u, v; LL w;
        scanf("%d%d%lld", &u, &v, &w);   // d_u + w >= d_v
        e[u].PB(Edge{v, w});
    }
    rep(i, md){
        int u, v; LL w;
        scanf("%d%d%lld", &u, &v, &w);   // d_u + w <= d_v  -->  d_v - w >= d_u
        e[v].PB(Edge{u, -w});
    }
    rep1(i, n) e[i + 1].PB(Edge{i, 0});
    LL ret = Bellman_Ford();
    printf("%lld\n", ret);
    return 0;
}