地鼠的困境SSL1333 最大匹配

题目大意

  地鼠家族一共有N个地鼠和M个鼠洞，每个都位于不同的(x, y)坐标中。假如有地鼠在发觉危险以后s秒内都没有回到鼠洞里的话，就可能成为老鹰的食物。当然了，一个鼠洞只能拯救一只地鼠的命运，所有地鼠都以相等的速度v移动。地鼠家族需要设计一种策略，使得老鹰来时，易受攻击的地鼠数量最少。

分析

  是一个最大匹配问题；

  输入时判断每一只地鼠分别可以去哪些洞，连线。

 

代码

var
a:array[1..5000,1..5000] of boolean;
v:array[1..30000] of boolean;
x,y:array[1..10000] of real;
st:array[1..30000] of longint;
i,j,k:longint;
x1,y1:real;
n,m,nm:longint;
z:longint;
v1,s:longint;

function find(r:longint):boolean;
var
i,j,k:longint;
begin
find:=true;
for i:=1 to n do
begin
if not v[i] and a[r,i]
then
begin
k:=st[i]; st[i]:=r; v[i]:=true;
if (k=0) or find(k) then exit;
st[i]:=k;
end;
end;
find:=false;
end;

procedure main;
var
i,j,k:longint;
z:boolean;
begin
for i:=1 to n do
begin
fillchar(v,sizeof(v),0);
z:=find(i);
end;
end;

begin
readln(nm);
for z:=1 to nm do
begin
fillchar(a,sizeof(a),false);
fillchar(st,sizeof(st),0);
readln(n,m,s,v1);
for i:=1 to n do
readln(x[i],y[i]);
for i:=1 to m do
begin
readln(x1,y1);
for j:=1 to n do
begin
if sqrt(sqr(x[j]-x1)+sqr(y[j]-y1))/v1<=s
then
a[j,i]:=true
end;
end;
main;
j:=0;
for i:=1 to m do
if st[i]<>0 then j:=j+1;
writeln(n-j);
end;
end.