P2756 飞行员配对方案问题

(color{#0066ff}{题目描述})

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

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

(color{#0066ff}{输入格式})

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

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

(color{#0066ff}{输出格式})

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

(color{#0066ff}{输入样例})

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

(color{#0066ff}{输出样例})

4
1 7
2 9
3 8
5 10 

(color{#0066ff}{数据范围与提示})

有spj

(color{#0066ff}{题解})

S向每个外籍飞行员连容量为1的边

题中给出的边，容量为1

每个英国飞行员向T连容量为1的边

跑Dinic

最后找边即可

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#include<algorithm>
#include<cmath>
#define _ 0
#define LL long long
inline LL in()
{
	LL x=0,f=1; char ch;
	while(!isdigit(ch=getchar()))(ch=='-')&&(f=-f);
	while(isdigit(ch)) x=x*10+(ch^48),ch=getchar();
	return x*f;
}
const int inf=0x7fffffff;
int n,m,s,t;
struct node{int dis,to,nxt;}e[105050];
int cnt=1;
std::queue<int> q;
int head[150],cur[150],dis[150],dep[150];
inline void add(int from,int to,int dis)
{
	cnt++;
	e[cnt].to=to;
	e[cnt].dis=dis;
	e[cnt].nxt=head[from];
	head[from]=cnt;
}
inline bool bfs()
{
	for(int i=0;i<=n+1;i++) dep[i]=0,cur[i]=head[i];
	q.push(s);
	dep[s]=1;
	while(!q.empty())
	{
		int tp=q.front();q.pop();
		for(int i=head[tp];i;i=e[i].nxt)
		{
			int go=e[i].to;
			if(!dep[go]&&e[i].dis>0)
			{
				dep[go]=dep[tp]+1;
				q.push(go);
			} 
		}
	}
	return dep[t];
}
inline int dfs(int x,int change)
{
	if(x==t||!change) return change;
	int flow=0,ls;
	for(int i=cur[x];i;i=e[i].nxt)
	{
		int go=e[i].to;
		cur[x]=i;
		if(dep[go]==dep[x]+1&&(ls=dfs(go,std::min(e[i].dis,change))))
		{
			change-=ls;
			flow+=ls;
			e[i].dis-=ls;
			e[i^1].dis+=ls;
			if(!change) break;
		}
	}
	return flow;
}
inline void maxflow()
{
	int flow=0;
	while(bfs()) flow+=dfs(s,inf);
	if(!flow) return (void)(printf("No Solution!"));
	printf("%d
",flow);
	for(int i=1;i<=m;i++)
	{
		for(int j=head[i];j;j=e[j].nxt)
		{
			int go=e[j].to;
			if(go>=m+1&&go<=n&&!e[j].dis) 
			{
				printf("%d %d
",i,go);
				break;
			}
		}
	}
}
int main()
{
	m=in(),n=in();
	s=0,t=n+1;
	int x,y;
	while(1)
	{
		x=in(),y=in();
		if(!(~x||~y)) break;
		add(x,y,1),add(y,x,0);
	}
	for(int i=1;i<=m;i++) add(s,i,1),add(i,s,0);
	for(int i=m+1;i<=n;i++)  add(i,t,1),add(t,i,0);
	maxflow();
	return 0;
}