hdu 2389(最大匹配bfs版）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2389

思路：纯裸的一个最大匹配题，不过悲摧的是以前一直用的dfs版一直过不了,TLE无数次啊，然后改成bfs就过了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
#include<queue>
using namespace std;
#define MAXN 3333
vector<int>map[MAXN];
int ly[MAXN],lx[MAXN];
int distx[MAXN],disty[MAXN];
int n,m;
struct Node {
    int x,y,s;
} p[MAXN];
struct Point {
    int x,y;
} q[MAXN];

bool bfs()
{
    bool flag=false;
    memset(distx,0,(n+2)*sizeof(int));
    memset(disty,0,(n+2)*sizeof(int));
    queue<int>Q;
    for(int i=1; i<=n; i++) {
        if(lx[i]==-1)Q.push(i);
    }
    while(!Q.empty()) {
        int u=Q.front();
        Q.pop();
        for(int i=0; i<map[u].size(); i++) {
            int v=map[u][i];
            if(disty[v]==0) {
                disty[v]=distx[u]+1;
                if(ly[v]==-1)flag=true;
                else {
                    distx[ly[v]]=disty[v]+1;
                    Q.push(ly[v]);
                }
            }
        }
    }
    return flag;
}

int dfs(int u)
{
    for(int i=0; i<map[u].size(); i++) {
        int v=map[u][i];
        if(disty[v]==distx[u]+1) {
            disty[v]=0;
            if(ly[v]==-1||dfs(ly[v])) {
                ly[v]=u;
                lx[u]=v;
                return 1;
            }
        }
    }
    return 0;
}


int MaxMatch()
{
    int res=0;
    memset(lx,-1,(n+2)*sizeof(lx[0]));
    memset(ly,-1,(n+2)*sizeof(ly[0]));
    while(bfs()) {
        for(int i=1; i<=n; i++) {
            if(lx[i]==-1)res+=dfs(i);
        }
    }
    return res;
}

int main()
{
    //  freopen("1.txt","r",stdin);
    int _case,time,t=1;
    scanf("%d",&_case);
    while(_case--) {
        scanf("%d%d",&time,&n);
        for(int i=1; i<=n; i++)map[i].clear();
        for(int i=1; i<=n; i++) {
            scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].s);
        }
        scanf("%d",&m);
        for(int i=1; i<=m; i++) {
            scanf("%d%d",&q[i].x,&q[i].y);
        }
        for(int i=1; i<=n; i++) {
            for(int j=1; j<=m; j++) {
                double xx=(p[i].x-q[j].x)*(p[i].x-q[j].x)*1.0;
                double yy=(p[i].y-q[j].y)*(p[i].y-q[j].y)*1.0;
                double dd=sqrt(xx+yy);
                if(dd<=1.0*time*p[i].s)map[i].push_back(j);
            }
        }
        int ans=MaxMatch();
        printf("Scenario #%d:
%d
",t++,ans);
        puts("");
    }
    return 0;
}