Hdu4292_Food(最大流)

题意:

有F种食物和D种饮料,每种食物或饮料数量有限,且每个人只享用一种食物和一种饮料。现在有n个顾客,每个人都有自己喜欢的食物和饮料,问最多能使几个人同时享用到自己喜欢的食物和饮料。

思路:

每个人只用一份food和drink就行了,所以把人拆点,中间连一条容量为1的边。
源点——>food——>人——>人——>drink——>汇点
源点——>food, 容量为食物的存量
food——>人,容量为INF
人——>人,容量为1
人——>drink,容量为INF
drink——>汇点,容量为drink的存量
跑最大流即可。

#include <cstdio>
#include <cstring>
#include <string>
#include <queue>
#include <algorithm>
const int INF = 0x3f3f3f3f;
const int maxn = 200+5;
using namespace std;

int food[maxn], drink[maxn];
char Gf[maxn][maxn], Gd[maxn][maxn];
// 图 
struct Edge{
	int u, v, cap, flow;
	Edge(int a, int b, int c, int d):u(a),v(b),cap(c),flow(d){}
};
vector<Edge> edges;
vector<int> G[maxn*5];

int dis[maxn*5];	// 分层的编号 
int cur[maxn*5]; // 当前弧,重要优化!!! 

void init(int a){
	for(int i = 0; i < a; ++i) G[i].clear();
	edges.clear();
}

void addEdge(int u, int v, int cap){
	edges.push_back(Edge(u,v,cap,0));
	edges.push_back(Edge(v,u,0,0));		// 反向弧
	int m = edges.size();
	G[u].push_back(m-2);
	G[v].push_back(m-1);
}

// 分层 
bool bfs(int s, int t){
	memset(dis, -1, sizeof(dis));
	dis[s] = 0;
	queue<int> Q;
	Q.push(s);
	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(dis[e.v] == -1&&e.cap > e.flow){
				dis[e.v] = dis[x] + 1;
				Q.push(e.v);
			}
		}
	}
	return dis[t] != -1;
}

int dfs(int s, int t, int cur_flow){
	if(s == t||cur_flow == 0) return cur_flow;
	int ans = 0;
	for(int& i = cur[s]; i < G[s].size(); ++i){
		int c = G[s][i];
		Edge e = edges[c];
		if(dis[e.v] == dis[s] + 1&&e.cap > e.flow){
			int a2 = min(cur_flow, e.cap-e.flow);
			int w = dfs(e.v, t, a2);
			edges[c].flow += w;
			edges[c^1].flow -= w;
			cur_flow -= w;
			ans += w;
			if(cur_flow <= 0) break;
		}
	}
	return ans;
}

// 最大流 
int Dinic(int s, int t){
	int ans = 0;
	while(bfs(s,t)){
		memset(cur, 0, sizeof(cur));
		ans += dfs(s,t,INF);
	}
	return ans;	
}

int main()
{
    freopen("in.txt","r",stdin);
    int N,F,D;
	while(scanf("%d%d%d",&N,&F,&D) == 3){
		for(int i = 1; i <= F; ++i) scanf("%d",&food[i]);
		for(int i = 1; i <= D; ++i) scanf("%d",&drink[i]);
		for(int i = 1; i <= N; ++i) scanf("%s",Gf[i]+1);
		for(int i = 1; i <= N; ++i) scanf("%s",Gd[i]+1);
		// 建图
		int s = 0, t = 2*N+F+D+1;
		init(2*N+F+D+2);
		for(int i = 1; i <= F; ++i) addEdge(s, i, food[i]);
		for(int i = 1; i <= N; ++i){
			addEdge(F+D+i, F+D+i+N, 1);
			for(int j = 1; j <= F; ++j)
				if(Gf[i][j] == 'Y') addEdge(j, F+D+i, INF);
			for(int j = 1; j <= D; ++j)
				if(Gd[i][j] == 'Y') addEdge(F+D+i+N, F+j, INF);
		}
		for(int i = 1; i <= D; ++i) addEdge(F+i, t, drink[i]);
		printf("%d\n", Dinic(s, t));
	}
    fclose(stdin);
	return 0;
}



  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值