[USACO5.4]奶牛的电信Telecowmunication题解

[USACO5.4]奶牛的电信Telecowmunication题解

题目描述

农夫约翰的奶牛们喜欢通过电邮保持联系，于是她们建立了一个奶牛电脑网络，以便互相交流。这些机器用如下的方式发送电邮：如果存在一个由 $c$ 台电脑组成的序列$a_1,a_2,\cdots ,a_c$，且 $a_1$ 与 $a_2$ 相连，$a_2$ 与 $a_3$ 相连，等等。那么电脑 $a_1$ 和 $a_c$ 就可以互发电邮。

很不幸，有时候奶牛会不小心踩到电脑上，农夫约翰的车也可能碾过电脑，这台倒霉的电脑就会坏掉。这意味着这台电脑不能再发送电邮了，于是与这台电脑相关的连接也就不可用了。

有两头奶牛就想：如果我们两个不能互发电邮，至少需要坏掉多少台电脑呢？请注意，$c_1,c_2$ 不能被破坏。请编写一个程序为她们计算这个最小值。

以如下网络为例：

   1*
  /
 3 - 2*

这张图画的是有 $2$ 条连接的 $3$ 台电脑。我们想要在电脑 $1$ 和 $2$ 之间传送信息。电脑 $1$ 与 $3$，$2$ 与 $3$ 直接连通。如果电脑 $3$ 坏了，电脑 $1$ 与 $2$ 便不能互发信息了。

输入格式

第一行：四个由空格分隔的整数：$N,M,c_1,c_2$。$N$ 是电脑总数，电脑由 $1$ 到 $N$ 编号。$M$ 是电脑之间连接的总数。后面的两个整数 $c_1$ 和 $c_2$ 是上述两头奶牛使用的电脑编号。连接没有重复且均为双向的（即如果 $c_1$ 与 $c_2$ 相连，那么 $c_2$ 与 $c_1$ 也相连）。两台电脑之间至多有一条连接。电脑 $c_1$ 和 $c_2$ 不会直接相连。

第 $2$ 到 $M+1$ 行：接下来的 $M$ 行中，每行包含两台直接相连的电脑的编号。

输出格式

一行，一个整数，表示使电脑 $c_1$ 和 $c_2$ 不能互相通信需要坏掉的电脑数目的最小值。

样例 #1

样例输入 #1

3 2 1 2
1 3
2 3

样例输出 #1

1

提示

对于 $100\%$ 的数据：$1\le N\le 100$，$1\le M\le 600$。

Solution

​ 观察题目，发现要求的是割去几个点能实现$c_1$与$c_2$不联通，但我们的最小割是割去边，因此考虑将割点转化为割边，考虑如何转化。

​ 对于每个点$x$，我们试图让它与一条能被割去的边对应，自然想到拆点，将每个点$x$拆为$2x-1$，$2x$两个点，两点间用一条代价为$1$的边连接。点$2x-1$作为$x$的入点，点$2x$作为$x$的出点，入点与原本连向$x$的点的出点连代价无穷大的边，出点与x原本向外连的点的入点连代价无穷大的边。这样一来，我们可以保证所有经过$x$的路径都必须经过$2x-1$与$2x$间代价为$1$的边。

​ 建图后跑最小割（最大流）即可

小贴士：记得网络流的起点和终点设成$c_1$的出点和$c_2$的入点

示例代码

//#pragma GCC optimize(3)
//#pragma GCC optimize("inline")
//#pragma GCC optimize("Ofast")
//#pragma GCC target("sse3","sse2","sse")
//#pragma GCC diagnostic error "-std=c++14"
//#pragma GCC diagnostic error "-fwhole-program"
//#pragma GCC diagnostic error "-fcse-skip-blocks"
//#pragma GCC diagnostic error "-funsafe-loop-optimizations"
//#pragma GCC optimize("fast-math","unroll-loops","no-stack-protector","inline")
#include<bits/stdc++.h>
#define ll long long
#define N 100005
#define inf 1000000007
#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fd(i,a,b) for(int i=a;i>=b;--i)
using namespace std;
int n,m,s,t,tot=1;
int cur[N],d[N];
struct graph{
	int st,to,nxt,cap;
}g[N];
int read(){
	int res=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-'&&((ch=getchar())>='0'&&(ch)<='9')) f=-1;
		else ch=getchar();
	}
	while(ch>='0'&&ch<='9')res=res*10+ch-'0',ch=getchar();
	return res*f;
}
void add(int u,int v,int w){
	g[++tot].nxt=g[u].st,g[tot].to=v,g[tot].cap=w,g[u].st=tot;
}
void Add(int u,int v,int w){
	add(u,v,w),add(v,u,0);
}
bool bfs(){
	queue<int> q;
	memset(d,-1,sizeof d);
	q.push(s);
	d[s]=0;
	while(!q.empty()){
		int u=q.front();
		q.pop();
		for(int i=g[u].st;i;i=g[i].nxt){
			int v=g[i].to;
			if(d[v]!=-1||g[i].cap==0)continue;
			d[v]=d[u]+1;
			q.push(v);
		}
	}
	return d[t]!=-1;
}
int dfs(int x,int y){
	if(x==t||y==0)return y;
	int flow=y;
	for(int &i=cur[x];i;i=g[i].nxt){
		int v=g[i].to;
		if(d[v]==d[x]+1&&g[i].cap){
			int z=dfs(v,min(flow,g[i].cap));
			g[i].cap-=z,g[i^1].cap+=z;
			flow-=z;
			if(flow==0)break;
		}
	}
	if(flow==y)d[x]=-1;
	return y-flow;
}
int maxflow(){
	int flow=0;
	while(bfs()){
		fo(i,1,2*n){
			cur[i]=g[i].st;
		}
		flow+=dfs(s,inf);
	}
	return flow;
}
int main(){
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	n=read(),m=read(),s=read(),t=read();
	s=2*s,t=2*t-1;
	fo(i,1,n)Add(2*i-1,2*i,1);
	fo(i,1,m){
		int u=read(),v=read();
		Add(2*u,2*v-1,inf),Add(2*v,2*u-1,inf);
	}
	printf("%d",maxflow());
}