码迷,mamicode.com
首页 > 其他好文 > 详细

POJ2536 Gopher II【二分图最大匹配】

时间:2015-03-17 10:34:25      阅读:134      评论:0      收藏:0      [点我收藏+]

标签:

题目链接:

http://poj.org/problem?id=2536


题目大意:

有N只鼹鼠和M个洞穴,如果鼹鼠在S秒内不能够跑到洞穴,就会被老鹰捉住吃掉。鼹鼠跑的速度为V米/秒。

已知一个洞穴只能容纳一只鼹鼠。给你鼹鼠和洞穴的坐标,那么问题来了:问最少有多少只鼹鼠被老鹰捉住

吃掉。


思路:

建立一个二分图,一边为鼹鼠,另一边为洞穴枚举求出每只鼹鼠到各个洞穴的距离,把能够在S秒内跑到该

穴(距离<=S*V)的进行连边。建好图后用匈牙利算法求出最多有多少只鼹鼠能够幸免于难( MaxMatch() ),

那么剩下的N - MaxMatch()就是最少有多少只鼹鼠被老鹰捉住吃掉。


AC代码:

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;

int N,M,S,V;

struct Node
{
    double x;
    double y;
}PA[330],PB[330];

int Map[330][330];
bool Mask[330];
int cx[330],cy[330];

int FindPath(int u)
{
    for(int i = 1; i <= M; ++i)
    {
        if(Map[u][i] && !Mask[i])
        {
            Mask[i] = 1;
            if(cy[i] == -1 || FindPath(cy[i]))
            {
                cy[i] = u;
                cx[u] = i;
                return 1;
            }
        }
    }
    return 0;
}

int MaxMatch()
{
    int res = 0;
    for(int i = 1; i <= N; ++i)
        cx[i] = -1;
    for(int i = 1; i <= M; ++i)
        cy[i] = -1;

    for(int i = 1; i <= N; ++i)
    {
        if(cx[i] == -1)
        {
            for(int j = 1; j <= M; ++j)
                Mask[j] = 0;
            res += FindPath(i);
        }
    }
    return res;
}

int main()
{
    while(~scanf("%d%d%d%d",&N,&M,&S,&V))
    {
        for(int i = 1; i <= N; ++i)
            scanf("%lf%lf",&PA[i].x,&PA[i].y);
        for(int i = 1; i <= M; ++i)
            scanf("%lf%lf",&PB[i].x,&PB[i].y);

        memset(Map,0,sizeof(Map));
        for(int i = 1; i <= N; ++i)
        {
            for(int j = 1; j <= M; ++j)
            {
                double x = PA[i].x - PB[j].x;
                double y = PA[i].y - PB[j].y;
                if(x*x+y*y <= S*V*S*V)
                    Map[i][j] = 1;
            }
        }
        printf("%d\n",N-MaxMatch());
    }

    return 0;
}


POJ2536 Gopher II【二分图最大匹配】

标签:

原文地址:http://blog.csdn.net/lianai911/article/details/44337225

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!