地址:点击打开链接
总结:
拆点,顾名思义,就是把一个点拆成两个点。因为在网络流的模型中,割特指边集而非点集,所以想要实现“可以被割掉的点”,就用用到拆点思想。另外,这种方法也可以被理解为:使点同边一样有一个容量上限。把点A分成两个,一个叫A1,表示A的“入点”;另一个叫A2,表示A的“出点”。A的所有出边都从A2连出,所有入边都连向A1。再连一条边从A1指向A2,其边权为点A的“点容量”。
思路:
本题求至少要删除多少个点使C1、C2不连通,问题的实质是一个最小割,跑一个最大流即可,因为只允许割点,不允许割边,所以原有的边边权都为INF。点权为1,表示每个点最多只能被删除一次。
代码:
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<queue>
#include<cstring>
using namespace std;
const int maxn=605,INF=0x3f3f3f3f;
//---------------------------Dinic开始
struct Dinic{
int n,s,t;
struct Edge{
int from,to,cap,flow;
};
vector<int>G[maxn];
vector<Edge>edges;
int cur[maxn],d[maxn];
bool vis[maxn];
void addedge(int f,int t,int c){ //加边
edges.push_back((Edge){f,t,c,0});
edges.push_back((Edge){t,f,0,0});
int m=edges.size();
G[f].push_back(m-2);
G[t].push_back(m-1);
}
bool BFS(){ //搜寻路径
memset(vis,0,sizeof(vis));
queue<int>Q;
Q.push(s);d[s]=0;vis[s]=1;
while(!Q.empty()){
int x=Q.front();Q.pop();
for(int i=0;i<G[x].size();i++){
Edge& e=edges[G[x][i]];
if(!vis[e.to] && e.cap>e.flow){
vis[e.to]=1;
d[e.to]=d[x]+1;
Q.push(e.to);
}
}
}
return vis[t];
}
int DFS(int x,int a){ //更新
if(x==t || a==0)return a;
int flow=0,f;
for(int& i=cur[x];i<G[x].size();i++){
Edge& e=edges[G[x][i]];
if(e.cap>e.flow && d[e.to]==d[x]+1 && (f=DFS(e.to,min(a,e.cap-e.flow)))>0){
e.flow+=f;
edges[G[x][i]^1].flow-=f;
flow+=f;
a-=f;
if(a==0)break;
}
}
return flow;
}
int MaxFlow(int S,int T){
s=S;t=T;int flow=0;
while(BFS()){
memset(cur,0,sizeof(cur));
flow+=DFS(S,INF);
}
return flow;
}
}dinic;
//--------------------------------------Dinic结束
int main(){
int N,M,c1,c2;scanf("%d%d%d%d",&N,&M,&c1,&c2);
for(int i=1;i<=N;i++){
dinic.addedge(i,i+N,1); //拆点,权值为1,表示只能被割一次
}
for(int i=1;i<=M;i++){
int x,y;scanf("%d%d",&x,&y);
dinic.addedge(x+N,y,INF);//x+N为出点
dinic.addedge(y+N,x,INF);
}
int maxx=dinic.MaxFlow(c1+N,c2); //注意是c1+n
printf("%d\n",maxx);
return 0;
}