POJ2536：Gopher II 匈牙利算法

题意：有n个地鼠和m个洞，有鹰飞来时，n个地鼠如果能在s秒从当前位置回到一个洞，就能不死，一个洞能容纳一个地鼠，它们的速度为v。求死的个数最小。

思路：匈牙利算法求最大匹配，n-最大匹配就是了,。。把能跑进洞里的老鼠和洞连起来。。

代码：

#include <iostream>
#include <math.h>
using namespace std;

#define MAX 110
int map[MAX][MAX];
int m,n,S,V,k;
int mk[MAX];
//从X集合中的顶点u出发用深度优先的策略寻找增广路
//(这种增广路只能使当前的匹配数增加1)
int nx,ny; //X和Y集合中顶点的个数
int cx[MAX],cy[MAX];
//cx[i]表示最终求得的最大匹配中与Xi匹配的Y顶点, cy[i]同理
int path(int u)
{
 for(int v=0; v<ny; v++) //考虑所有Yi顶点v
 {
  if(map[u][v]&&!mk[v])
  {
   mk[v]=1;
   //如果v没有匹配,或者如果v已经匹配了，
   //但从y[v]出发可以找到一条增广路
   if(cy[v]==-1|| path(cy[v]))
   {
    cx[u] = v; //把v匹配给u
    cy[v] = u; //把u匹配给v
    return 1; //找到可增广路
   }
  }
 }
 return 0 ; //如果不存在从u出发的增广路
}
int MaxMatch() //求二部图最大匹配的匈牙利算法
{
 int res=0;
 memset(cx,0xff,sizeof(cx)); //从0匹配开始增广
 memset(cy,0xff,sizeof(cy));
 for(int i=0; i<nx; i++)
 {
  if(cx[i]==-1) //从每个未盖点出发进行寻找增广路
  {
   memset(mk,0,sizeof(mk));
   res+=path(i); //每找到一条增广路，可使得匹配数加1
  }
 }
 return res;
}
int main()
{
 int i,j,k;
 double dx[101],dy[101];
 double dx2[101],dy2[101];
 while(cin>>n>>m>>S>>V)
 {
  nx=n;ny=m;
  memset(map,0,sizeof(map));
  for(i=0;i<n;i++)
  {
   cin>>dx[i]>>dy[i];
  }
  for(i=0;i<m;i++)
  {
   cin>>dx2[i]>>dy2[i];
  }
  for(i=0;i<n;i++)
   for(j=0;j<m;j++)
   {
    if(sqrt((dx[i]-dx2[j])*(dx[i]-dx2[j])+(dy[i]-dy2[j])*(dy[i]-dy2[j]))<=S*V)//这里又犯低级的错了。。
    {
     map[i][j]=1;//刚开始上面的dy2中的2没写进去。。导致了WA好几次。。
    }
   }
   int max=MaxMatch();
  // cout<<max<<endl;
   cout<<n-max<<endl;
 }
 return 0;
}