SSL1333 地鼠的困境（最大匹配）

地鼠的困境

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

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

分析：建个二分图然后最大匹配。

代码

#include <cstdio>
#include <cmath>
using namespace std;

bool f[200][200],cover[200];
double a[200][2],b[200][2];
int T,n,m,s,v,ans,link[200];

void init()
{
scanf("%d%d%d%d",&n,&m,&s,&v);
for (int i=1;i<=n;i++)
scanf("%lf%lf",&a[i][1],&a[i][2]);
for (int i=1;i<=m;i++)
scanf("%lf%lf",&b[i][1],&b[i][2]);
int s3=v*s;
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
{
double s1=a[i][1]-b[j][1],s2=a[i][2]-b[j][2];
if (sqrt(s1*s1+s2*s2)<=s3) f[i][j]=true;
}
}

bool find(int i)
{
for (int k=1;k<=m;k++)
if (f[i][k]&&(!cover[k]))
{
int q=link[k];
link[k]=i;
cover[k]=true;
if ((q==0)||(find(q))) return true;
link[k]=q;
}
return false;
}

void work()
{
for (int i=1;i<=n;i++)
{
for (int j=1;j<=100;j++)
cover[j]=false;
if (find (i)) continue;
ans++;
}
}

int main()
{
scanf("%d",&T);
for (int ii=1;ii<=T;ii++)
{
init();
work();
printf("%d\n",ans);
ans=0;
for (int i=1;i<=100;i++)
{
link[i]=0;
for (int j=1;j<=100;j++)
f[i][j]=false;
}
}
return 0;
}