PowerOJ1736: 飞行员配对方案问题

Description

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

Input

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

程序运行结束时，将最佳飞行员配对方案输出到文件output.txt 中。第1 行是最佳飞行
员配对方案一次能派出的最多的飞机数M。接下来M 行是最佳飞行员配对方案。每行有2
个正整数i和j，表示在最佳飞行员配对方案中，飞行员i和飞行员j 配对。
如果所求的最佳飞行员配对方案不存在，则输出‘No Solution!’。
Sample Input
Raw
5 10
1 7
1 8
2 6
2 9
2 10
3 7
3 8
4 7
4 8
5 10
-1 -1
Sample Output
Raw
4
1 7
2 9
3 8
5 10
Hint

王晓东《线性规划和网络流24题》系列题目不需要管output.txt，是标准输入输出，多文件单数据。
所有题目中n的值一般小于150.

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <stack>
#include <queue>
#include <vector>
using namespace std;

const int INF = 0x3f3f3f3f;
const int maxn = 200+2;
const int maxm = 400000+2;

int head[maxn],level[maxn],fa[maxn],no,s,e;
int m,n;
struct Point{
    int u,v,flow,next;
    Point(){};
    Point(int sta,int end,int z,int w):u(sta),v(end),next(z),flow(w){};
}p[maxm];

inline void add(int x,int y,int z){
    p[no] = Point(x,y,head[x],z);
    head[x] = no++;
    p[no] = Point(y,x,head[y],0);
    head[y] = no++;
}

inline void Init(){
    memset(head,-1,sizeof(head));
    memset(fa,-1,sizeof(fa));
    no = 0;
    scanf("%d %d",&m,&n);
    int readi,readj;
    while(scanf("%d %d",&readi,&readj))
        if(readi == -1 && readj == -1 )
            break;
        else
            add(readi,readj,1);             //建图

    s = 0;
    e = n + 1;
    for(int i = 1;i <= n;++i){
        if(i <= m)
            add(s,i,1);
        else 
            add(i,e,1);
    }
    return;
}

bool bfs(int s,int e){
    int tempx,tempy;
    memset(level,0,sizeof(level));
    level[s] = 1;
    queue <int> q;
    q.push(s);
    while(!q.empty()){
        tempx = q.front();q.pop();
        if(tempx == e)
            return true;
        for(int i = head[tempx];i != -1; i = p[i].next){
            if(p[i].flow && !level[tempy = p[i].v]){
                level[tempy] = level[tempx] + 1;
                q.push(tempy);
            }
        }
    }
    return false;
}

int dfs(int u,int maxf,int e){
    if(u == e)
        return maxf;
    int ret = 0;
    for(int i = head[u];i != -1; i = p[i].next){
        int v = p[i].v,f = p[i].flow;
        if(level[u] + 1 == level[v] && f){
            int Min = min(maxf - ret,f);
            f = dfs(v,Min,e);
            if(f) 
                fa[u] = p[i].v;
            p[i].flow -= f;
            p[i^1].flow += f;   
            ret += f;
            if(ret == maxf)
                return ret;
        }
    }
    return ret;
}

inline void Dinic(){
    int ans = 0;
    while(bfs(s,e))
        ans += dfs(s,INF,e);
    if(ans){
        printf("%d\n",ans);
        for(int i = 1;i <= m;++i)
            if(fa[i] != -1 && fa[i] != s)
                printf("%d %d\n",i,fa[i]);
    }
    else
        printf("No Solution!\n");
    return;
}

int main()
{
    Init();
    Dinic();
    return 0;
}