网络流二十四题 飞行员配对方案问题

P2756 飞行员配对方案问题

题目背景

第二次世界大战时期..

题目描述

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

对于给定的外籍飞行员与英国飞行员的配合情况，编程找出一个最佳飞行员配对方案，使皇家空军一次能派出最多的飞机。

输入输出格式

输入格式：

第 1 行有 2 个正整数 m 和 n。n 是皇家空军的飞行员总数(n<100)；m 是外籍飞行员数(m<=n)。外籍飞行员编号为 1~m；英国飞行员编号为 m+1~n。

接下来每行有 2 个正整数 i 和 j，表示外籍飞行员 i 可以和英国飞行员 j 配合。最后以 2个-1 结束。

输出格式：

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

输入输出样例

输入样例#1： 

5 10
1 7
1 8
2 6
2 9
2 10
3 7
3 8
4 7
4 8
5 10
-1 -1

输出样例#1： 

4
1 7
2 9
3 8
5 10 

---

思路

配对数量直接建立二分图套用dinic,数量为0输出No Solution!
输出路径的话，判断是否还有流量，即反边是否为真，连向s和t的边直接continue，
（s边和t边最后建的话，直接break）

样例都没过的代码却AC得代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
using namespace std;

int n,m,s,t;
struct node {
    int v,next,cap;
} e[6100];
int head[6100],tot=1;
int dis[6100];

void add_node(int a,int b,int c) {
    e[++tot]=(node) {b,head[a],c};
    head[a]=tot;
}

bool bfs() {
    memset(dis,-1,sizeof(dis));
    queue<int> q;
    q.push(s);
    dis[s]=0;
    while(!q.empty()) {
        int x=q.front();
        q.pop();
        for(int i=head[x]; i; i=e[i].next) {
            int v=e[i].v;
            if(e[i].cap&&dis[v]==-1) {
                dis[v]=dis[x]+1;
                q.push(v);
            }
        }
    }
    return dis[t]!=-1;
}

int dfs(int now,int f) {
    if(now==t) return f;
    int rest=f;
    for(int i=head[now]; i; i=e[i].next) {
        int v=e[i].v;
        if(e[i].cap>0&&rest>0&&dis[v]==dis[now]+1) {
            int t=dfs(v,min(rest,e[i].cap));
            if(!t) dis[v]=0;
            e[i].cap-=t;
            e[i^1].cap+=t;
            rest-=t;
        }
    }
    return f-rest;
}

void dinic() {
    int ans=0;
    while(bfs()) ans+=dfs(s,inf);
    if(!ans) {
        printf("No Solution!");
        return;
    }
    printf("%d
",ans);
    for(int i=2;i<=tot;i+=2)
    {
        int v=e[i].v;
        int v1=e[i^1].v;
        if(v==s||v1==s||v==t||v1==t) break;
        if(e[i^1].cap!=0)
        {
            printf("%d %d
",e[i^1].v,e[i].v);    
        }        
    }
}

int main() {
    scanf("%d%d",&m,&n);
    int x,y;
    while(scanf("%d%d",&x,&y)==2) {
        if(x==-1&&y==-1) break;
        add_node(x,y,1);
        add_node(y,x,0);
    }
    s=0;t=n+1;
    for(int i=1; i<=m; ++i) {
        add_node(s,i,1);
        add_node(i,s,0);
    }
    for(int i=m+1; i<=n; ++i) {
        add_node(i,t,1);
        add_node(t,i,0);
    }
    dinic();
}