51Nod-2006-飞行员配对（二分图最大匹配）

ACM模版

描述

题解

这个题的来源是网络流24题，貌似是一个十分不错的网络流习题集，暑假抽空做做吧！

越学习，越发现自己是一个大大的菜鸡，心痛啊~~~我网络流都不会……

这个题很明显是一个二分图的最大匹配问题（不要问我怎么看出来的），所以既可以用匈牙利算法解（代码 One），也可以使用构建网络流模型来解（代码 Two）。匈牙利的问题就不多说了，裸算法，（0级题都是裸算法模版题），而这里如果用网络流模型来解，就涉及到了一个建图的问题，首先我们先建立一个超级源点通向外籍飞行员，权值为 $1$，接着建立一个超级汇点，由所有英国飞行员通向它，权值同样为 $1$，剩下的就是外籍和英国飞行员之间的连通关系，权值一样都是 1 <script type="math/tex" id="MathJax-Element-3">1</script>，然后此图的最大流，即为二分图的最大匹配数，很有趣的题，这里使用的是 Dinic 算法来求最大流。

增广路相关算法真是有趣至极，今天看了好久，总算是有所入门了，待我暑假，网络流专题走起！！！

代码

One：

//  匈牙利算法
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

const int MAXN = 105;
const int MAXE = 1e5 + 10;

struct edge
{
    int v;
    int nt;
};

int n, m;
int ans, tot;
int head[MAXN];
int match[MAXN];
bool flag[MAXN];
edge map[MAXE];

void add_edge(int x, int y)
{
    tot++;
    map[tot].v = y;
    map[tot].nt = head[x];
    head[x] = tot;
}

void init()
{
    memset(head, 0, sizeof(head));
    memset(match, 0, sizeof(match));

    ans = tot = 0;
}

bool dfs(int x)
{
    for (int e = head[x]; e; e = map[e].nt)
    {
        int v = map[e].v;
        if (!flag[v])
        {
            flag[v] = true;
            if (match[v] == 0 || dfs(match[v]))
            {
                match[v] = x;
                return true;
            }
        }
    }
    return false;
}

void hungary()
{
    for (int i = 1; i <= m; i++)
    {
        memset(flag, false, sizeof(flag));
        if (dfs(i))
        {
            ans++;
        }
    }
}

int main()
{
    init();

    scanf("%d%d", &m, &n);

    int x, y;
    while (scanf("%d%d", &x, &y), x != -1  &&  y != -1)
    {
        add_edge(x, y);
    }

    hungary();

    printf("%d\n", ans);

    return 0;
}

Two：

//  网络流模型
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <stack>

using namespace std;

const int MAXN = 105;
const int MAXE = 1e5 + 10;
const int INF = 0x3f3f3f3f;

struct edge
{
    int v;
    int next, flow, op; //  opposite
};

int n, m;
int ans, tot;
int st, ed;     //  超级源点与超级汇点
int dis[MAXN];
int head[MAXN];
int path[MAXN];
bool vis[MAXN];
edge map[MAXE];
queue<int> que;
stack<int> sta;

void add_edge(int x, int y, int flow)
{
    tot++;
    map[tot].v = y;
    map[tot].flow = flow;
    map[tot].next = head[x];
    head[x] = tot;
    map[tot].op = tot + 1;
    tot++;
    map[tot].v = x;
    map[tot].flow = 0;
    map[tot].next = head[y];
    head[y] = tot;
    map[tot].op = tot - 1;
}

void init()
{
    ans = tot = 0;
    memset(head, 0, sizeof(head));

    st = 0;
    ed = n + 1;
    for (int i = 1; i <= m; i++)
    {
        add_edge(st, i, 1);
    }
    for (int i = m + 1; i <= n; i++)
    {
        add_edge(i, ed, 1);
    }
}

bool SPFA()
{
    memset(vis, false, sizeof(vis));
    memset(dis, 0x3f, sizeof(dis));

    vis[st] = true;
    dis[st] = 0;
    while (!que.empty())
    {
        que.pop();
    }

    que.push(st);

    while (!que.empty())
    {
        int u = que.front();
        que.pop();
        vis[u] = false;

        for (int t = head[u]; t; t = map[t].next)
        {
            int v = map[t].v;
            if (dis[u] + 1 < dis[v] && map[t].flow)
            {
                dis[v] = dis[u] + 1;
                if (!vis[v])
                {
                    vis[v] = true;
                    que.push(v);
                }
            }
        }
    }

    if (dis[ed] < INF)
    {
        return true;
    }
    return false;
}

void Dinic()
{
    while (!sta.empty())
    {
        sta.pop();
    }
    sta.push(st);

    int x, last = 1;
    while (last)
    {
        while (last != sta.size())
        {
            sta.pop();
        }
        x = sta.top();
        if (x != ed)
        {
            int t = head[x];
            while (t && (dis[x] + 1 != dis[map[t].v] || map[t].flow == 0))
            {
                t = map[t].next;
            }
            if (t == 0)
            {
                dis[x] = INF;
                sta.pop();
                last--;
            }
            else
            {
                path[last] = t;
                sta.push(map[t].v);
                last++;
            }
        }
        else
        {
            int min_flow = INF;
            for (int i = 1; i < last; i++)
            {
                if (map[path[i]].flow < min_flow)
                {
                    min_flow = map[path[i]].flow;
                }
            }
            ans += min_flow;

            int tmp = 0;
            for (int i = last - 1; i > 0; i--)
            {
                map[path[i]].flow -= min_flow;
                map[map[path[i]].op].flow += min_flow;
                if (map[path[i]].flow == 0)
                {
                    tmp = i;
                }
            }
            last = tmp;
        }
    }
}

void solve()
{
    while (SPFA())
    {
        Dinic();
    }
}

int main()
{
    scanf("%d%d", &m, &n);

    init();

    int x, y;
    while (scanf("%d%d", &x, &y), x != -1 && y != -1)
    {
        add_edge(x, y, 1);
    }

    solve();

    if (ans == 0)
    {
        printf("No Solution!\n");
    }
    else
    {
        printf("%d\n", ans);
    }

    return 0;
}