网络流24题——飞行员配对方案

 题目链接：https://www.luogu.org/problemnew/show/P2756

【问题分析】

二分图最大匹配问题。

【建模方法】

在二分图的基础上增加源S和汇T。
1、S向X集合中每个顶点连一条容量为1的有向边。
2、Y集合中每个顶点向T连一条容量为1的有向边。
3、XY集合之间的边都设为从A集合中的点到B集合之中的点，容量为1的有向边。

求网络最大流，流量就是匹配数，所有满流边是一组可行解。

【建模分析】

基本的二分图最大匹配，可以直接用匈牙利算法或Hopcroft_Karp算法解决，更一般的方法是网络最大流。

网络最大流写法

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 2000;
const int MAXM = 11000;
const int INF = 0x3f3f3f3f;
struct Edge1
{
	int from,to,cap,flow;
};
struct Dinic
{
	int n,m,s,t;
	vector<Edge1> edges;
	vector<int> G[MAXN];
	bool vis[MAXN];
	int d[MAXN];
	int cur[MAXN];
	void init(int n)
	{
		this -> n = n;
		for(int i = 0; i <= n + 1; i++){
			G[i].clear();
		}
		edges.clear();
	}
	void AddEdge(int from,int to,int cap)
	{
		edges.push_back((Edge1){from,to,cap,0});
		edges.push_back((Edge1){to,from,0,0});
		m = edges.size();
		G[from].push_back(m - 2);
		G[to].push_back(m - 1);
	}
	bool BFS()
	{
		memset(vis,0,sizeof(vis));
		queue<int> Q;
		Q.push(s);
		d[s] = 0;
		vis[s] = 1;
		while(!Q.empty()) {
			int x = Q.front();
			Q.pop();
			for(int i = 0; i < G[x].size(); i++) {
				Edge1& e = edges[G[x][i]];
				if(!vis[e.to] && e.cap > e.flow) {
					vis[e.to] = 1;
					d[e.to] = d[x] + 1;
					Q.push(e.to);
				}
			}
		}
		return vis[t];
	}
	int DFS(int x,int a)
	{
		if(x == t || a == 0) return a;
		int flow = 0,f;
		for(int& i = cur[x]; i < G[x].size(); i++) {
			Edge1& e = edges[G[x][i]];
			if(d[x] + 1 == d[e.to] && (f = DFS(e.to,min(a,e.cap - e.flow))) > 0) {
				e.flow += f;
				edges[G[x][i] ^ 1].flow -= f;
				flow += f;
				a -= f;
				if(a == 0) break;
			}
		}
		return flow;
	}
	int Maxflow(int s,int t) {
		this -> s = s,this -> t = t;
		int flow = 0;
		while(BFS()) {
			memset(cur,0,sizeof(cur));
			flow += DFS(s,INF);
		}
		return flow;
	}
}din;
int main(void)
{
    int m,n,u,v;
    int S,T;
    while(scanf("%d %d",&m,&n) != EOF) {
        din.init(n + 2);
        S = 0,T = n + 1;
        while((scanf("%d %d",&u,&v) != EOF) && !(u == -1 && v == -1 )) {
            din.AddEdge(u,v,1);
        }
        for(int i = 1; i <= m; i++) din.AddEdge(S,i,1);
        for(int i = m + 1; i <= n; i++) din.AddEdge(i,T,1);
        printf("%d\n",din.Maxflow(S,T));
        for(int i = 0; i < din.edges.size(); i++) {
            Edge1 e = din.edges[i];
            if(e.from != 0 && e.to != n + 1 && e.flow && e.from < e.to) {
                printf("%d %d\n",e.from,e.to);
            }
        }
    }
    return 0;
}

匈牙利找最大匹配写法

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 500;
int match[MAXN];
vector<int> G[MAXN];
bool used[MAXN];
int n,m;
bool find(int i)
{
	for(int j = 0; j < G[i].size(); j++) {
		int v = G[i][j];
		if(!used[v]) {
			used[v] = 1;
			if(match[v] == -1 || find(match[v])) {
				match[v] = i;
				return true;
			}
		}
	}
	return false;
}
int max_match()
{
	memset(match,-1,sizeof(match));
	int ans = 0;
	for(int i = 1; i <= m; i++) {
        memset(used,0,sizeof(used));
        if(find(i)) {
            ans++;
        }
	}
	return ans;
}
int main(void)
{
    int u,v;
    while(scanf("%d %d",&m,&n) != EOF ) {
        while(scanf("%d %d",&u,&v) != EOF && !(u == -1 && v == -1)) {
            G[u].push_back(v);
        }
        printf("%d\n",max_match());
        for(int i = m + 1; i <= n; i++) {
            if(match[i] != -1) {
                printf("%d %d\n",match[i],i);
            }
        }
    }
    return 0;
}