洛谷：P1613 跑路（倍增 + Floyd维护最短路）

洛谷：跑大路

一开始着手题意有点难懂，理解后可以发现是很有意思的一道题

给定一张图，边权相同都为正（可以视为1），而且观察数据可以发现此图有可能有环存在。

然后就是 跑路机 的定义：

跑路鸡 跑的距离是固定的，一次（一秒）只能跑 $2^k（2的整数倍）$的距离

如果裸求最短路，即使 1 到 n 的最短路为3，用 跑路鸡 也要跑两秒；如果有一条路径长度为4，显然 跑路鸡 跑一秒就到了，更快！

因此：

我们用倍增的思想（毕竟关系到$2^k$）去维护 某点 i 到 某点 j 之间是否能 一秒到达（用 跑路鸡 的前提下）

这样肯定要遍历所有路径，dp维护，倍增本质还是dp嘛

看到数据 $n<=50$ ，显然 Floyd 算法很适合此处的 dp 维护

在 Floyd 维护好 倍增路径 后，再在新图的倍增路径上跑一次最短路，这个最短路就能满足题目的 跑路鸡 求解了

代码：

#include<bits/stdc++.h>
#include<unordered_set>
#include<unordered_map>
#define mem(a,b) memset(a,b,sizeof a)
#define cinios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
#define sca scanf
#define pri printf
#define ul (u << 1)
#define ur (u << 1 | 1)
#define fx first
#define fy second
//#pragma GCC optimize(2)
//[博客地址](https://blog.csdn.net/weixin_51797626?t=1) 
using namespace std;

typedef long long ll;
typedef pair<int, int> PII;

const int N = 60, M = 110, MM = 3000010;
int INF = 0x3f3f3f3f, mod = 100003;
ll LNF = 0x3f3f3f3f3f3f3f3f;
int n, m, k, T, S, D;
bool g[N][N][65];//布尔g[i][j][k]记录 i点 到 j点 是否有一条长度为 2^k(千米) 的路径
int dist[N][N];//最终 i 到 j 的距离（单位千米）

void floyd() {
	//题目数据范围在long int之内
	for (int d = 1; d <= 64; d++)//Floyd框架下的 dp
		for (int k = 1; k <= n; k++)
			for (int i = 1; i <= n; i++)
				for (int j = 1; j <= n; j++)
					if (g[i][k][d - 1] && g[k][j][d - 1]) {

						//如果 i ~ k 、k ~ j 都有一条长度为 2^(d-1) 的路径
						//显然 i ~ j 就有一条长度为 2^d 的路径
						g[i][j][d] = true;//标记
						dist[i][j] = 1;//更新 i ~ j 的跑路距离，显然一步就能到
					}

	for (int k = 1; k <= n; k++)
		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= n; j++)
				//最后再跑一遍 floyd 维护 跑路鸡 状态下的最短路即可
				dist[i][j] = min(dist[i][j], dist[i][k] + dist[k][j]);
}

int main() {
	cinios;

	cin >> n >> m;
	mem(dist, 0x3f);//初始化

	while (m--)
	{
		int a, b;
		cin >> a >> b;//a ~ b 有 一千米 长度的边，也就是 跑路鸡 跑 2^0(千米) 
		g[a][b][0] = true;
		dist[a][b] = 1;
	}
	
	floyd();
	
	cout << dist[1][n];

	return 0;
}