飞行员配对方案[网络流]

传送门


直接网络流模板

注意输出方案,我们发现,但我们简单方向边不为0时,边的两头就是选了的


 代码

#include<bits/stdc++.h>
#define N 105
#define M N*N*4 
#define inf 0x3fffffff
using namespace std;
int first[N],next[M],to[M],w[M],tot=1;
int m,n,ed,st,x,y,dis[N],vis[N];
void add(int x,int y,int z){
	next[++tot]=first[x],first[x]=tot,to[tot]=y,w[tot]=z;
}
bool bfs(){
	memset(dis,0,sizeof(dis));
	memset(vis,0,sizeof(vis));
	queue<int> q; q.push(st); vis[st]=1;
	while(!q.empty()){
		int x=q.front(); q.pop();
		for(int i=first[x];i;i=next[i]){
			int t=to[i]; if(!w[i] || vis[t]) continue;
			dis[t] = dis[x] + 1; vis[t] = 1;
			if(t==ed) return true;
			q.push(t);
		}
	}return false;
}
int dfs(int u,int flow){
	if(u==ed) return flow;
	int ans=0;
	for(int i=first[u];i;i=next[i]){
		int t=to[i]; 
		if(dis[t] == dis[u] + 1){
			int delta = dfs(t , min(flow , w[i]));
			flow -= delta , ans += delta;
			w[i] -= delta , w[i^1] += delta;
			if(flow==0) break;
		}
	}
	if(flow) dis[u] = -1;
	return ans;
}
int dinic(){
	int ans=0;
	while(bfs()){
		if(int x=dfs(st,inf)) ans+=x;
	}return ans;
}
int main(){
	scanf("%d%d",&m,&n); st=0 , ed=n+1;
	for(int i=1;i<=m;i++) add(st,i,1),add(i,st,0);
	for(int i=m+1;i<=n;i++) add(i,ed,1),add(ed,i,0);
	while(1){
		scanf("%d%d",&x,&y); if(x==-1) break;
		add(x,y,1),add(y,x,0);
	}
	int x=dinic(); if(x==0){printf("No Solution");return 0;}
	printf("%d\n",x); 
	for(int i=2;i<=tot;i+=2)
		if(to[i]!=st && to[i^1]!=st && to[i]!=ed && to[i^1]!=ed && w[i^1])
			printf("%d %d\n",to[i^1],to[i]);
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值