题意
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;
}