POJ3281最大流

这里建图十分简单,不够很常用,把f和d分开,分成两个板块,分别连接可以连接的牛。
然后建立两边的连接,每一头牛分别对应自己,表示自己喜欢的f和d,
然后建立超级源点和超级汇点,跑最大流量。看最大的流量是多少就是多少了。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<climits>
#include<stack>
#include<vector>
#include<queue>
#include<set>
#include<map>
//#include<regex>
#include<cstdio>
#define up(i,a,b)  for(int i=a;i<b;i++)
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const long long INF = 0x3f3f3f3f;
const int inf = 1e9;
using namespace std;
typedef pair<int, int> pir;
int n, f, d;
struct edge { int to, cap, res; };
vector<edge>g[1000];
int level[1000];
int iter[1000];
void addedge(int s, int t)
{
	g[s].push_back(edge{ t,1,(int )g[t].size() });
	g[t].push_back(edge{ s,0,(int )g[s].size() - 1 });
}//邻接表
void bfs(int  s)//bfs算层数,然后后面会用这里求得层数,进行最短路径的寻找
{
	memset(level, -1, sizeof(level));
	queue<int >que;
	level[s] = 0;
	que.push(s);
	while (!que.empty())
	{
		int v = que.front();
		que.pop();
		up(i, 0, g[v].size())
		{
			edge &e = g[v][i];
			if (e.cap > 0 && 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 = iter[v]; i < g[v].size(); i++)//这里是算已经用过的,当前节点的
	//的边,因为这里找到一条增广路就退出了,但其实可能不止一条增广路,还会重新dfs,那么就要去除已经走过的边了 
	{
		edge &e = g[v][i];
		if (e.cap > 0 && level[v] < level[e.to])
		{
			int d = dfs(e.to, t, min(f, e.cap));
			if (d > 0)
			{
				e.cap -= d;
				g[e.to][e.res].cap += d;
				return d;
			}
		}
	}
	return 0;
}
int max_flow(int s, int t) {
	int flow = 0;
	while (1)
	{
		bfs(s);
		if (level[t] < 0)return flow;
		memset(iter, 0, sizeof(iter));
		int f;
		while ((f = dfs(s, t,inf)) > 0)//多次增广,直到不能为止
		{
			flow += f;
		}
	}
}//计算最大流
int main()
{
	cin >> n >> f >> d;
	int x, y,temp;
	up(i, 0, n)
	{
		cin >> x >> y;
		up(j, 0, x)
		{
			cin >> temp;
			temp--;
			addedge(temp, f + i);
		}
		up(j, 0, y)
		{
			cin >> temp;
			temp--;
			addedge(f + n + i, f+2*n+temp);
		}
		addedge(f + i, f + n + i);
	}
	up(i, 0, f)
	{
		addedge(f + 2 * n + d, i);
	}
	up(i, 0, d)
	{
		addedge( f + 2*n + i, f + 2 * n + d + 1);
	}//以上全是建图加边的操作
	int res = max_flow(f + 2 * n + d, f + 2 * n + d + 1);
	cout << res << endl;
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值