[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\ast 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--;	
	}
}