【HDU】2295 Radar

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXM 110
#define MAXN 50000
#define EPS 1e-8
#define INF 0x7FFFFFFF
using namespace std;
int n, m, k, size;
int L[MAXN], R[MAXN], U[MAXN], D[MAXN], C[MAXN];
int H[MAXN], S[MAXN];
bool vis[MAXM];
struct Point {
    int x, y;
};
Point city[MAXM], radar[MAXM];
inline int dbcmp(double x, double y) {
    if (fabs(x - y) < EPS)
        return 0;
    return x > y ? 1 : -1;
}
double DIS(int i, int j) {
    double x, y;
    x = city[i].x - radar[j].x;
    y = city[i].y - radar[j].y;
    return x * x + y * y;
}
void Init() {
    int i;
    for (i = 0; i <= m; i++)
        H[i] = -1;
    for (i = 0; i <= n; i++) {
        R[i] = i + 1;
        L[i + 1] = i;
        U[i] = D[i] = i;
        S[i] = 0;
    }
    R[n] = 0;
    size = n + 1;
}
void Link(int r, int c) {
    D[size] = D[c];
    U[size] = c;
    U[D[c]] = size;
    D[c] = size;
    if (H[r] < 0)
        H[r] = L[size] = R[size] = size;
    else {
        L[size] = H[r];
        R[size] = R[H[r]];
        L[R[H[r]]] = size;
        R[H[r]] = size;
    }
    S[c]++;
    C[size++] = c;
}
int A() {
    int i, j, k, res;
    memset(vis, false, sizeof(vis));
    for (res = 0, i = R[0]; i; i = R[i]) {
        if (vis[i])
            continue;
        res++;
        for (j = D[i]; j != i; j = D[j]) {
            for (k = R[j]; k != j; k = R[k])
                vis[C[k]] = true;
        }
    }
    return res;
}
void Remove(int c) {
    int i;
    for (i = D[c]; i != c; i = D[i]) {
        L[R[i]] = L[i];
        R[L[i]] = R[i];
    }
}
void Resume(int c) {
    int i;
    for (i = D[c]; i != c; i = D[i])
        L[R[i]] = R[L[i]] = i;
}
bool Dance(int now) {
    if (R[0] == 0)
        return true;
    else if (now + A() <= k) {
        int i, j, temp, c;
        for (temp = INF,i = R[0]; i; i = R[i]) {
            if (temp > S[i]) {
                temp = S[i];
                c = i;
            }
        }
        for (i = D[c]; i != c; i = D[i]) {
            Remove(i);
            for (j = R[i]; j != i; j = R[j])
                Remove(j);
            if (Dance(now + 1))
                return true;
            for (j = L[i]; j != i; j = L[j])
                Resume(j);
            Resume(i);
        }
    }
    return false;
}
bool OK(double mid) {
    int i, j;
    Init();
    mid *= mid;
    for (i = 1; i <= n; i++) {
        for (j = 1; j <= m; j++) {
            if (dbcmp(DIS(i, j), mid) <= 0)
                Link(j, i);
        }
    }
    return Dance(0);
}
int main() {
    int T;
    int i;
    double low, high, mid;
    scanf("%d", &T);
    while (T--) {
        scanf("%d%d%d", &n, &m, &k);
        for (i = 1; i <= n; i++)
            scanf("%d%d", &city[i].x, &city[i].y);
        for (i = 1; i <= m; i++)
            scanf("%d%d", &radar[i].x, &radar[i].y);
        for (low = 0, high = 1e5; dbcmp(low, high) < 0;) {
            mid = (low + high) / 2;
            if (OK(mid))
                high = mid;
            else
                low = mid;
        }
        printf("%.6lf\n", high);
    }
    return 0;
}