USACO Training Section 5.4-Telecowmunication 奶牛的电信【最小割边集字典序最小】woj1791

题意

n个点m条边的无向图，问最少破坏多少个点可以使得S,T不连通，按字典序最小输出方案

解法

对于第一个问，显然的一个拆点技巧，将一个点i拆分为2个点i,i’,i->i’,容量为1
原图的边<u,v>，容量设为inf，因为每个点拆分为2个点，所以建边<u’ ,v ,inf>,<v’ , u ,inf>
然后求出最大流即最小割
对于第二问，因为最大流后求出的最小割是任意最小割。所以在这里，我们可以采取依次枚举每个点，将它拆开后的边容量变为0，然后跑最大流，看最大流是否改变，如果改变了则输出。

参考代码

#include<iostream>
#include<vector>
#include<queue>
#include<cstdio>
#include<cstring>
using namespace std;
#define ll long long
#define in rad()
inline int rad(){
	int x=0,f=1;char c=getchar();while(c>'9'||c<'0')c=getchar();
	if(c=='-')f=-1,c=getchar();while(c>='0'&&c<='9')x=x*10+c-48,c=getchar();
	return x*f;
}
const int maxn=1e3+10;
const int maxm=4e4+100;
const int inf=1e9;
struct Edge {
	int u,v,w,nxt;
} e[maxm],ee[maxm];
int first[maxn]= {},cnt=1,S,T;
inline void add(int u,int v,int w) {
	e[++cnt].v=v;e[cnt].u=u;e[cnt].w=w;e[cnt].nxt=first[u];first[u]=cnt;
}
inline void ins(int u,int v,int w){
	add(u,v,w);add(v,u,0);
}
int n,m;
namespace D{
	int dis[maxn],cur[maxn];
	int bfs(int s,int t){
		memset(dis,-1,sizeof(dis));
		queue<int>q;
		dis[s]=0;q.push(s);
		while(!q.empty()){
			int u=q.front();q.pop();
			for(int i=first[u];i;i=e[i].nxt){
				int v=e[i].v;
				if(dis[v]!=-1||e[i].w==0)continue;
				dis[v]=dis[u]+1;
				q.push(v);
				if(v==t)return 1;
			}
		}
		return 0;
	} 
	int dfs(int u,int t,int f){
		if(!f||u==t)return f;
		int w,used=0;
		for(int &i=cur[u];i;i=e[i].nxt){
			int v=e[i].v;
			if(dis[v]!=dis[u]+1||e[i].w==0)continue;
			w=dfs(v,t,min(f,e[i].w));
			used+=w;f-=w;
			e[i].w-=w;e[i^1].w+=w;
			if(f==0)break;
		}
		return used;
	}
	int dinic(int s,int t){
		int ret=0;
		while(bfs(s,t)){
			memcpy(cur,first,sizeof(first));
			ret+=dfs(s,t,inf);
		}
		return ret;
	}
}
 
int main(){
 	n=in;m=in;S=in;T=in;
 	for(int i=1;i<=n;i++)
 		ins(i,i+n,1);//拆点
	for(int i=1;i<=m;i++){
		int u=in,v=in ;
		ins(u+n,v,inf);ins(v+n,u,inf);
	} 
	memcpy(ee,e,sizeof(e));
	S+=n;
	int ans=D::dinic(S,T); 
	printf("%d\n",ans);
	for(int i=2;i<=2*n;i+=2){		 
		ee[i].w=0;
		memcpy(e,ee,sizeof(e));
		int tmp=D::dinic(S,T);
		if(tmp==ans) //如果最大流没有改变，则不是割边 
			ee[i].w=1;
		else
			printf("%d ",i/2),ans=tmp;//最大流改变，这条边就是在割集中 
	} 
	return 0;
}