【网络流24题】【COGS14】飞行员搭配问题

【问题描述】

    飞行大队有若干个来自各地的驾驶员，专门驾驶一种型号的飞机，这种飞机每架有两个驾驶员,需一个正驾驶员和一个副驾驶员。由于种种原因，例如相互配合的问题，有些驾驶员不能在同一架飞机上飞行，问如何搭配驾驶员才能使出航的飞机最多。

 

如图，假设有10个驾驶员，如图中的V1，V2，…，V10就代表达10个驾驶员,其中V1，V2，V3，V4，V5是正驾驶员,V6，V7，V8，V9，V10是副驾驶员。如果一个正驾驶员和一个副驾驶员可以同机飞行，就在代表他们两个之间连一条线,两个人不能同机飞行，就不连。例如V1和V7可以同机飞行，而V1和V8就不行。请搭配飞行员，使出航的飞机最多。注意:因为驾驶工作分工严格,两个正驾驶员或两个副驾驶员都不能同机飞行.
 

【输入格式】

输入文件有若干行
第一行，两个整数n与n1，表示共有n个飞行员(2<=n<=100),其中有n1名飞行员是正驾驶员.
下面有若干行,每行有2个数字a,b。表示正驾驶员a和副驾驶员b可以同机飞行。

注:正驾驶员的编号在前,即正驾驶员的编号小于副驾驶员的编号.

【输出格式】

输出文件有一行
第一行，1个整数，表示最大起飞的飞机数。

【输入输出样例】

输入文件名： flyer.in

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

输出文件名：flyer.out

4 

裸的最大流啊= =

就是拿来练模板的

加一个超级源点超级汇点就好了

但是调了半天TAT果然太弱

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 21050
#define MAXINT 2100000000
using namespace std;
struct edge
{
	int f,c,v;
	edge *next,*rev;
	edge(){};
	edge(int ver,int cap,edge *nex):v(ver),c(cap),f(0),next(nex),rev(NULL){};
	void* operator new(size_t,void *p)
	{
		return p;
	}
}*s[MAXN];
int dis[MAXN],num[MAXN];
int n,m;
int src,des;
int u[10001],v[10001],w[10001];
int un;
int a[10001],b[10001];
int nu;
void init()
{
	int queue[MAXN],head=0,tail=0;
	for (int i=0;i<=des;i++)
	    dis[i]=MAXN,num[i]=0;
	queue[tail++]=des;
	dis[des]=0;
	num[0]=1;
	while (head!=tail)
	{
		int v=queue[head++];
		for (edge *e=s[v];e;e=e->next)
		{
			if ((e->rev&&e->rev->c==0)||dis[e->v]<MAXN) continue;
			dis[e->v]=dis[v]+1;
			++num[dis[e->v]];
			queue[tail++]=e->v;
		}
	}
}
int maxflow()
{
	int st=src,ret=0;
	edge *E[MAXN],*rep[MAXN];
	for (int i=0;i<=n+1;i++) E[i]=s[i];
	while (dis[src]<des)
	{
		if (st==des)
		{
			int delta=MAXINT;
 			for (int i=src;i!=des;i=E[i]->v)
			    delta=min(delta,E[i]->c);
			for (int i=src;i!=des;i=E[i]->v)
			{
				E[i]->c-=delta;
				E[i]->rev->c+=delta;
				E[i]->f+=delta;
				E[i]->rev->f-=delta;
			}
			ret+=delta;
			st=src;
		}
		edge *e;
		for (e=E[st];e;e=e->next)
		    if (e->c>0&&dis[st]==dis[e->v]+1) break;
		if (e)
		{
			E[st]=e;
			rep[e->v]=e->rev;
			st=e->v;
		}
		else
		{
			if ((--num[dis[st]])==0) break;
			E[st]=s[st];    
			int mind=n;
			for (edge *t=s[st];t;t=t->next)
			    if (t->c>0)
			        mind=min(mind,dis[t->v]);
			dis[st]=mind+1;
			++num[dis[st]];
			if (st!=src) st=rep[st]->v;
		}
	}
	return ret;
}
int main()
{
	freopen("flyer.in","r",stdin);
	freopen("flyer.out","w",stdout);
	scanf("%d%d",&n,&m);
	src=0,des=n+1;
	while (~scanf("%d%d",&u[un],&v[un]))w[un++]=1;
	u[un]=0;v[un]=0;
	for (int i=1;i<=m;i++)
	    u[un]=0,v[un]=i,w[un++]=1;
	for (int i=m+1;i<=n;i++)
	    u[un]=i,v[un]=n+1,w[un++]=1;
	edge *buffer=new edge[un*2];
	edge *data=buffer;
	for (int i=0;i<un;i++)
	{
		s[u[i]]=new((void*) data++)edge(v[i],w[i],s[u[i]]);
		s[v[i]]=new((void*) data++)edge(u[i],0,s[v[i]]);
		s[u[i]]->rev=s[v[i]];
		s[v[i]]->rev=s[u[i]];
	}
	init();
	printf("%d",maxflow());
	
}