洛谷P1345 奶牛的電信Telecowmunication 最小割

給出n1500n\leq1500個點,然後有兩個點c1,c2c_{1},c_{2}會通過某些通路相連,然後你可以去除掉某些點,如果某個點被去掉。那麼和這個點相連接的邊也會被去除。然後求問去掉最少的點,使得兩個點不連通。
拆點,中間連接一條容量爲11的邊,然後點與點的連接從下面的點連接到上面的點。求最小割。
對於那些認爲不能拆開的點,一般連infinf的流量。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;

namespace Dinic {
	const int N=1e4+3,M=1e5+3; 
	int front[N],head[N],dep[N],tot=-1,n,S,T;
	struct edge{ int v,nxt;ll w; }e[2*M];
	void init() {
		tot=-1;
		for(int i=1;i<=n;i++) front[i]=-1;
	}
	void add(int u,int v,ll w) {
		e[++tot].v=v,e[tot].w=w,e[tot].nxt=front[u],front[u]=tot;
		e[++tot].v=u,e[tot].w=0,e[tot].nxt=front[v],front[v]=tot;
	}
	bool bfs() {
		for(int i=1;i<=n;i++) dep[i]=0;
		queue<int> q;
		q.push(S);
		dep[S]=1;
		while(!q.empty()) {
			int u=q.front();
			q.pop();
			for(int i=front[u];i!=-1;i=e[i].nxt) {
				int v=e[i].v;
				if(dep[v]==0&&e[i].w) {
					dep[v]=dep[u]+1;
					q.push(v); 
				}
			}
		}
		return dep[T];
	}
	ll dfs(int cur,ll dis=1e18) {
		if(cur==T) return dis;
		ll flow=0,sum=0;
		for(int &i=head[cur];i!=-1;i=e[i].nxt) {
			int v=e[i].v;
			if(dep[v]==dep[cur]+1&&e[i].w) {
				flow=dfs(v,min(dis,e[i].w));
				dis-=flow,sum+=flow;
				e[i].w-=flow,e[i^1].w+=flow;
				if(!dis) break;
			}
		}
		return sum;
	}
	ll solve() {
		ll ans=0;
		while(bfs()) {
			for(int i=1;i<=n;i++) head[i]=front[i];
			ans+=dfs(S);
		}
		return ans;
	}
}

int main() {
	int n,m,c1,c2;
	scanf("%d%d%d%d",&n,&m,&c1,&c2);
	Dinic::n=2*n;
	int S=Dinic::S=c1;
	int T=Dinic::T=c2+n;
	Dinic::init();
	for(int i=1;i<=n;i++) {
		if(i!=c1&&i!=c2) Dinic::add(i,i+n,1);
		else Dinic::add(i,i+n,1e9);
	}
	for(int i=1;i<=m;i++) {
		int u,v;
		scanf("%d%d",&u,&v);
		Dinic::add(u+n,v,1e9);
		Dinic::add(v+n,u,1e9);
	}
	int ans=Dinic::solve();
	printf("%d\n",ans);  
	return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章