洛谷 p2756 飞行员配对方案问题（最大流，二分图匹配）

题目

题目背景
第二次世界大战期间，英国皇家空军从沦陷国征募了大量外籍飞行员。由皇家空军派出的每一架飞机都需要配备在航行技能和语言上能互相配合的两名飞行员，其中一名是英国飞行员，另一名是外籍飞行员。在众多的飞行员中，每一名外籍飞行员都可以与其他若干名英国飞行员很好地配合。

题目描述
一共有 nn 个飞行员，其中有 mm 个外籍飞行员和 (n - m)(n−m) 个英国飞行员，外籍飞行员从 11 到 mm 编号，英国飞行员从 m + 1m+1 到 nn 编号。 对于给定的外籍飞行员与英国飞行员的配合情况，试设计一个算法找出最佳飞行员配对方案，使皇家空军一次能派出最多的飞机。

输入格式
输入的第一行是用空格隔开的两个正整数，分别代表外籍飞行员的个数 mm 和飞行员总数 nn。
从第二行起到倒数第二行，每行有两个整数 u, vu,v，代表外籍飞行员 uu 可以和英国飞行员 vv 配合。
输入的最后一行保证为 -1 -1，代表输入结束。

输出格式
本题存在 Special Judge。
请输出能派出最多的飞机数量，并给出一种可行的方案。
输出的第一行是一个整数，代表一次能派出的最多飞机数量，设这个整数是 kk。
第 22 行到第 k + 1k+1 行，每行输出两个整数 u, vu,v，代表在你给出的方案中，外籍飞行员 uu 和英国飞行员 vv 配合。这 kk 行的 uu 与 vv 应该互不相同。

思路

二分图最大匹配，我们用dinic跑一边就好了。

代码实现

#include<cstdio>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
#define rep(i,f_start,f_end) for (int i=f_start;i<=f_end;++i)
#define per(i,n,a) for (int i=n;i>=a;i--)
#define MT(x,i) memset(x,i,sizeof(x) )
#define rev(i,start,end) for (int i=start;i<end;i++)
#define inf 0x3f3f3f3f
#define mp(x,y) make_pair(x,y)
#define lowbit(x) (x&-x)
#define MOD 1000000007
#define exp 1e-8
#define N 1000005 
#define fi first 
#define se second
#define pb push_back
typedef long long ll;
typedef vector <int> VI;
typedef pair<int ,int> PII;
typedef pair<int ,PII> PIII;
ll gcd (ll a,ll b) {return b?gcd (b,a%b):a; }
inline int read() {
    char ch=getchar(); int x=0, f=1;
    while(ch<'0'||ch>'9') {
        if(ch=='-') f = -1;
        ch=getchar();
    } 
    while('0'<=ch&&ch<='9') {
        x=x*10+ch-'0';
        ch=getchar();
    }   return x*f;
}

const int maxn=1e5+10;
struct edge {
    int v,flow,next;
}e[maxn];
int head[maxn],cnt=-1;

inline void add (int u,int v,int flow) {
    e[++cnt]=(edge) {v,flow,head[u]};
    head[u]=cnt;
}
inline void add_edge (int u,int v,int flow) {
    add (u,v,flow);
    add (v,u,0);
}

int dis[maxn],cur[maxn];
int s,t,n,m,ans,x,y;
int bfs () {
    MT (dis,0);
    queue<int > q;
    q.push (s);
    dis[s]=1;
    while (q.size ()) {
        int x=q.front (); q.pop ();
        for (int i=head[x];i!=-1;i=e[i].next) {
            if (dis[e[i].v]==0&&e[i].flow) {
                dis[e[i].v]=dis[x]+1;
                q.push (e[i].v);
            }
        }
    }
    return dis[t];
}

int dfs (int now,int nowflow) {
    if (now==t) return nowflow;
    for (int &i=cur[now];~i;i=e[i].next) {
        if (dis[e[i].v]==dis[now]+1&&e[i].flow) {
            int canflow=dfs (e[i].v,min (nowflow,e[i].flow));
            if (canflow>0) {
                e[i].flow-=canflow;
                e[i^1].flow+=canflow;
                return canflow;
            }
        }
    }
    return 0;
}

void Dinic () {
    while (bfs ()) {
       for (int i=s;i<=t;i++) cur[i]=head[i];
       while (int val=dfs(s,inf)) ans+=val;
    }
}
int main () {
    cin>>m>>n;
    MT (head,-1);
    s=0,t=n+1;
    rep (i,1,m) add_edge (s,i,1);
    rep (i,m+1,n) add_edge (i,t,1);

    while (scanf ("%d%d",&x,&y)&&x!=-1&&y!=-1) add_edge (x,y,1);
    Dinic ();
    if (!ans) {
        cout<<"no answer"<<endl;
    }
    else {
        cout<<ans<<endl;
        rep (i,1,m) 
            for (int j=head[i];~j;j=e[j].next) 
                if (e[j].flow==0&&e[j^1].flow==1&&i!=s&&i!=t&&e[j].v!=s&&e[j].v!=t) 
                   printf ("%d %d
",i,e[j].v); 
    }
    return 0;
}