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。

题解

对于这么一个非常稠密的有向图，他的边权全为1。我们可以这么考虑：先找出可以在一步内到达的点，然后其他的点之间无论怎么走都不能约掉二进制上的位（证明可以较为简单的脑补==）。
如果i->m的距离为2^k，m->j的距离为2^k，那么i->j的距离就会成为2^k+1，也就是说，i->j也可以一步到达了；
这样我们构造出新图，如果可以一步到达的两点，我们将他们之间连一条长为1的边，其余的点将边权设为正无穷（本体中64就够了==）。得到新图后，对新图进行最短路算法，就可以得到答案了；

代码

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#define maxn  100
using namespace std;

long long map[maxn][maxn];
int f[maxn][maxn][70];
int n, m;

int main() {
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= n; j++)
            map[i][j] = 100;

    for(int i = 1; i <= m; i++) {
        int u, v;
        scanf("%d%d", &u, &v);
        f[u][v][0] = 1;
        map[u][v] = 1;
    }
    for(int k = 1; k <= 63; k++)
        for(int i = 1; i <= n; i++)
            for(int w = 1; w <= n; w++) 
                if(f[i][w][k-1])
                    for(int j = 1; j <= n; j++) 
                        if(f[w][j][k-1]) {
                            f[i][j][k] = 1;
                            map[i][j] = 1;
                        }

    for(int k = 1; k <= n; k++)
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= n; j++)
               map[i][j] = min(map[i][j], map[i][k] + map[k][j]);
     cout << map[1][n];
}