地鼠的困境(SSL_1333)

来源:互联网 发布:南风知我意琰阙百度云 编辑:程序博客网 时间:2024/05/17 23:36

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.
原创粉丝点击