网络流一:最大流EK算法(bfs)

//POJ 1273  //https://blog.csdn.net/A_Comme_Amour/article/details/79356220 
//网络流最大流EK算法 
// 时间复杂度上限为nm^2 
#include <iostream>
#include <cstdio>
#include <queue>
#include <cstring>

using namespace std;

#define maxn 0x7fffffff 

int n, m;
int gra[210][210] = {0};
int pre[210] = {0};
int flow[210] = {0}; //记录源点到这个点的流量 
queue<int>qu;
int flowMax = 0; // 累加最大流 

int read()
{
	int x = 0, w = 1;
	char ch = getchar();
	while(ch <='0' || ch > '9')
	{
		if(ch == '-') w = -1;
		ch = getchar();
	}
	while(ch <= '9' && ch >= '0')
	{
		x = x * 10 + ch - '0';
		ch = getchar();
	}
	return x * w;
}

int bfs(int s, int t) // 寻找增广路 
{
	while(!qu.empty())
	{
		qu.pop();
	} // 代价小 直接循环清 
	for(int i = 1; i <= n; i++)
	{
		pre[i] = -1;
	}  
	pre[s] = 0; // 考虑存在一条路重新回到s 
	flow[s] = maxn; // 便于更新增广路上的flow 
	qu.push(s);
	while(!qu.empty())
	{
		int now = qu.front();
		qu.pop();
		if(now == t) break;
		for(int i = 1; i <= n; i++)
		{
			if(gra[now][i] > 0 && pre[i] == -1)
			{
				pre[i] = now;
				flow[i] = min(flow[now], gra[now][i]);
				qu.push(i);
			}
		}
	}
	if(pre[t] == -1)
	{
		return -1;
	}
	return flow[t];
}

void EK(int s, int t)// 通过pre更新
{
	flowMax = 0; 
	memset(flow, 0, sizeof(flow));
	int increase = 0;// 本次增广路上的最小流flow 
	while(increase = bfs(1, n))
	{
		if(increase == -1)
		{
			break; // 增广路全部寻找完 
		}
		int now = t; 
		while(now != s)
		{
			gra[pre[now]][now] -= increase;
			gra[now][pre[now]] += increase; // 建立反向边 保证最优 
			now = pre[now];
		}
		flowMax += increase;
	}
}

int main()
{
	while(scanf("%d %d", &m, &n) != EOF)
	{
		memset(gra, 0, sizeof(gra));
		for(int i = 1; i <= m; i++)
		{
			int a = read();
			int b = read();
			int c = read();
			gra[a][b] += c; //题目特定 
		}
	
		EK(1, n);  
		
		printf("%d\n", flowMax);     
	}
	
	
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值