洛谷 P1613 跑路
Description
- 小A的工作不仅繁琐,更有苛刻的规定,要求小A每天早上在6:00之前到达公司,否则这个月工资清零。可是小A偏偏又有赖床的坏毛病。于是为了保住自己的工资,小A买了一个十分牛B的空间跑路器,每秒钟可以跑2^k千米(k是任意自然数)。当然,这个机器是用longint存的,所以总跑路长度不能超过maxlongint千米。小A的家到公司的路可以看做一个有向图,小A家为点1,公司为点n,每条边长度均为一千米。小A想每天能醒地尽量晚,所以让你帮他算算,他最少需要几秒才能到公司。数据保证1到n至少有一条路径。
Input
第一行两个整数n,m,表示点的个数和边的个数。
接下来m行每行两个数字u,v,表示一条u到v的边。
Output
- 一行一个数字,表示到公司的最少秒数。
Sample Input
4 4
1 1
1 2
2 3
3 4
Sample Output
1
Data Size
50%的数据满足最优解路径长度<=1000;
100%的数据满足n<=50,m<=10000,最优解路径长度<=maxlongint。
题解:
- 妙题,妙不可言。
- 首先2的k次方,想到什么?倍增!
- 然后n <= 50,想到什么?Floyd!
- 我们树上倍增时是怎么做的:dp(i, j) = dp(dp(i, j - 1), j - 1)
- 那么在这里就可以设dp(i, j, p)表示从i到j走2的p次方步能不能走到
- 转移就像Floyd那样枚举中转点k:
- if(dp(i, k, p - 1) == 1 && dp(k, j, p - 1) == 1) dp(i, j, p) = 1;
- 然后开个dis数组如果dp(i, j, p)=1,dis(i, j)=1,记录一下距离。
- 最后跑一遍最短路就行了。
#include <iostream>
#include <cstdio>
#include <cstring>
#define N 105
using namespace std;
int n, m;
int dis[N][N];
bool dp[N][N][N];
int main()
{
memset(dis, 0x3f, sizeof(dis));
cin >> n >> m;
for(int i = 1; i <= m; i++)
{
int u, v; cin >> u >> v;
dp[u][v][0] = 1;
dis[u][v] = 1;
}
for(int p = 1; p <= 64; p++)
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
for(int k = 1; k <= n; k++)
if(dp[i][j][p - 1] && dp[j][k][p - 1])
dp[i][k][p] = 1, dis[i][k] = 1;
for(int k = 1; k <= n; k++)
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
if(i != j && j != k && i != k)
dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
cout << dis[1][n];
return 0;
}