地鼠的困境 最大匹配

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