ACwing 853 - 有边数限制的最短路(Bellman - Ford)

给定一个n个点m条边的有向图,图中可能存在重边和自环, 边权可能为负数。

请你求出从1号点到n号点的最多经过k条边的最短距离,如果无法从1号点走到n号点,输出impossible。

注意:图中可能 存在负权回路 。

输入格式

第一行包含三个整数n,m,k。

接下来m行,每行包含三个整数x,y,z,表示存在一条从点x到点y的有向边,边长为z。

输出格式

输出一个整数,表示从1号点到n号点的最多经过k条边的最短距离。

如果不存在满足条件的路径,则输出“impossible”。

数据范围
1≤n,k≤500,
1≤m≤10000,
任意边长的绝对值不超过10000。

输入样例:

3 3 1
1 2 1
2 3 1
1 3 3

输出样例:

3

题目大意:

给出一个图,最多能走 k 条边,询问是否能到达 n 点, 如果可以到达则输出最短距离,如果不能到达则输出 “impossible” 。

解题思路:

既然存在负权边,则两种dijkstra都不能用了,考虑bellm - ford 和 spfa,但这道题有边数限制,只能走不超过 k 条边,询问能否到达n点,只能用bellm - ford 了。先看下bellm - ford 的核心算法

for (int i = 0; i < n; i ++)
    {
        for (int j = 1; j <= m; j ++)
        {
            int a = e[j].a, b = e[j].b, w = e[j].w;
            dis[b] = min(dis[b], dis[a] + w);
        }
    }

我们发现有内外两重循环,从作用上来看,外循环是走n条边(其实最多走n - 1 条边就够了,因为n 个点,最远n - 1条边),内循环则是枚举所有边进行松弛,既然只能走 k 次,那么改一下外循环,改成

for (int i = 0; i < k;  i++)

即可,因为我们最多走k条边,即最多走n 条边,所以改外循环就可以,内循环也要稍加修改,这里要借助一个存档数组backup,存dis上一次的状态,每次用上一次的状态去更新,如果不借助backup则会串联,举个例子:

  • 1 2 1
  • 2 3 1
  • 1 3 4

n = 3, k = 1,很好想,最后dis[n] 应该是 4 ,但是如果不借助backup数组存档,在第一次跑内循环的时候 dis[2] = min(dis[2], dis[1] + 1) = 1 这里是正确的,再往下走,dis[3] = min(dis[3], dis[2] + 1) = 2,这里就出错了,因为只能走一条边,这里却走了两条,所以我们应该用上一次 dis 的状态去更新,借助一个 backup 数组即可。

Code:

#include <iostream>
#include <cstring>

using namespace std;

const int N = 1e4 + 50;

int n, m, k;
int dis[N], backup[N];

struct edge { int a, b, w; } e[N];

void bellman()
{
    memset(dis, 0x3f, sizeof dis);
    dis[1] = 0;
    
    for (int i = 0; i < k; i ++)//最多经过 k 条边
    {
        memcpy(backup, dis, sizeof dis);//每次用上一次的状态去更新
        for (int j = 1; j <= m; j ++)
        {
            int a = e[j].a, b = e[j].b, w = e[j].w;
            dis[b] = min(dis[b], backup[a] + w);
        }
    }
}

int main()
{
    cin >> n >> m >> k;
    
    for (int i = 1; i <= m; i ++)
        cin >> e[i].a >> e[i].b >> e[i].w;
    
    bellman();
    
    if (dis[n] > 0x3f3f3f3f / 2) puts("impossible");
    else cout << dis[n] << endl;
    
    return 0;
}
posted @ 2020-08-22 18:25  Hayasaka  阅读(63)  评论(0编辑  收藏  举报