Gym - 101666E -- Easter Eggs

版权声明: https://blog.csdn.net/moon_sky1999/article/details/81545392

题目来源:http://codeforces.com/gym/101666/attachments

由于是求最小距离的最大值,因此很容易想到二分答案。

在验证每一个值的过程中,我们可以将距离<mid的边保留,>mid的边删除,因此变成了求二分图的最大独立集的问题。

相关概念及定理:https://blog.csdn.net/moon_sky1999/article/details/81331795

由定理知:|最大独立集|+|最小顶点覆盖|=|V|

对于二分图:|最大匹配|=|最小顶点覆盖|。

因此可以对处理后的图进行二分图匹配,求得最大匹配。即可验证。

代码:

#include <bits/stdc++.h>
using namespace std;
const double eps = 1e-9;
int n, b, r, from[510];
double xb[510], yb[510], xr[510], yr[510], g[510][510], e[510][510];
bool vis[510];

bool match(int x) {
    for (int i = 1; i <= n; ++i) {
        if (e[x][i] > eps && vis[i] == 0) {
            vis[i] = 1;
            if (from[i] == -1 || match(from[i])) {
                from[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int hungry() {
    int tot = 0;
    memset(from, -1, sizeof(from));
    for (int i = 1; i <= b; ++i) {
        memset(vis, 0, sizeof(vis));
        if (match(i))
            ++tot;
    }
    return tot;
}

bool can(double x) {
    memset(e, 0, sizeof(e));
    for (int i = 1; i <= b; ++i) {
        for (int j = 1; j <= r; ++j)
            if (g[i][j] < x)
                e[i][j] = g[i][j];
    }
    int tot = hungry();
    if (b + r - tot >= n)return 1;
    return 0;
}

int main() {
    scanf("%d%d%d", &n, &b, &r);
    for (int i = 1; i <= b; ++i)
        scanf("%lf%lf", &xb[i], &yb[i]);
    for (int i = 1; i <= r; ++i)
        scanf("%lf%lf", &xr[i], &yr[i]);
    for (int i = 1; i <= b; ++i)
        for (int j = 1; j <= r; ++j)
            g[i][j] = sqrt((xb[i] - xr[j]) * (xb[i] - xr[j]) + (yb[i] - yr[j]) * (yb[i] - yr[j]));
    double le = 0.0, ri = 10000.0;
    while (ri - le > eps) {
        double mid = (le + ri) / 2;
        if (can(mid))le = mid;
        else ri = mid;
    }
    printf("%.9f\n", ri);
    return 0;
}

猜你喜欢

转载自blog.csdn.net/moon_sky1999/article/details/81545392