hdu2389（HK算法）

传送门：Rain on your Parade

题意：t个单位时间后开始下雨，给你N个访客的位置（一维坐标平面内）和他的移动速度，再给M个雨伞的位置，问在下雨前最多有多少人能够拿到雨伞（两个人不共用一把伞）。

分析：这题匈牙利算法撸不过，只好去学习Hopcroft-Carp算法，复杂度为O(sqrt(V)*E)，该算法预先找好增广路径集，避免了许多没不要的匹配.

#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstdlib>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define LL long long
#define mod 100000000
#define inf 0x3f3f3f3f
#define eps 1e-6
#define N 10010
#define FILL(a,b) (memset(a,b,sizeof(a)))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define PII pair<int,int>
using namespace std;
int dx[N],dy[N],matchx[N],matchy[N];
int vis[N],dis,n,m,t;
vector<int>g[N];
struct node
{
    int x,y,z;
}s[N];
bool search_path()
{
    queue<int>que;
    dis=inf;
    FILL(dx,-1);FILL(dy,-1);
    for(int i=1;i<=n;i++)
    {
        if(matchx[i]==-1)
        {
            que.push(i);
            dx[i]=0;
        }
    }
    while(!que.empty())
    {
        int u=que.front();
        que.pop();
        if(dx[u]>dis)break;
        for(int i=0,sz=g[u].size();i<sz;i++)
        {
            int v=g[u][i];
            if(dy[v]==-1)
            {
                dy[v]=dx[u]+1;
                if(matchy[v]==-1)dis=dy[v];
                else
                {
                    dx[matchy[v]]=dy[v]+1;
                    que.push(matchy[v]);
                }
            }
        }
    }
    return dis!=inf;
}
int dfs(int u)
{
    for(int i=0,sz=g[u].size();i<sz;i++)
    {
        int v=g[u][i];
        if(!vis[v]&&dy[v]==dx[u]+1)
        {
            vis[v]=1;
            if(matchy[v]!=-1&&dy[v]==dis)continue;
            if(matchy[v]==-1||dfs(matchy[v]))
            {
                matchy[v]=u;
                matchx[u]=v;
                return 1;
            }
        }
    }
    return 0;
}
int HK()
{
    int res=0;
    FILL(matchx,-1);FILL(matchy,-1);
    while(search_path())
    {
        FILL(vis,0);
        for(int i=1;i<=n;i++)
            if(matchx[i]==-1&&dfs(i))res++;
    }
    return res;
}
bool judge(int a,int b,int x,int y,int z)
{
    return (a-x)*(a-x)+(b-y)*(b-y)<=(z*t)*(z*t);
}
int main()
{
    int T,cas=1;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&t,&n);
        for(int i=1;i<=n;i++)g[i].clear();
        for(int i=1;i<=n;i++)
            scanf("%d%d%d",&s[i].x,&s[i].y,&s[i].z);
        scanf("%d",&m);
        for(int i=1;i<=m;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            for(int j=1;j<=n;j++)
            {
                if(judge(x,y,s[j].x,s[j].y,s[j].z))
                    g[j].push_back(i);
            }
        }
        int ans=HK();
        printf("Scenario #%d:
%d
",cas++,ans);
        puts("");
    }
}