题目描述
小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
/*
本题涉及到 倍增 + 最短路 ,所以还是可以用
floyd方法求最短路,只不过此处k不是边的限制了,
而是最原始经典floyd的点。
注意点:
1.要有floyd()经典数组 d[][]
2.要有判断能不能到达的数组: (bool)f[36][N][N] (f[i][u][v])
表示 u走 2^i 秒 能不能到达j
3.注意了:
不能只在 输入m条边的时候,初始化f[][][],d[][]
而还要应该在 枚举:阶段i,发现那些可以跳跃到达的
点也初始化 f[][][] ,d[][]
(因为你跳跃着从一个点u到达另外一个点v,边长d[u][v]也可以看成1)
4.由于本题是maxlongint: INT_MAX,为2^32,
所以不能用之前的18
*/
#include <bits/stdc++.h>
using namespace std;
const int N=100;
int d[N][N];
bool f[40][N][N];
int main()
{
int n,m;
cin>>n>>m;
memset(d,0x3f,sizeof d);
for(int i=1;i<=n;i++)
d[i][i]=0;
while(m--)
{
int u,v;
cin>>u>>v;
d[u][v]=min(1,d[u][v]);
f[0][u][v]=true; //走一步就能到达
}
for(int i=1;i<=32;i++)
for(int u=1;u<=n;u++)
for(int v=1;v<=n;v++)
for(int k=1;k<=n;k++)
{
if(f[i-1][u][k]&&f[i-1][k][v]) //如果两半都能到达:一个是u-->k,一个是k-->v
{
f[i][u][v]=true;
d[u][v]=1;
}
}
for(int k=1;k<=n;k++)
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]<<endl;
return 0;
}