在某一个国家,那儿有n个城市,他们通过m条双向道路相连。城市从1到n编号。如果城市a和b通过一条道路直接相连,那么他们之间的距离就是一个小时。这个国家的道路网络可以允许你从任意一个城市到达另外的城市。
现在你要破坏尽可能多的道路,但是要保证从城市s1到t1不超过l1小时,并且从城市s2到t2不超过l2小时。
输出最多可以破坏的道路数目,如果没有解,请输出-1
Input
单组测试数据。
第一行有两个整数n,m(1 ≤ n ≤ 3000, n-1 ≤ m ≤ min(3000,n*(n-1)/2) )。
接下来m行,每行有两个整数 ai, bi (1 ≤ ai, bi ≤ n, ai ≠ bi),表示ai和bi之间有一条道路。
输入保证是一个连通图。
最后两行每行有三个整数s1, t1, l1和 s2, t2, l2, (1 ≤ si, ti ≤ n, 0 ≤ li ≤ n)。
Output
输出一个整数,表示最多可以破坏的道路数目,如果没有解,输出-1。
Input示例
5 4
1 2
2 3
3 4
4 5
1 3 2
3 5 2
Output示例
0
思路:
先用bfs计算出任意两点间的最短距离。再考虑去掉s1-t1和s2-t2的路径重合部分。
#include <iostream>
#include <algorithm>
#include <string.h>
#include <queue>
#include <vector>
using namespace std;
const int MAXN = 3005;
int n, m;
int s1, t1, l1;
int s2, t2, l2;
vector<int> line[MAXN];
int dis[MAXN][MAXN];
bool vis[MAXN];
void bfs()
{
for (int i = 1; i <= n; i++)
{
memset(vis, 0, sizeof(vis));
queue<int> buf;
buf.push(i);
vis[i] = true;
while (!buf.empty())
{
int temp = buf.front();
buf.pop();
int size = line[temp].size();
for (int j = 0; j < size; j++)
{
int v = line[temp][j];
if (!vis[v])
{
vis[v] = true;
dis[i][v] = dis[i][temp] + 1;
buf.push(v);
}
}
}
}
}
int main ()
{
cin >> n >> m;
memset(dis, 0, sizeof(dis));
int a, b;
for (int i = 1; i <= m; i++)
{
cin >> a >> b;
line[a].push_back(b);
line[b].push_back(a);
}
cin >> s1 >> t1 >> l1;
cin >> s2 >> t2 >> l2;
bfs();
int result = dis[s1][t1] + dis[s2][t2];
if (dis[s1][t1] > l1 || dis[s2][t2] > l2)
{
cout << "-1" << endl;
}
else
{
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= n; j++)
{
if (dis[s1][i] + dis[i][j] + dis[j][t1] <= l1 && dis[s2][i] + dis[i][j] + dis[j][t2] <= l2)
{
result = min(result, dis[s1][i] + dis[i][j] + dis[j][t1] + dis[s2][i] + dis[j][t2]);
}
if (dis[t1][i] + dis[i][j] + dis[j][s1] <= l1 && dis[t2][i] + dis[i][j] + dis[j][s2] <= l2)
{
result = min(result, dis[t1][i] + dis[i][j] + dis[j][s1] + dis[t2][i] + dis[j][s2]);
}
if (dis[s1][i] + dis[i][j] + dis[j][t1] <= l1 && dis[t2][i] + dis[i][j] + dis[j][s2] <= l2)
{
result = min(result, dis[s1][i] + dis[i][j] + dis[j][t1] + dis[t2][i] + dis[j][s2]);
}
if (dis[t1][i] + dis[i][j] + dis[j][s1] <= l1 && dis[s2][i] + dis[i][j] + dis[j][t2] <= l2)
{
result = min(result, dis[t1][i] + dis[i][j] + dis[j][s1] + dis[s2][i] + dis[j][t2]);
}
}
}
cout << m - result << endl;
}
return 0;
}