洛谷P1613 跑路

题目描述

小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的倍数，说明可以一秒直接到达，将这两个点的距离设为1。用bool数组表示两个点的状态 g[i,j,k] ：i到j要k步。三个点直接的距离为k-1，则可以合并为k。将所有能一步走到的点全部标记后，就可以用最短路解题

代码

#include<iostream>
#include<cstring>
#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef pair<ll,int> PII;
const int N=55,INF=1e16;
int h[N],ne[N],e[N],idx;
ll d[N][N];
bool g[N][N][66];
int n,m;
void add(int a,int b)
{
	e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
int main()
{
	memset(d,10,sizeof d);
	memset(h,-1,sizeof h);
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int a,b;
		cin>>a>>b;
		add(a,b);
		g[a][b][0]=true;  //  初始状态 
		d[a][b]=1;
	}
	
	for(int u=0;u<=64;u++)  // 枚举所有u  (2^u) 2^64次方足以 
	for(int i=1;i<=n;i++)  //  暴力枚举所有点 
	for(int k=1;k<=n;k++)
	for(int j=1;j<=n;j++)
	{
		if(g[i][k][u] && g[k][j][u])  //  如果这两个路径可以合并为一步走完，就合并 
		{
			g[i][j][u+1]=true;   //  维护这两个点的状态 
			d[i][j]=1;           //  更新距离，用于找最短路 
		}
	}
	
	for(int k=1;k<=n;k++){       //  floyd 找最短路 
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++)
			{
				d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
			}
		}
	}
	cout<<d[1][n];
	
	return 0;
}