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