Codeforces 238E：Meeting Her

【大意】
有一个n 个结点的有向图，边权均为1。Urapl 想从a 出发去b。有p 个公交车公司。在每
一秒的开始，第i 个公司的公交车随机选择一条从si 到ti 的最短路径然后走这条路径。如果
一个公交车经过Urpal 所在的交叉点，则Urpal 可以上这辆公交车，他可以在中途任意一个结
点下车。
在任何时刻Urpal 只知道他自己的位置和约会地点。当他上了公交车时他只知道这辆公交
车属于第几个公司。当然Urpal 知道城市地图和每个公司的(si; ti)。
求最坏情况下他需要乘坐公交车的次数。
【解答】
先Floyd跑出任意两点的最短路，用这个来求每个巴士的必经点，显然只有在该巴士的必经点才能选择上这个巴士。
每个巴士在每个点会随机走一个方向，为了保证答案最坏所以应该走的是花费最大的那一条路。
每一个状态可以用（位置，巴士）这样的二元组来表示，共有$O(N^2)$种状态，但是无法转移。用类似最短路的方法求答案，时间复杂度用Bellman_Ford类似的方式计算是$O(N^4)$

#include <bits/stdc++.h>
#define N 105
using namespace std;

int n,m,p,S,T;
int dis[N][N],tmp[N],vis[N],tme;
int ans[N],s[N],t[N],c[N][N],mark[N];

inline int rd() {
    int x=0,f=1;char ch=getchar();
    while (ch>'9'||ch<'0') {if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int dfs(int x,int y) {
    if (vis[x] == tme) return tmp[x];
    vis[x] = 1;
    int cur = -1;
    for (int i=1;i<=n;i++)
        if (dis[x][i] == 1 && 1+dis[i][y] == dis[x][y])
            cur = max(cur, dfs(i,y));
    if (cur==-1) cur = 1<<29;
    cur = min(cur, ans[x]+1);
    return tmp[x] = cur;
}

int main() {
    n = rd(), m = rd(), S = rd(), T = rd();
    memset(dis,127/3,sizeof(dis));
    for (int i=1;i<=n;i++) dis[i][i] = 0;
    for (int _=1;_<=m;_++) {
        int a = rd(), b = rd();
        dis[a][b] = 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]);

    p = rd();
    for (int i=1;i<=p;i++)
        s[i] = rd(), t[i] = rd();

    for (int i=1;i<=p;i++) {
        memset(mark,0,sizeof(mark));
        int x = s[i], y = t[i];
        for (int j=1;j<=n;j++)
            if (dis[x][j] <= n && dis[x][y] == dis[x][j]+dis[j][y])
                mark[ dis[x][j] ]++;
        for (int j=1;j<=n;j++)
            if (dis[x][j] <= n && dis[x][y] == dis[x][j]+dis[j][y]) {
                if (mark[ dis[x][j] ] == 1) c[i][j] = 1;
            }
    }

    memset(ans,127/3,sizeof(ans));
    ans[T] = 0;
    for (;;) {
        int flag = 0;
        for (int i=1;i<=p;i++)
            for (int j=1;j<=n;j++) if (c[i][j]) {
                tme++;
                int x = dfs(j,t[i]);
                x<ans[j] ? flag = 1, ans[j]=x :0;
            }
        if (!flag) break;
    }
    ans[S]>n ? puts("-1") : printf("%d\n",ans[S]);
    return 0;
}