题目描述
小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的边。
输出格式:
一行一个数字,表示到公司的最少秒数。
【数据范围】
50%的数据满足最优解路径长度<=1000;
100%的数据满足n<=50,m<=10000,最优解路径长度<=maxlongint。
题目分析
题目中2^k很显然就是要和倍增扯上关系
dis[ ][ ]
d
i
s
[
]
[
]
数组一开始先保存原图
map[u,v,k]
m
a
p
[
u
,
v
,
k
]
表示u到v是否存在一条总长度为
2k
2
k
的路径
可以用类似Floyd的做法预处理出map数组
若有
map[u,v,k]==1
m
a
p
[
u
,
v
,
k
]
==
1
,则令
dis[u,v]=1
d
i
s
[
u
,
v
]
=
1
最后用dis数组跑Floyd
最后答案就是
dis[1,n]
d
i
s
[
1
,
n
]
#include<iostream>
#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
int read()
{
int f=1,x=0;
char ss=getchar();
while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
return f*x;
}
int n,m;
int dis[100][100],map[100][100][100];
int main()
{
n=read();m=read();
memset(dis,33,sizeof(dis));
for(int i=1;i<=n;++i) dis[i][i]=0;
for(int i=1;i<=m;++i)
{
int u=read(),v=read();
dis[u][v]=map[u][v][0]=1;
}
for(int k=1;k<=64;++k)
for(int t=1;t<=n;++t)
for(int i=1;i<=n;++i)
for(int j=1;j<=n;++j)
{
if(map[i][t][k-1]&map[t][j][k-1])
map[i][j][k]=dis[i][j]=1;
}
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]);
cout<<dis[1][n];
return 0;
}