地鼠的困境 最大匹配

题意/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：

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

 

题解/solution：

        修改一下读入，然后就是最大匹配。（有多组数据，注意清零）

 

代码/Code：

type
  arr=record
    x,y:real;
  end;
var
  n,m,num,toto:longint;
  v:array [0..10001] of boolean;
  link:array [0..10001] of longint;
  a:array [0..501,0..501] of boolean;
function find(i:integer):boolean;
var
  k,q:integer;
begin
  find:=true;
  for k:=1 to n do
    if a[i,k] and not v[k] then
      begin
        q:=link[k]; link[k]:=i; v[k]:=true;
        if (q=0) or find(q) then exit;
        link[k]:=q;
    end;
  find:=false;
end;

procedure main;
var
  i:longint;
begin
  num:=0;
  fillchar(link,sizeof(link),0);
  for i:=1 to n do
    begin
      fillchar(v,sizeof(v),false);
      if find(i) then inc(num);
    end;
end;

procedure init;
var
  i,j,o,p:longint;
  xy,yx:array [0..10001] of arr;
begin
  fillchar(a,sizeof(a),0);
  readln(n,m,o,p);
  for i:=1 to n do
    readln(xy[i].x,xy[i].y);
  for i:=1 to m do
    readln(yx[i].x,yx[i].y);
  for i:=1 to n do
    for j:=1 to m do
      if sqrt(sqr(xy[i].x-yx[j].x)+sqr(xy[i].y-yx[j].y))<=o*p then
        a[i,j]:=true;
end;

begin
  readln(toto);
  while toto<>0 do
    begin
      init;
      main;
      writeln(n-num);
      dec(toto);
    end;
end.