【网络流】飞行员配对方案问题

题目描述
第二次世界大战时期， 英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的2 名飞行员， 其中1 名是英国飞行员，另1名是外籍飞行员。在众多的飞行员中，每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。 如何选择配对飞行的飞行员才能使一次派出最多的飞机。对于给定的外籍飞行员与英国飞行员的配合情况，试设计一个算法找出最佳飞行员配对方案，使皇家空军一次能派出最多的飞机。

输入
第1 行: 有 2个正整数m和 n。n是皇家空军的飞行员总数(n＜100)；m是外籍飞行员数。外籍飞行员编号为 1~m；英国飞行员编号为 m+1~n。 接下来每行有2个正整数 i 和 j ，表示外籍飞行员 i 可以和英国飞行员 j 配合。 文件最后以 2 个 -1 结束

输出
第 1 行是最佳飞行员配对方案一次能派出的最多的飞机数 M。 接下来 M 行是最佳飞行员配对方案。每行有 2个正整数 i 和 j ，表示在最佳飞行员配对方案中，飞行员 i 和飞行员 j 配对。 如果所求的最佳飞行员配对方案不存在，则输出‘No Solution!’ 。

样例输入
5 10
1 7
1 8
2 6
2 9
2 10
3 7
3 8
4 7
4 8
5 10
-1 -1
样例输出
4

这道题太裸了，源点连边到外国人， 外国人连边到飞行员， 飞行员连到汇点 。。。。。。

代码:

#include <cstdio>
#include <algorithm>
#include <cstring>

#define MAXN 1000
#define MAXM 500000
#define INF 1234567890
using namespace std;

struct node
{
    int v, c;
    node *next, *back;
};

//使用前要初始化
struct ISAP
{
    int n, s, t;
    node Edge[2*MAXM+10], *adj[MAXN+10], *ecnt;
    int d[MAXN+10];
    int num[MAXN+10];

    void init()
    {
        memset(adj, 0, sizeof adj);
        ecnt = Edge;
    }

    void init(int nn, int ss, int tt)
    {
        n = nn;
        s = ss;
        t = tt;
        memset(adj, 0, (n+1)*sizeof adj[0]);
        ecnt = Edge;
    }

    void addedge(int u, int v, int c)
    {
        node* p1 = ++ecnt;
        p1->v = v;
        p1->c = c;
        p1->next = adj[u];
        adj[u] = p1;

        node* p2 = ++ecnt;
        p2->v = u;
        p2->c = 0;
        p2->next = adj[v];
        adj[v] = p2;

        p1->back = p2;
        p2->back = p1;
    }

    int aug(int i, int augco)
    {
        int augc = 0;
        int mind = n-1;
        if (i == t)
            return augco;

        for (node* p = adj[i]; p; p=p->next)
        {
            int j = p->v, c = p->c;
            if (c > 0)
            {
                if (d[i] == d[j] + 1)
                {
                    int delta = min(c, augco - augc);
                    delta = aug(j, delta);
                    augc += delta;
                    p->c -= delta;
                    p->back->c += delta;
                    if (d[s] >= n)
                        return augc;
                    if (augc == augco)
                        return augc;
                }
                mind = min(mind, d[j]);
            }
        }

        if (augc == 0)
        {
            num[d[i]]--;
            if (num[d[i]] == 0) d[s] = n;
            d[i] = mind + 1;
            num[d[i]]++;
        }
        return augc;
    }

    int max_flow()
    {
       int f = 0;
       memset(d, 0, (n+1)*sizeof d[0]);
       memset(num, 0, (n+1)*sizeof num[0]);

       while (d[s] < n)
       {
           f += aug(s, INF);
       }

       return f;
    }
};

ISAP isap;

int GET()
{
    char ch;
    do{ch=getchar();}while(ch==' '||ch=='\n'||ch=='\t');
    int ret=0;
    while(ch>='0'&&ch<='9'){
        ret=ret*10+ch-'0';
        ch=getchar();}
    return ret;
}

void read()
{
    int m, n, u, v, c;
    scanf("%d%d", &m, &n);
    isap.init(n+2, 1, n+2);
    for (int i=1; i<=m; i++)
        isap.addedge(1, i+1, 1);
    while(scanf("%d %d", &u, &v)!=EOF && u != -1 && v != -1)
        isap.addedge(u+1, v+1, 1);
    for (int i=m+1; i<=n;i++)
        isap.addedge(i+1, n+2, 1);
}

int main()
{
    read();
    int ans = isap.max_flow();
    printf("%d", ans);
}