HDU 3605 Escape【二分图多重匹配】

这道题目和POJ2584很相似，我的博客也有详细的解释。点我传送(。・`ω´・)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<functional>
#include<algorithm>
using namespace std;
int maps[100005][20];//邻接矩阵
int vis[20];
int cy[20][100005];//点集Y的匹配情况
int numy[20],limit[20];//numy代表这个星球的匹配数量
int n,m;

int path(int u)//找增广路
{
	int i,j;
	for(i=0;i<m;i++)
	{
		if(maps[u][i]&&!vis[i])
		{
			vis[i]=1;
			if(numy[i]<limit[i])
			{
				cy[i][numy[i]++]=u;
				return 1;
			}
			else
			{
				for(j=0;j<limit[i];j++)
					if(path(cy[i][j]))
				{
					cy[i][j]=u;
					return 1;
				}
			}
		}
	}
	return 0;
}

int MulMatch()//二分图多重匹配
{
	int ans=0;
	memset(cy,0,sizeof(cy));
	memset(numy,0,sizeof(numy));
	for(int i=0;i<n;i++)
	{
		memset(vis,0,sizeof(vis));
		if(!path(i))
			return 0;
	}
	return 1;
//	if(ans==n)//开始这样写会超时，优化成上面就可以提前退出，就可以过
//		printf("YES\n");
//	else
//		printf("NO\n");
}
int main()
{
	int i,j,a;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		memset(maps,0,sizeof(maps));
		for(i=0;i<n;i++)
			for(j=0;j<m;j++)
				scanf("%d",&maps[i][j]);
		for(i=0;i<m;i++)
		scanf("%d",&limit[i]);
		if(MulMatch())
			printf("YES\n");
		else
			printf("NO\n");
	}
	return 0;
}