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;
}