题目
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
题解
二分图最大匹配,具体做法:
把每个地鼠和它能及时逃进的地洞连边,求最大匹配。
注意最终答案求有多少只危险的地鼠,所以答案是n-ans,而且有多组数据,注意数组和变量的清零(我就是提交了n遍才发现忘了清零的~)。
记得double要用%lf(%llf?)读入
代码
var
t,n,m,s,v,i,j,k,l,ans:longint;
q:string;
x,y,b,c:array[1..101]of real;
a:array[1..101,1..101]of boolean;
cover:array[1..101]of boolean;
link:array[1..101]of longint;
function find(i:longint):boolean;
var
q,k:longint;
begin
find:=true;q:=0;
for k:=1 to m do
if a[i,k] and cover[k] then
begin
q:=link[k];link[k]:=i;cover[k]:=false;
if (q=0)or find(q) then exit;
link[k]:=q;
end;
exit(false);
end;
begin
readln(t);
for l:=1 to t do
begin
fillchar(a,sizeof(a),false);
readln(n,m,s,v);
for i:=1 to n do
readln(x[i],y[i]);
ans:=0;
for i:=1 to m do
readln(b[i],c[i]);
for i:=1 to n do
for j:=1 to m do
if sqrt(sqr(x[i]-b[j])+sqr(y[i]-c[j]))<=s*v then a[i,j]:=true;
fillchar(link,sizeof(link),0);
for i:=1 to n do
begin
fillchar(cover,sizeof(cover),true);
if find(i) then inc(ans);
end;
writeln(n-ans);
end;
end.
那么再来一个C++的?
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
int t,n,m,s,v,cnt;
int ls[105],ne[10005],y[10005],link[105];
bool cover[105];
struct node{
double x,y;
}mice[105],hole[105];
bool dfs(int k){
for (int i=ls[k];i;i=ne[i])
if (!cover[y[i]]){
int t=link[y[i]];
link[y[i]]=k;
cover[y[i]]=1;
if (t==0||dfs(t)) return 1;
link[y[i]]=t;
}
return 0;
}
int main(){
scanf("%d",&t);
for (int tt=1;tt<=t;tt++){
scanf("%d%d%d%d",&n,&m,&s,&v);
for (int i=1;i<=n;i++)
scanf("%lf%lf",&mice[i].x,&mice[i].y);
cnt=0;
memset(link,0,sizeof(int)*(m+3));
memset(ls,0,sizeof(int)*(n+3));
for (int i=1;i<=m;i++){
double xz,yz;
scanf("%lf%lf",&xz,&yz);
for (int j=1;j<=n;j++){
double ccc=mice[j].x-xz,ddd=mice[j].y-yz;
ccc=ccc*ccc;ddd=ddd*ddd;
double dis=sqrt(ccc+ddd);
if (dis<=s*v){
ne[++cnt]=ls[j];ls[j]=cnt;y[cnt]=i;
}
}
}
int ans=n;
for (int i=1;i<=n;i++){
memset(cover,0,sizeof(bool)*(m+3));
if (dfs(i)) ans--;
}
printf("%d\n",ans);
}
}