GYM 101147 C.The Wall（二分匹配-hungary）

Description
一个n*m的区域，有n+m个防御工事，x=i(0<=i<=n-1)可以杀死i<=x < i+1里的所有敌人，以(n/2,y)(0<=y<=m-1)为圆心,n/2为半径的上半圆弧会对杀死前面距该圆弧距离不超过1的敌人(注意不能跨圆弧杀人)，现在给出p个敌人的坐标，问至少需要开放多少防御工事可以消灭多少敌人
Input
第一行一整数T表示用例组数，每组用例首先输入三个整数n,m,p分别表示区域的范围和敌人人数，之后p行每行两个实数x,y表示该敌人的位置(1<=n,m<=100, 1<=p<=1e4, 0 < x < n, 0 < y < m+n/2)
Output
输出消灭所有敌人最少需要开放多少防御工事
Sample Input
1
4 4 9
3.1 1.9
0.4 1.4
1.4 2.6
3.2 2.2
2.2 5.6
2.7 4.3
1.2 4.0
1.2 4.5
2.8 3.6
Sample Output
3
Solution
把n个竖直的防御工事看作一排点，把m个弧形的防御工事看作一排点，一个敌人必然会被一个竖直的防御工事x或者一个弧形的防御工事y消灭，那么对于每个敌人对应的x和y连边，最后得到一张二分图，边表示一名敌人，要的是用最少的点来覆盖所有的边，即求该二分图的最小点覆盖，也即为二分图最大匹配
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 111
#define eps 1e-8
int uN,vN;
int g[maxn][maxn];
int linker[maxn];
bool used[maxn];
bool dfs(int u)
{
    for(int v=0;v<vN;v++)
        if(g[u][v]&&!used[v])
        {
            used[v]=1;
            if(linker[v]==-1||dfs(linker[v]))
            {
                linker[v]=u;
                return 1;
            }    
        }  
    return 0;  
}    
int hungary()
{
    int res=0;
    memset(linker,-1,sizeof(linker));
    for(int u=0;u<uN;u++)
    {
        memset(used,0,sizeof(used));
        if(dfs(u))res++;
    } 
    return res;   
}  
int sign(double x)
{
    if(fabs(x)<eps)return 0;
    if(x>eps)return 1;
    return -1;
}
int T,n,m,p;
bool check(double x,double y,double z)
{
    double r=0.5*n;
    if(sign(y-z)<0)return 0;
    double dis=sqrt((r-x)*(r-x)+(z-y)*(z-y));
    if(sign(dis-r)>0&&sign(dis-r-1)<0)return 1;
    return 0; 
}
int main()
{
    freopen("wall.in","r",stdin);
    scanf("%d",&T);
    while(T--)
    {
        memset(g,0,sizeof(g));
        scanf("%d%d%d",&n,&m,&p);
        uN=n,vN=m;
        for(int i=1;i<=p;i++)
        {
            double x,y;
            int a=-1,b=-1;
            scanf("%lf%lf",&x,&y);
            a=(int)x;
            for(int j=0;j<m;j++)
                if(check(x,y,1.0*j))b=j;
            if(a!=-1&&b!=-1)g[a][b]=1;
        }
        int ans=hungary();
        printf("%d\n",ans);
    }
    return 0;
}