hdu6214--Smallest Minimum Cut(最小割边)

题目地址

I’m

题目大意

就是让你求出最小割应该割几条边

解题思路

利用dinic算法来计算。首先解题过程是给每个边的权值value变为value*(e+1)+1 (e为图的边数,乘的数字只要比边数大就行),然后跑一遍最大流,求最小割的边数就flow%(e+1)。
解释一下,先说加一是为啥吧!如果一个边是最小割的边,它一定是满流的,所以你容量加了1,它提供的流量肯定也加一,而其他的边肯定是不变的。所以最后原本的flow其实是变为了flow*(e+1)+p,p就是每个最小割边提供的那个1。那么对这个最大流模(e+1)就剩下了p,也就是最小割的边。
在说一下为什么乘(e+1)呢?首先看如果你对于每个边都直接加1,那最后的flow就变为了flow+p,这个p就没有办法求了,最大流也废了。那可能有人就想到了能不能跑两边,我第一遍先求原本的flow,然后给各个边加1再求flow+p然后相减。这个对错我也不太清楚,想尝试的可以试一试。我们有更简便的方法就是对边值乘上一个大于边数的数。这样的话,求出来的是flow*(e+1)+p,我们乘的数字是大于边数的数字,如果所有边都是最小割的边,那么p就等于m, 式子变为了flow*(e+1)+m如果e+1小于边数的话,在取模(e+1)就会影响到m的值。因此必须要乘以一个大于边数的数字才不会影响最小割的值。

AC代码

#include <iostream>
#include <queue>
#include <vector> 
#include <cstring>

using namespace std;

const int INF = 0x3f3f3f3f;
const int MAX_M = 1055;
const int MAX_V = 300;

int level[MAX_V];

struct edge
{
	int to, cap, rev;
};

vector<edge> G[MAX_V];

void add_edge(int from, int to, int cost)
{
	G[from].push_back((edge){to, cost, G[to].size()});
	G[to].push_back((edge){from, 0, G[from].size()-1});
}

void bfs(int s)
{
	memset(level, -1, sizeof(level));
	queue<int> que;
	level[s] = 0;
	que.push(s);
	
	while (!que.empty())
	{
		int v = que.front();
		que.pop();
		
		for (int i=0; i<G[v].size(); i++)
		{
			edge &e = G[v][i];
			
			if (e.cap && level[e.to] < 0)
			{
				level[e.to] = level[v]+1;
				que.push(e.to);
			}
		}
	}
}

int dfs(int v, int t, int f)
{
	if (v == t)
	{
		return f;
	}
	
	for (int i=0; i<G[v].size(); i++)
	{
		edge &e = G[v][i];
		
		if (e.cap && level[v]+1==level[e.to])
		{
			int d = dfs(e.to, t, min(f, e.cap));
			
			if (d > 0)
			{
				e.cap-=d;
				G[e.to][e.rev].cap+=d;
				return d;
			}
		}
	}
	
	return 0;
}

int dinic(int s, int t)
{
	int flow = 0;
	
	while (1)
	{
		bfs(s);
		
		if (level[t]<0)
		{
			return flow;
		}
		
		int f;
		
		while ((f = dfs(s, t, INF)) > 0)
			flow+=f;
	}
}

int main()
{
	int t;
	
	scanf("%d", &t);
	
	while (t--)
	{
		int n, m;
		int s, t;
		
		scanf("%d %d", &n, &m);
		scanf("%d %d", &s, &t);
		
		for (int i=0; i<=n; i++)
		{
			G[i].clear();
		}
		
		for (int i=0; i<m; i++)
		{
			int x, y, cost;
			
			scanf("%d %d %d", &x, &y, &cost);
			add_edge(x, y, cost*MAX_M+1);
		}
		
		printf("%d\n", dinic(s, t)%MAX_M);
	}
	
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值