洛谷1613 跑路【倍增+floyd】

题目描述：

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

输入输出格式

解题思路：

由于距离是2的幂的距离的亮点可以1秒到达，所以我们可以用g[x][y][k]表示x到y的路径中是否用距离为$2^k$的，倍增+floyd处理即可，如果有，就从x到y连一条长度为1的边，最后跑一边最短路即可。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=55,INF=0x3f3f3f3f;
int n,m,f[N][N];
bool g[N][N][65];

int main()
{
    //freopen("lx.in","r",stdin);
    int x,y;
    memset(f,INF,sizeof(f));
    n=getint(),m=getint();
    while(m--)
    {
        x=getint(),y=getint();
        f[x][y]=g[x][y][0]=1;
    }
    for(int p=1;p<=64;p++)
        for(int k=1;k<=n;k++)
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++)
                {
                    g[i][j][p]=g[i][j][p]||g[i][k][p-1]&&g[k][j][p-1];
                    if(g[i][j][p])f[i][j]=1;
                }
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    cout<<f[1][n];
    return 0;
}