分层图最短路学习笔记

一、分层图

分层图是指有很多个平行的图，各个平行的图之间有特殊的连接边。

用分层图的几种情况

• 有$k$个不同集合的边，将每个集合内的边建成一张图，再建立第$k+1$个图，是一个虚层，用这个虚层将这$k$张图连接起来。每次可以通过虚层转移到另一个集合的图中。
• 有$k$个机会使得走当前此边不花费代价或花费特殊的代价，可以建立$k$张相同的该图，每张图之间用有边的点连接起来，其代价是$0$或是特殊的值。每向下走一层，就代表用了一次机会，使得当前的路花费为$0$，最多可以走$k$次。如例2.飞行路线。
• 有$k$个机会逆向行驶，我们可以建$k$张相同的该图，将每层图之间有边的两个点用的逆向的边连接。每向下走一层，就代表用了一次机会逆向走了一次，最多可以走$k$次。

二、经典例题

例1. 小雨坐地铁

题意： 小雨所在的城市一共有$m$条地铁线，分别标号为$1$号线，$2$号线，$\dots$，$m$号线。整个城市一共有$n$个车站，编号为$1\sim n$。其中坐$i$号线需要花费$a_i$的价格，每坐一站就需要多花费$b_i$的价格。$i$号线有$c_i$个车站，而且这$c_i$个车站都已知，如果某一站有多条地铁线经过，则可以在这一站换乘到另一条地铁线，并且能多次换乘。现在小雨想从第$s$个车站坐地铁到第$t$个车站，地铁等待时间忽略不计，最最少花费的价格，若不能到达输出$-1$。
思路： 这就是建分层图的第一种情况，有$k$条地铁线路，我们就建$k$张图，然后再建第$k+1$张虚层，将各个可以中转线路的点连接起来，并设定从虚层到地铁需要乘该线的代价，从地铁到虚层代价为$0$，这里的虚层就相当于地铁站，这就实现了地铁的转线。
代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>

using namespace std;

typedef pair<int, int> pii;
const int N = 510010, inf = 0x3f3f3f3f;

int n, m, S, T;
int h[N], e[N], ne[N], w[N], idx;
int dist[N];
bool st[N];

void add(int a, int b, int c)
{
    e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
}

void dijkstra()
{
    memset(dist, 0x3f, sizeof(dist));
    dist[n * m + S] = 0;
    priority_queue<pii, vector<pii>, greater<pii> > heap;
    heap.push({0,n * m + S});
    while(heap.size()) {
        auto t = heap.top();
        heap.pop();
        int ver = t.second, distance = t.first;
        for(int i = h[ver]; ~i; i = ne[i]){
            int j = e[i];
            if(dist[j] > distance + w[i]) {
                dist[j] = distance + w[i];
                heap.push({dist[j], j});
            }
        }
    }
}

int main()
{
    int price, ad, cur, pre, num;
    memset(h, -1, sizeof(h));
    cin >> n >> m >> S >> T;
    for(int i = 1; i <= m; i ++) {
        cin >> price >> ad >> num;
        for(int j = 0; j < num; j ++) {
            cin >> cur;
            if(j) add(pre + (i - 1) * n, cur + (i - 1) * n, ad), add(cur + (i - 1) * n, pre + (i - 1) * n, ad);
            add(n * m + cur, cur + (i - 1) * n, price), add(cur + (i - 1) * n, n * m + cur, 0);
            pre = cur;
        }
    }
    dijkstra();
    int ans = dist[n * m + T];
    if(ans == inf) cout << -1 << endl;
    else cout << ans << endl;
    return 0;
}

例2. 飞行路线

题意： 求最短路，但是路线中有$k$条边的边权可以变为$0$
思路： 分层图套路，最重要的就是建图，建$k$张相同的该图，各层内部正常连边，各层之间建边时，建一条从上到下花费为$0$的边，代表免费乘坐一次。跑一遍从$s$到$n\ast k+t$的最短路即可。
代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>

using namespace std;

typedef pair<int,int> pii;
const int N = 10010, M = 100010, K = 15;

int n, m, k;
int S, T;
int h[N * K], e[2 * M * K], ne[2 * M * K], w[2 * M * K], idx;
int dist[N * K];
bool st[N * K];

void add(int a, int b, int c)
{
    e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
}

void dijkstra()
{
    for(int i = 1; i <= N * K; i ++) dist[i] = 0x3f3f3f3f;
    dist[S] = 0;
    priority_queue<pii, vector<pii>, greater<pii> > heap;
    heap.push({0, S});
    while(heap.size()) {
        auto t = heap.top();
        heap.pop();
        int ver = t.second, distance = t.first;
        if(st[ver]) continue;
        st[ver] = true;
        for(int i = h[ver]; ~i; i = ne[i]) {
            int j = e[i];
            if(dist[j] > distance + w[i]) {
                dist[j] = distance + w[i];
                heap.push({dist[j], j});
            }
        }
    }
}

int main()
{
    scanf("%d%d%d", &n, &m, &k);
    memset(h, -1, sizeof(h));
    scanf("%d%d",&S, &T);
    S ++, T ++;
    for(int i = 0; i < m; i ++) {
        int a, b, c;
        scanf("%d%d%d", &a, &b, &c);
        a ++, b ++;
        add(a, b, c), add(b, a, c);
        for(int j = 1; j <= k; j ++) {
            add(a + j * n, b + j * n, c), add(b + j * n, a + j * n, c);
            add(a + (j - 1) * n, b + j * n, 0), add(b + (j - 1) * n, a + j * n, 0);
        }
    }
    dijkstra();
    int ans = 0x3f3f3f3f;
    for(int i = 0; i <= k; i ++) ans = min(ans, dist[T + i * n]);
    printf("%d\n", ans);
    return 0;
}

转载并改编自大佬博客
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45735851/article/details/108219481