匈牙利算法

匈牙利算法，n平方解决二分图匹配问题

【链式前向星好啊】

例（shui）题：

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

说是网络流怎么可以用匈牙利呢。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define t (dis[i])
using namespace std;
const int maxn=110,maxm=510<<1;
void init();
void work();
bool dfs(int);
int main() {
	init(),work();
	return 0;
}
inline int gi();
int n,n1,m,fir[maxn],nxt[maxm],dis[maxm];
inline void adde(int a,int b,int id) {
	nxt[id]=fir[a],fir[a]=id,dis[id]=b;
}
bool vis[maxn];
int T[maxn];
void work() {
	memset(T,-1,sizeof(T));
	int ans=0;
	for(int i=1; i<=n1; i++) {
		memset(vis,0,sizeof(vis));
		ans+=dfs(i);
	}
	printf("%d\n",ans);
}
void init() {
	scanf("%d%d",&n,&n1);
	int a,b;
	while(scanf("%d%d",&a,&b)==2)adde(a,b,++m),adde(b,a,++m);
}
bool dfs(int now) {
	for(int i=fir[now]; i; i=nxt[i])
		if(!vis[t]) {
			vis[t]=1;
			if(T[t]==-1||dfs(T[t])){
				T[t]=now;return 1;
			}
		}
	return 0;
}

顺便再贴上【网络流24题】】【水题】骑士共存的代码（网络流？？？为什么要用网络流）

#include<cstdio>
#include<algorithm>
#include<cstring>
#define check(x,y) (x>0&&x<=n&&y>0&&y<=n)
using namespace std;
const int X[]= {99999,-1,1,2,2};
const int Y[]= {888888888,2,2,1,-1};
int id[201][201],fir[40001],nxt[40001*8*2],dis[40001*8*2];
inline int gi() {
	int x=0,f=1;
	char ch=getchar();
	while(ch>'9'||ch<'0')f=(ch=='-')?-1:f,ch=getchar();
	while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
	return x*f;
}
int Id=0;
inline void adde(int a,int b) {
	nxt[++Id]=fir[a],fir[a]=Id,dis[Id]=b;
}
int N=0,n,m;
inline void init() {
	int a,b;
	n=gi(),m=gi();
	for(int i=1; i<=m; i++)a=gi(),b=gi(),id[a][b]=-1;
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			if(id[i][j]!=-1)id[i][j]=++N;
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)if(id[i][j]>0)
			for(int k=1; k<5; k++)
				if(check((i+X[k]),(j+Y[k]))&&id[i+X[k]][j+Y[k]]>0) {
					adde(id[i][j],id[i+X[k]][j+Y[k]]);
					adde(id[i+X[k]][j+Y[k]],id[i][j]);
				}
}
int match[40001],vis[40001];
inline bool dfs(int now) {
	for(int i=fir[now]; i; i=nxt[i])
		if(!vis[dis[i]]) {
			int t=dis[i];
			vis[t]=1;
			if(match[t]==-1||dfs(match[t])) {
				match[t]=now;
				return 1;
			}
		}
	return 0;
}
inline void hun() {
	int ans=0;
	memset(match,-1,sizeof match);
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++)
			if((i+j)&1) {
				memset(vis,0,sizeof vis);
				ans+=dfs(id[i][j]);
			}
	printf("%d",N-ans);
}
int main() {
	init(),hun();
	return 0;
}