「网络流 24 题」[1] 飞行员配对方案问题

题意

有m个外籍飞行员和n个英国飞行员，中间存在两两可以合作的关系，一架飞机需要一个外籍飞行员和一个英国飞行员。

求最多能飞出去多少架飞机，并且输出任意一种方案。

https://www.oj.swust.edu.cn/problem/show/1736

(似乎因为这个有spj...)

(洛谷上也有spj啊...)

(LOJ上也有spj啊...)

(打出两个大凶好烦啊...)

分析

因为是第一次敲网络流所以说还是说一下怎么整...

简单来说就是对于一个有向图，边上会有一个容量限制，表示这条边最多能流多少，然后求出最后某个点能流到另外一个点多少流量。

然后学习的是dinic算法。

模板→ http://note.youdao.com/noteshare?id=eff621a0649ed0b0c32ab9457e52b20b&sub=A3CF0D061D5A40F4B333B0AB74B64EF4

按照大佬们的说法，网络流的问题困难的地方主要在建图。

但是对于这道题来说还是比较好解决的...

按照二分图原来的样子建图，所有的边的容量都是1，然后建一个源点s，连向二分图中左边所有的点，建一个汇点t，将二分图右边所有的点连向

t。然后求s到t的最大流就好了。

code

#include<bits/stdc++.h>
#define inf 100000
using namespace std;
#define M 605
void read(int &x){
	x=0; char c=getchar(); int f=1;
	for (;c<48;c=getchar())if (c=='-')f=-1;
	for (;c>47;c=getchar())x=(x<<1)+(x<<3)+(c^48);
	x*=f;
}
struct ed{
	int x,cap,nx;
}e[M*M];
int ecnt,nx[M];

void add(int x,int y,int cap){
	e[ecnt]=(ed){y,cap,nx[x]};
	nx[x]=ecnt++;
	e[ecnt]=(ed){x,0,nx[y]};
	nx[y]=ecnt++;
}

struct Dinic{
	int nnx[M];
	int Q[M],level[M],s,t;
	bool bfs(int x){
		memset(level,0,sizeof(level));
		int l=0,r=0,i;
		level[x]=1;
		Q[r++]=x;
		for (;l<r;){
			x=Q[l++];
			for (i=nx[x];~i;i=e[i].nx){
				if (e[i].cap>0&&!level[e[i].x]){
					level[e[i].x]=level[x]+1;
					Q[r++]=e[i].x;
					if (e[i].x==t)return 1;
				}
			}
		}
		return level[t]>0;
	}
	int dfs(int x,int f){
		if (x==t)return f;
		int sum=0;
		for (int &i=nnx[x];~i;i=e[i].nx)if (e[i].cap>0&&level[x]+1==level[e[i].x]){
			int d=dfs(e[i].x,min(e[i].cap,f-sum));
			e[i].cap-=d;
			e[i^1].cap+=d;
			sum+=d;
			if (sum==f)return f;
		}
		if (sum==0)level[x]=0;
		return sum;
	}
	int max_flow(int S,int T){
		s=S; t=T;
		int res=0;
		for (;bfs(s);){
			memcpy(nnx,nx,sizeof(nx));
			res+=dfs(s,inf); 
		}
		return res;
	}
}dinic;

int main(){
//	freopen("1.txt","r",stdin);
	int m,n,i,x,y;
	read(n); read(m);
	memset(nx,-1,sizeof(nx));
	for(i=1;i<=m;i++)add(0,i,1);
	for(i=m+1;i<=n+m;i++)add(i,n+m+1,1);
	for(;~scanf("%d%d",&x,&y);)add(x,y,1);
	int res=dinic.max_flow(0,n+m+1);
	printf("%d\n",res);
	return 0;
}