P1613 跑路 倍增，最短路Floyd算法

题目描述

小A的工作不仅繁琐，更有苛刻的规定，要求小A每天早上在6：00之前到达公司，否则这个月工资清零。可是小A偏偏又有赖床的坏毛病。于是为了保住自己的工资，小A买了一个十分牛B的空间跑路器，每秒钟可以跑2^k千米（k是任意自然数）。当然，这个机器是用longint存的，所以总跑路长度不能超过maxlongint千米。小A的家到公司的路可以看做一个有向图，小A家为点1，公司为点n，每条边长度均为一千米。小A想每天能醒地尽量晚，所以让你帮他算算，他最少需要几秒才能到公司。数据保证1到n至少有一条路径。

输入格式

第一行两个整数n，m，表示点的个数和边的个数。

接下来m行每行两个数字u，v，表示一条u到v的边。

输出格式

一行一个数字，表示到公司的最少秒数。

输入输出样例

输入 #1

4 4
1 1
1 2
2 3
3 4

输出 #1

1

说明/提示

【样例解释】

1->1->2->3->4，总路径长度为4千米，直接使用一次跑路器即可。

【数据范围】

50%的数据满足最优解路径长度<=1000；

100%的数据满足n<=50，m<=10000，最优解路径长度<=maxlongint。

思路：每秒可以跑2^k千米，等价于如果i到j的距离为2^k千米,那么只需1s即可到达。其次，如果i到t的距离为2的k-1次方，t到j的距离为2的k-1次方，那么i到j的距离仍满足为2的k次方，所以也只需1s.预处理完G和dis后，再跑一遍最短路，最后求出1到n需要多少秒。

代码如下：

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
int n,m,dis[60][60];
bool G[60][60][65];
void work(){//对G和dis预处理 
	for(int k=1;k<=64;k++)
	for(int i=1;i<=n;i++)
	for(int t=1;t<=n;t++)
	for(int j=1;j<=n;j++)
	if(G[i][t][k-1]&&G[t][j][k-1]){//如果满足i到t只需1秒，t到j只需1秒，距离都满足2的k-1次方，那么i到j满足距离为2的k次方，那么i到j也只需1s 
		G[i][j][k]=true;//表明i与j的距离为2^k 
		dis[i][j]=1;//只需1s 
	}
}
void Floyd(){//最短路 
	for(int k=1;k<=n;k++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
			}
		}
	}
}
int main(){
	memset(dis,INF,sizeof(dis));//初始化 
	memset(G,false,sizeof(G));
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		int u,v;
		cin>>u>>v;
		dis[u][v]=1;//初始化 ,u到v可以一秒到达 
		G[u][v][0]=true;//u到v的距离是1，是2的0次方，满足 
	}
	work();//初始化G[i][j][k] 
	Floyd();//最短路径入口 
	cout<<dis[1][n]<<endl;//输出1到n所需时间 
	return 0;
}