SSLOJ 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

分析

    构建二分图：需要用到两点间距离公式 sqrt((x1-x2)^2+(y1-y2)^2)

   只要满足就连边好了

代码

   

#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
double n,m,s,v;
double a[10000][2],b[10000][2];
bool cover[10001];
int link[10001];
vector <int> f[10001];
bool find(int x)  //匹配
{
    for (int i=0;i<f[x].size();i++)
    {
        if (!cover[f[x][i]])
        {
            cover[f[x][i]]=true;
            int q=link[f[x][i]];
            link[f[x][i]]=x;
            if (q==0||find(q)) return true;
            link[f[x][i]]=q;
        }
    }
    return false;
}
int main () 
{
    int T;
    cin>>T;
    for (int ii=1;ii<=T;ii++)  
    {
        memset(link,0,sizeof(link));
        for (int i=0;i<10000;i++) //清空向量
           f[i].clear();
        cin>>n>>m>>s>>v;
        for (int i=1;i<=n;i++)
          cin>>a[i][0]>>a[i][1];
        for (int j=1;j<=m;j++)
          cin>>b[j][0]>>b[j][1];
        for (int i=1;i<=n;i++)
        {
            for (int j=1;j<=m;j++)
            {
                double x1=a[i][0],y1=a[i][1],x2=b[j][0],y2=b[j][1];
                double h=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)); //两点间距离公式
                if (h/v<=s)  //t=s/v
                  f[i].push_back(j);
            }
        }
        int ans=0;
        for (int i=1;i<=n;i++)
        {
            memset(cover,false,sizeof(cover));
            ans+=find(i);
        }
        cout<<n-ans<<endl;
    }
}