POJ 3041 Asteroids(二分匹配模板题)

/*
poj 3041
行和列建立二分图
一个点上有东西则建立一条边。
题目就相当于求最小点覆盖数=最大二分匹配数
*/

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;


/* **************************************************************************
//二分图匹配(匈牙利算法的DFS实现)
//初始化:g[][]两边顶点的划分情况
//建立g[i][j]表示i->j的有向边就可以了,是左边向右边的匹配
//g没有边相连则初始化为0
//uN是匹配左边的顶点数,vN是匹配右边的顶点数
//调用:res=hungary();输出最大匹配数
//优点:适用于稠密图,DFS找增广路,实现简洁易于理解
//时间复杂度:O(VE)
//***************************************************************************/
//顶点编号从0开始的
const int MAXN=510;
int uN,vN;//u,v数目
int g[MAXN][MAXN];
int linker[MAXN];
bool used[MAXN];
bool dfs(int u)//从左边开始找增广路径
{
 int v;
 for(v=0;v<vN;v++)//这个顶点编号从0开始,若要从1开始需要修改
 if(g[u][v]&&!used[v])
 {
 used[v]=true;
 if(linker[v]==-1||dfs(linker[v]))
 {//找增广路,反向
 linker[v]=u;
 return true;
 }
 }
 return false;//这个不要忘了,经常忘记这句
}
int hungary()
{
 int res=0;
 int u;
 memset(linker,-1,sizeof(linker));
 for(u=0;u<uN;u++)
 {
 memset(used,0,sizeof(used));
 if(dfs(u)) res++;
 }
 return res;
}
//******************************************************************************/



int main()
{
 int k;
 int n;
 int u,v;
 while(scanf("%d%d",&n,&k)!=EOF)
 {
 uN=vN=n;
 memset(g,0,sizeof(g));
 while(k--)
 {
 scanf("%d%d",&u,&v);
 u--;
 v--;
 g[u][v]=1;
 }
 printf("%d\n",hungary());
 }
 return 0;
}

你可能感兴趣的:(二分匹配)