poj 2536 Gopher II 二分匹配应用 匈牙利算法

题意：有某种动物n个，以及m个洞穴。给定n个动物坐标，以及m个洞穴坐标。每个动物有相同的运动速度v，单位是m/s。在老鹰来时，如果不能在s秒内（含s秒）躲进洞里，将会被吃，求怎么分配动物及洞穴，使得尽量少的动物被抓。

思路：建立结点，左边为动物，右边为洞穴，对于左边任意一个动物，若能够在s秒内躲进右边任意一个洞穴，则之间连一条边。接下来就是求最大匹配，用总数减去最大匹配数，就是最少被抓数量。

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<iostream>
using namespace std;
const int maxn = 400;
struct edge
{
	int from,to;
};
struct Pt
{
	double x,y;
}pt[maxn];
vector<edge> edges;
vector<int> g[maxn];
int match[maxn],check[maxn];
int n,m,s,v;
bool dfs(int u)
{
	for(int i=0;i<g[u].size();i++)
	{
		int v=edges[g[u][i]].to;
		if(!check[v])
		{
			check[v]=1;
			if(match[v]==-1||dfs(match[v]))
			{
				match[u]=v;
				match[v]=u;
				return true;
			}
		}
		
	}
	return false;
}
int Hungarian()
{
	int ans=0;
	memset(match,-1,sizeof(match));
	for(int i=1;i<=n;i++)
	{
		if(match[i]==-1)
		{
			memset(check,0,sizeof(check));
			if(dfs(i))ans++;
		}
		
	}
	return ans;
}
double dist(int i,int j)
{
	return sqrt((pt[i].x-pt[j].x)*(pt[i].x-pt[j].x)+
	(pt[i].y-pt[j].y)*(pt[i].y-pt[j].y));
}
void addedge(int from,int to)
{
	edges.push_back(edge{from,to});
	edges.push_back(edge{to,from});
	int m=edges.size();
	g[from].push_back(m-2);
	g[to].push_back(m-1);
}
int main()
{
	//freopen("in.txt","r",stdin);
	while(scanf("%d%d%d%d",&n,&m,&s,&v)!=EOF)
	{
		double x,y;
		edges.clear();
		for(int i=1;i<=n+m;i++)g[i].clear();
		int i=1;
		for(;i<=n;i++)
		{
			cin>>x>>y;
			pt[i].x=x;
			pt[i].y=y;
		}
		m=n+m;
		for(i=n+1;i<=m;i++)
		{
			cin>>x>>y;
			pt[i].x=x;
			pt[i].y=y;
		}
		for(int i=1;i<=n;i++)
		{
			for(int j=n+1;j<=m;j++)
			{
				double dt=dist(i,j)/(v*1.0);
				if(dt<=s)
				{
					addedge(i,j);
				}
			}
		}
		int ans=Hungarian();
		printf("%d\n",n-ans);
	}
}