poj3281(#网络流最大流)

/*
translation:
    n头牛,f种食物,d种饮料。每头牛有各自喜欢的食物和饮料。而每种饮料和食物只能分配给一头牛。最多能够有多少头
    牛可以同时得到喜欢的食物和饮料。
solution:
    最大流
    为了保证一头牛只享有一份食物和饮料。必须把牛拆成两点。一边和对应喜欢的食物连接,一边和对应喜欢的饮料连接。
    然后每条边的容量是1.求最大流即可。
note:
    # 一开始想错了思路,以为用二部图匹配能够做出来,思路是这样的:将牛喜欢的各自食物和饮料连起来。然后求最多
      匹配。但是这样求出来的结果其实是食物和饮料的最大匹配,和最多能够分配多少头牛并没有什么联系。
*/
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>

using namespace std;
const int maxn = 500;
const int INF = 0x3f3f3f3f;

struct Edge
{
    int to, cap, rev;
    Edge(int to_, int cap_, int rev_):to(to_),cap(cap_),rev(rev_){}
};
vector<Edge> G[maxn];
vector<int> food[maxn];
vector<int> drink[maxn];
int level[maxn], iter[maxn];
int n, f, d;

void add_edge(int from, int to, int cap)
{
    G[from].push_back(Edge(to, cap, 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 > 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++) {
        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.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}

int dinic(int s, int t)
{
    int flow = 0;
    for(;;) {
        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()
{
    //freopen("in.txt", "r", stdin);
    while(~scanf("%d%d%d", &n, &f, &d)) {
        for(int i = 0; i < maxn; i++) {
            G[i].clear();
            food[i].clear();
            drink[i].clear();
        }

        int fn, dn, fi, di;
        for(int i = 0; i < n; i++) {
            scanf("%d%d", &fn, &dn);
            for(int j = 0; j < fn; j++) {
                scanf("%d", &fi);   fi--;
                food[i].push_back(fi);
            }
            for(int j = 0; j < dn; j++) {
                scanf("%d", &di);   di--;
                drink[i].push_back(di);
            }
        }

        for(int i = 0; i < n; i++) {
            for(int j = 0; j < food[i].size(); j++) {
                for(int k = 0; k < drink[i].size(); k++) {
                    add_edge(food[i][j], i + f + d, 1);
                    add_edge(i + f + d + n, drink[i][k] + f, 1);
                }
            }
        }

        int s = 2 * n + f + d, t = s + 1;
        for(int i = 0; i < f; i++)  add_edge(s, i, 1);
        for(int i = 0; i < d; i++)  add_edge(i + f, t, 1);
        for(int i = 0; i < n; i++)  add_edge(f + d + i, f + d + n + i, 1);

        printf("%d\n", dinic(s, t));
    }
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值