地鼠的困境_ssl1333_匹配



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

思路：

枚举两两的洞和鼹鼠判断，可以在规定时间回到洞里的就连边，二分图求最大匹配。

注意要用总数减去最大匹配，这才是题目要求的解

源代码/pas：

 

type
  edge=record
    x,y,next:longint;
  end;
  point=record
    x,y:real;
  end;
var
  v:array[1..100]of boolean;
  link,ls:array[1..100]of longint;
  e:array[1..1000]of edge;
  maxE,n,m,a,b,g:longint;
procedure add(x,y:longint);
begin
  inc(maxE);
  e[maxE].x:=x;
  e[maxE].y:=y;
  e[maxE].next:=ls[x];
  ls[x]:=maxE;
end;
function find(x:longint):boolean;
var
  i,k:longint;
begin
  find:=true;
  i:=ls[x];
  while i>0 do
  with e[i] do
  begin
    if not v[y] then
    begin
      v[y]:=true;
      k:=link[y];
      link[y]:=x;
      if (k=0)or find(k) then exit;
      link[y]:=k;
    end;
    i:=next;
  end;
  find:=false;
end;
procedure main;
var
  i,ans:longint;
begin
  ans:=0;
  fillchar(link,sizeof(link),0);
  for i:=1 to n do
  begin
    fillchar(v,sizeof(v),false);
    if find(i) then inc(ans);
  end;
  writeln(n-ans);
end;
procedure init;
var
  i,j:longint;
  x,y:array[1..100]of point;
begin
  maxE:=0;
  fillchar(x,sizeof(x),0);
  fillchar(y,sizeof(y),0);
  fillchar(ls,sizeof(ls),0);
  fillchar(e,sizeof(e),0);
  readln(n,m,a,b);
  for i:=1 to n do readln(x[i].x,x[i].y);
  for i:=1 to m do readln(y[i].x,y[i].y);
  for i:=1 to n do
  for j:=1 to m do
  if sqrt(sqr(x[i].x-y[j].x)+sqr(x[i].y-y[j].y))/a<=b then
  add(i,j);
end;
begin
  readln(g);
  repeat
    init;
    main;
    dec(g);
  until g=0;
end.