nyoj677 网络流 最大流最小割 EK算法

碟战

时间限制： 2000 ms  |  内存限制： 65535 KB

难度： 4

描述

知己知彼，百战不殆！在战争中如果被敌人掌握了自己的机密，失败是必然的。K国在一场战争中屡屡失败，就想到自己的某些城市可能会有敌方的间谍。

在仔细调查后，终于得知在哪些城市存在间谍。当然这个消息也被敌方间谍得知，所以间谍们开始撤离，试图到达K国唯一机场，然后抢夺飞机回国。由于城市内部比较复杂，K国领导人决定封锁道路，阻止所有间谍到达机场。城市编号为1~N，两个城市有不超过1条双向道路相连。机场在N号城市，不会有间碟。

由于要节约兵力，至少要封锁多少条道路才能阻止所有间谍到达机场？

输入

第一行包含一个整数T（T <= 100)，为测试数据组数。
接下来每组测试数据第一行包含三个整数n,m,p（2<= n <= 200,1< m < 20000,1 <= p < n),分别表示城市数量，道路数量，存在间谍的城市的数量。
接下来的一行包含p个整数x（1 <= x < n),表示存在间谍城市的编号。
接下来的m行，每行包含两个整数i,j，表示城市i与城市j有道路相通。

输出

输出“Case #i: ans”(不含引号),i为第i组数据，ans为需要封锁道路的条数。

样例输入

2
4 4 2
1 2
1 2
2 4
1 3
3 4
4 3 2
1 2
2 3
3 4
2 4

样例输出

Case #1: 2
Case #2: 2

将0作为源点，与所有有间谍的城市相连，容量无穷大，然后城市之间的道路容量为1，汇点为N，注意是无向图，所以输入map[x][y]=map[y][x]=1;求最大流用最简单的EK算法，直接套模板，不懂的可以见上一篇文章EK算法，BFS求增广路；

#include<iostream>
#include<queue>
#include<cstring> 
#include<cstdio>
using namespace std;
const int maxn=205;
const int inf=0x7fffffff;

int map[maxn][maxn]; //残留网络，初始化为原图
bool visit[maxn];
 int pre[maxn];
 //int m,n;
int n;

bool bfs(int s,int t)  //寻找一条从s到t的增广路，若找到返回true
 {
     int p;
     queue<int > q;
     memset(pre,-1,sizeof(pre));
     memset(visit,false,sizeof(visit));

    pre[s]=s;
     visit[s]=true;
     q.push(s);
     while(!q.empty())
     {
         p=q.front();
         q.pop();
         for(int i=1;i<=n;i++)
         {
             if(map[p][i]>0&&!visit[i])
             {
                 pre[i]=p;
                 visit[i]=true;
                 if(i==t) return true;
                 q.push(i);
             }
         }
     }
     return false;
 }

int EdmondsKarp(int s,int t)
 {
    int flow=0,d,i;
    while(bfs(s,t))
    {
        //d=inf;  因为所有道路都是容量都是1，所以不用找最小值
        //for(i=t;i!=s;i=pre[i])
        //   d=d<r[pre[i]][i]? d:r[pre[i]][i];
        for(i=t;i!=s;i=pre[i])
        {
            map[pre[i]][i]-=1;
            //map[i][pre[i]]+=1;本来需要反向容量加一
            //map[pre[i]][i]-=1;但是这里反向边不需要。。。可能是双向边的关系
	}
        flow+=1;
    }
    return flow;
 }


 int main()
 {
 	int T,t,x,y,m,k;
 	int count=1;
	 scanf("%d",&T);
	 while(T--){
	 	scanf("%d%d%d",&n,&m,&k);
	 	memset(map,0,sizeof(map));
	 	for (int i=0;i<k;i++){
	 		scanf("%d",&t);
	 		map[0][t]=inf;
		 }
		 for (int i=0;i<m;i++){
		 	scanf("%d%d",&x,&y);
		 	map[x][y]=1;
		 	map[y][x]=1;
		 }
         printf("Case #%d: %d\n",count++,EdmondsKarp(0,n));
     }
     return 0;
 }