YbtOJ 跑路上班 &luogu1613 (基于floyd倍增表+最短路)

题目

仔细想一想，对于从 $i$ 走 $2^k(k\in N且k\le 64)$ 到 $j$ ，对答案的贡献为1，等效于此时$i$ 与 $j$ 间有一条边。那我们就可以用基于floyd思想的倍增表去维护这种边。

具体的，设$m{p}_{k,i,j,ind}$ 表示只有前 $k$ 个点时, $i$ 到 $j$ 之间走$2^{ind}$ 是否能到达。类比floyd,有转移（floyd）

$m{p}_{k,i,j,ind}=m{p}_{k-1,i,j,ind}\Vert (m{p}_{k-1,i,k,ind-1}\&m{p}_{k-1,k,j,ind-1})$

类比floyd，其中 $k$ 这一维可以滚动掉,而 $ind$ 作为阶段放在最外层，从而维护出一个倍增表，可以维护可达性。

之后随便一种最短路都行，由于这里floyd写起来方便，而且数据很小，就直接floyd了。

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxn=55;
bool mp[maxn][maxn][65];
int n,m,dis[maxn][maxn];
signed main()
{
    memset(dis,0x3f,sizeof(dis));
    cin>>n>>m;
    int t=64;
    for(int i=1,x,y;i<=m;i++)
        cin>>x>>y,mp[x][y][0]=1,dis[x][y]=1;
    for(int ind=1;ind<=t;ind++)
        for(int k=1;k<=n;k++)
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++)
                    mp[i][j][ind]|=mp[i][k][ind-1]&mp[k][j][ind-1];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            for(int k=1;k<=64;k++)
                if(mp[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]<<endl;
    return 0;
}//能过,但是大常数qwq