51nod - 1444 破坏道路

在某一个国家，那儿有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;
}