[ssl 1333] 地鼠的困境

题目

Description
　　地鼠家族面临着一个新的威胁——猎食者。
　　地鼠家族一共有N个地鼠和M个鼠洞，每个都位于不同的(x, y)坐标中。假如有地鼠在发觉危险以后s秒内都没有回到鼠洞里的话，就可能成为老鹰的食物。当然了，一个鼠洞只能拯救一只地鼠的命运，所有地鼠都以相等的速度v移动。地鼠家族需要设计一种策略，使得老鹰来时，易受攻击的地鼠数量最少。

Input
本题有多组数据。第1行为测试数据组数T(T<=50)。
对于每组数据，第一行4个整数n, m, s和v(n, m <= 100)。以后n行为地鼠的坐标，以后m行为鼠洞的坐标。距离的单位是m，时间的单位是s，速度的单位是m/s。

Output
对于每组数据输出一行，为易受攻击的地鼠的数量。

Sample Input
1
2 2 5 10
1.0 1.0
2.0 2.0
100.0 100.0
20.0 20.0

Sample Output
1

---

解题思路

匈牙利算法，求其距离小于s*v的最小匹配值

---

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
struct node{
    int x,y,next;
}a[10001];
int ans,add,ww,n,m,s,v,link[101],list[101],vt[101];
double w[101],e[101],w1,w2;
double cy(double c) { return c*c;}
bool find(int t)
{
   int g=list[t]; 
   while (g){
    if (!vt[a[g].y]){
        vt[a[g].y]=1; 
        int u=link[a[g].y]; 
        link[a[g].y]=t; 
        if (!u||find(u)) return true; 
        link[a[g].y]=u; 
    }
    g=a[g].next; 
   }
   return false; 
}
int main()
{
    scanf("%d",&ww); 
    while (ww)
    {
     add=0; 
     memset(list,0,sizeof(list)); 
     memset(link,0,sizeof(link)); 
     scanf("%d%d%d%d",&n,&m,&s,&v); ans=0;
     for (int i=1;i<=n;i++)
      scanf("%lf%lf",&w[i],&e[i]); 
     for (int i=1;i<=m;i++) 
      {
        scanf("%lf%lf",&w1,&w2); 
        for (int j=1;j<=n;j++)
        {
          double q=sqrt(cy(w[j]-w1)+cy(e[j]-w2));
          if (q<=s*v) {
            a[++add].x=j; a[add].y=i; 
            a[add].next=list[j]; list[j]=add; 
          }
        }
      }
      for (int i=1;i<=n;i++)
      {
        memset(vt,0,sizeof(vt)); 
        if (!find(i)) ans++; 
      }
      printf("%d\n",ans); 
     ww--;  
    }
}