[最大独立集]地鼠的困境

题目描述

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

输入输出格式

输入格式

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

输出格式

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

输入输出样例

输入样例#1：

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

输出样例#1：

1

题目解析

枚举两两的洞和鼠判断，可以在规定时间回到洞里的就连边，二分图求最大匹配。
注意要用总数减去最大匹配，这才是题目要求的解

代码

#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
#include<cmath>
using namespace std;
int t,n,m,s,v,link[105],ans;
bool flag[105];
vector <int> a[105];
double mus[105][4],hole[105][4];
bool find(int x)
{
	for(int i=0;i<a[x].size();i++)
	 if(!flag[a[x][i]])
	 {
	   int j=a[x][i];
	   flag[j]=1;
	   int q=link[j];
	   link[j]=x;
	   if(!q||find(q)) return true;
	   link[j]=q;
	 }
	return false;
}//最大匹配
int main()
{
	cin>>t;
	while(t>0)
	{
	  memset(link,0,sizeof(link));
	  for(int i=1;i<=n;i++)
	   a[i].clear();
	  ans=0;
	  t--;
	  cin>>n>>m>>s>>v;
	  for(int i=1;i<=n;i++)
	   cin>>mus[i][1]>>mus[i][2];
	  for(int i=1;i<=m;i++)
	   cin>>hole[i][1]>>hole[i][2];
	  for(int i=1;i<=m;i++)
	   for(int j=1;j<=n;j++)
	    if((mus[j][1]-hole[i][1])*(mus[j][1]-hole[i][1])+(mus[j][2]-hole[i][2])*(mus[j][2]-hole[i][2])<s*v*s*v)
	   //判断第j只鼠是否可以到第i个洞
		 a[j].push_back(i);//连边
	  for(int i=1;i<=n;i++)
	  {
	  	memset(flag,0,sizeof(flag));
	  	ans+=find(i);
	  }
	  cout<<n-ans<<endl;题目的解为总数减去最大匹配
	}
}