poj 2536 Gopher II 最大匹配

该博客主要介绍了POJ 2536 Gopher II问题,通过建立二分图并计算最大匹配来求解如何最大化存活人数。作者提供了题目的来源链接和详细的题解思路,以及利用图的最大匹配规模确定最少死亡人数的算法。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目

题目链接:http://poj.org/problem?id=2536

题目来源:http://www.cnblogs.com/vongang/archive/2012/02/21/2361882.html

题解

最小化死掉的人,也就是最大化存活的人。

人和洞建二分图,根据题中信息算出某人对于某个洞是否能够到达,能到达就连上边。

图上的最大匹配的规模就是最多的存活的人, n|M| 就是答案。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <stack>
#include <queue>
#include <string>
#include <vector>
#include <set>
#include <map>
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
typedef pair<double, double> pdd;
// head
const int N = 105;
const int M = N;
const int E = N * M;

struct Gragh {
    struct Edge {
        int to, nxt;
        Edge(int to, int nxt) : to(to), nxt(nxt) {}
        Edge() {}
    };

    Edge e[E];
    int head[N], match[M], ec;
    bool vis[N];

    void addEdge(int from, int to) {
        e[ec] = Edge(to, head[from]);
        head[from] = ec++;
    }
    void init(int n, int m) {
        ec = 0;
        for (int i = 0; i <= n; i++) head[i] = -1;
        for (int i = 0; i <= m; i++) match[i] = -1;
    }
    bool dfs(int x) {
        vis[x] = true;
        for (int i = head[x]; ~i; i = e[i].nxt) {
            int to = e[i].to, w = match[to];
            if (w == -1 || (!vis[w] && dfs(w))) {
                match[to] = x;
                return true;
            }
        }
        return false;
    }
    int bipatiteMatch(int n) {
        int ans = 0;
        for (int i = 0; i < n; i++) {
            memset(vis, 0, sizeof vis);
            if (dfs(i)) ans++;
        }
        return ans;
    }
};

inline double sq(double x) {
    return x * x;
}
double dis(pdd &a, pdd &b) {
    return sqrt(sq(a.fi - b.fi) + sq(a.se - b.se));
}

Gragh g;
pdd a[N], b[N];

int main() {
    int n, m, s, v;
    while (scanf("%d%d%d%d", &n, &m, &s, &v) == 4) {
        s *= v;
        g.init(n, m);
        for (int i = 0; i < n; i++) {
            scanf("%lf%lf", &a[i].fi, &a[i].se);
        }
        for (int i = 0; i < m; i++) {
            scanf("%lf%lf", &b[i].fi, &b[i].se);
        }
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < m; j++) {
                if (dis(a[i], b[j]) > s) continue;
                g.addEdge(i, j);
            }
        }
        printf("%d\n", n - g.bipatiteMatch(n));
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值