P1345 最小割

由题可知,已知无限图及起止点,求最小割点。
回忆紫书上的知识:最大流最小割,这里的割指的是割边
对于这个题,我们把每一个点拆成左右相连的两部分就等同于把最小割点问题转化成了最小割边问题。
证明:
如A图,假设一个点有两个乃至更多个相邻的割边,那么这个图的割边可以简化为图B所示的形式。以此推理,最后所有的割边都可以转换为左右点间的边。
图A
图B
至于无向图转有向图的技巧可以看这张图片

#include <queue>
#include <cstring>
#include <vector>
#include <iostream>
using namespace std;
const int maxn = 10000 + 10;
const int inf = 0x3f3f3f3f;
struct Edge
{
	int from, to, cap, flow;
};
struct Dinic
{
	int m = 0, s = 0, t = 0;
	vector<Edge> edges;
	vector<int> g[maxn];
	bool vis[maxn];
	int d[maxn];
	int cur[maxn];
	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 maxflow = 0, f;
		for (int &i = cur[x]; i < g[x].size(); i++)
		{
			Edge &e = edges[g[x][i]];
			if (d[x] + 1 == d[e.to] && (f = dfs(e.to, min(a, e.cap - e.flow))) > 0)
			{
				e.flow += f;
				edges[g[x][i] ^ 1].flow -= f;
				maxflow += f;
				a -= f;
				if (a == 0)
					break;
			}
		}
		return maxflow;
	}
	int Maxflow(int start, int end)
	{
		s = start;
		t = end;
		int flow = 0;
		while (bfs())
		{
			memset(cur, 0, sizeof(cur));
			flow += dfs(s, inf);
		}
		return flow;
	}
	void clear(int n)
	{
		edges.clear();
		for (int i = 0; i <=n; i++)
			g[i].clear();
	}
	void addedge(int from, int to, int cap)
	{
		edges.push_back((Edge){from, to, cap, 0});
		edges.push_back((Edge){to, from, 0, 0});
		m = edges.size();
		g[from].push_back(m - 2);
		g[to].push_back(m - 1);
	}
};
Dinic dn;
int main()
{
	int n,m,s,t;
	int u,v;
	scanf("%d%d%d%d",&n,&m,&s,&t);
	s--;
	t--;
	s+=n;
	dn.clear(2*n+2);
	for(int i=0;i<n;i++){
		if(i==s||i==t) continue;
		dn.addedge(i,i+n,1);
	}
	for(int i=0;i<m;i++){
		scanf("%d%d",&u,&v);
		u--;
		v--;
		dn.addedge(u+n,v,1);
		dn.addedge(v+n,u,1);
	}
	int ans=dn.Maxflow(s,t);
	printf("%d",ans);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值