地鼠的困境(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.