MYF

POJ 2536 Gopher II

题目链接

POJ 2536

方法:二分图

题目分析

题目大意

给出n只老鼠的坐标和m个洞的坐标,每个洞只能躲一只老鼠,还告诉了老鼠奔跑的速度以及老鹰到达的时间,问最少有几只老鼠会找不到洞被老鹰吃掉。

解析

二分图的题目。

n个老鼠直接跑到m个洞里,每个老鼠都有m种选择,当然还要除去那些根本来不及跑到老鹰就来了的洞穴。如果尽可能的保证躲在洞里的老鼠足够多,那么就相当于左子集为老鼠的下标,右子集为洞洞下标,然后跑一次最大匹配。根据距离来加边,然后跑二分图即可。

于是无脑的用二分图的模版即可。

代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
#include <set>
#include <map>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <string>
#include <vector>
#include <iomanip>
#include <bitset>
#include <cstring>
#include <iostream>
#include <deque>
#include <algorithm>
#define Memset(a,val) memset(a,val,sizeof(a))
#define PI acos(-1)
#define PB push_back
#define MP make_pair
#define rt(n) (i == n ? '\n' : ' ')
#define hi printf("Hi----------\n")
#define IN freopen("input.txt","r",stdin);
#define OUT freopen("output.txt","w",stdout);
#define debug(x) cout<<"Debug : ---"<<x<<"---"<<endl;
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
typedef pair<int,int> PII;
typedef long long ll;
const int mod=1000000007;
const int INF=0x3f3f3f3f;
const double eps=1e-8;
const int MAXN=1000;
int uN,vN; //u,v数目
double g[MAXN][MAXN];//编号是0~n-1的
int linker[MAXN];
bool used[MAXN];
struct Point {
double x, y;
void in(){
scanf("%lf%lf",&x,&y);
}
}rat[MAXN],hole[MAXN];
double dis(Point x,Point y){
return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));

}
bool dfs(int u)
{
int v;
for(v=1;v<=vN;v++)
if(g[u][v]&&!used[v])
{
used[v]=true;
if(linker[v]==-1||dfs(linker[v]))
{
linker[v]=u;
return true;
}
}
return false;
}
int hungary()
{
int res=0;
int u;
memset(linker,-1,sizeof(linker));
for(u=1;u<=uN;u++)
{
memset(used,0,sizeof(used));
if(dfs(u)) res++;
}
return res;
}
int main(){
int n,m,s,v;
while (scanf("%d%d%d%d",&n,&m,&s,&v)!=EOF) {
Memset(g, 0);
uN=n,vN=m;
for (int i=1; i<=n; i++) {
rat[i].in();
}
for (int i=1; i<=m; i++) {
hole[i].in();
}
for (int i=1; i<=n; i++) {
for (int j=1; j<=m; j++) {
if (dis(rat[i],hole[j])<s*v) {
g[i][j]=1;
}
}
}
cout<<n-hungary()<<endl;
}
return 0;
}