地鼠的困境(SSL_1333)

Time Limit:1000MS Memory Limit:65536K

Total Submit:225 Accepted:59

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

思路：先求出在S秒后地鼠可以跑到的最大路长（w），再枚举每只地鼠，判断可以到哪个洞，如果可以就在map里面放true，最后套匈牙利算法

var
cover:array[0..1001] of boolean;
link:array[0..1001] of longint;
map:array[0..1001,0..1001] of boolean;
i,k,n,s,v,t,m,j,nn,w,ans:longint;
x1,y1,x,y:array[0..1001]of real;
ss,x2,y2,q,p:real;
function find(i:longint):boolean;
var
k,q,s:longint;
begin
find:=true;
for k:=1 to n do
if (map[i,k])and(not(cover[k])) then
begin
q:=link[k];
cover[k]:=true;
link[k]:=i;
if (q=0)or(not(find(q))) then exit;
link[k]:=q;
end;
exit(false);
end;

begin
readln(nn);
for k:=1 to nn do
begin
readln(n,m,s,v);
w:=s*v;
ans:=0;
fillchar(link,sizeof(link),0);
fillchar(cover,sizeof(cover),false);
fillchar(map,sizeof(map),false);
fillchar(x,sizeof(x),0);
fillchar(y,sizeof(y),0);
for i:=1 to n do
readln(x[i],y[i]);
for i:=1 to m do
begin
readln(q,p);
for j:=1 to n do
begin
x2:=q-x[j];
y2:=p-y[j];
ss:=sqrt(sqr(x2)+sqr(y2));
if w>=ss then map[j,i]:=true;
end;
end;
for j:=1 to n do
begin
fillchar(cover,sizeof(cover),false);
find(j);
end;
for j:=1 to n do
if link[j]>0 then inc(ans);
writeln(n-ans);
end;
end.