1333 地鼠的困境 最大匹配

题目

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

题解

二分图最大匹配，具体做法：
把每个地鼠和它能及时逃进的地洞连边，求最大匹配。
注意最终答案求有多少只危险的地鼠，所以答案是n-ans，而且有多组数据，注意数组和变量的清零（我就是提交了n遍才发现忘了清零的~)。

记得double要用%lf（%llf?)读入

代码

var
  t,n,m,s,v,i,j,k,l,ans:longint;
  q:string;
  x,y,b,c:array[1..101]of real;
  a:array[1..101,1..101]of boolean;
  cover:array[1..101]of boolean;
  link:array[1..101]of longint;

function find(i:longint):boolean;
var
  q,k:longint;
begin
  find:=true;q:=0;
  for k:=1 to m do
    if a[i,k] and cover[k] then
      begin
        q:=link[k];link[k]:=i;cover[k]:=false;
        if (q=0)or find(q) then exit;
        link[k]:=q;
      end;
  exit(false);
end;

begin
  readln(t);
  for l:=1 to t do
    begin
      fillchar(a,sizeof(a),false);
      readln(n,m,s,v);
      for i:=1 to n do
        readln(x[i],y[i]);
      ans:=0;
      for i:=1 to m do
        readln(b[i],c[i]);
      for i:=1 to n do
        for j:=1 to m do
          if sqrt(sqr(x[i]-b[j])+sqr(y[i]-c[j]))<=s*v then a[i,j]:=true;
      fillchar(link,sizeof(link),0);
      for i:=1 to n do
        begin
          fillchar(cover,sizeof(cover),true);
          if find(i) then inc(ans);
        end;
      writeln(n-ans);
    end;
end.

那么再来一个C++的？

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>

using namespace std;

int t,n,m,s,v,cnt;
int ls[105],ne[10005],y[10005],link[105];
bool cover[105];
struct node{
	double x,y;
}mice[105],hole[105];

bool dfs(int k){
	for (int i=ls[k];i;i=ne[i])
	if (!cover[y[i]]){
		int t=link[y[i]];
		link[y[i]]=k;
		cover[y[i]]=1;
		if (t==0||dfs(t)) return 1;
		link[y[i]]=t;
	}
	return 0;
}

int main(){
	scanf("%d",&t);
	for (int tt=1;tt<=t;tt++){
		scanf("%d%d%d%d",&n,&m,&s,&v);
		for (int i=1;i<=n;i++)
			scanf("%lf%lf",&mice[i].x,&mice[i].y);
		cnt=0;
		memset(link,0,sizeof(int)*(m+3));
		memset(ls,0,sizeof(int)*(n+3));
		for (int i=1;i<=m;i++){
			double xz,yz;
			scanf("%lf%lf",&xz,&yz);
			for (int j=1;j<=n;j++){
				double ccc=mice[j].x-xz,ddd=mice[j].y-yz;
				ccc=ccc*ccc;ddd=ddd*ddd;
				double dis=sqrt(ccc+ddd);
				if (dis<=s*v){
					ne[++cnt]=ls[j];ls[j]=cnt;y[cnt]=i;
				}
			}
		}
		int ans=n;
		for (int i=1;i<=n;i++){
			memset(cover,0,sizeof(bool)*(m+3));
			if (dfs(i)) ans--;
		}
		printf("%d\n",ans);
	}
}

……而你是不可奢求的梦，刹那欢喜，已经足够——改自某书