題目連結：點我啊╭(╯^╰)╮

題目大意：

     $n$ 個城市、 $m$ 個雷達，最多可以用 $k$

解題思路：

    很清晰是一道DLX重複覆蓋的題，本題的精度要精確到 $1e-8$

程式碼思路：

    主要還是重複覆蓋的模板題，注意是用了 A* 的優化，估價函式 $f$ 意義為從當前狀態最好情況下需要新增幾條邊才能覆蓋

核心：二分 + DLX 的簡單題，同時注意精度

$PS：$本題的 dfs 返回值可以為 $bool$，因為可以用當前深度與 $k$ 進行比較

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

#define N 3020
#define M 55
struct DancingLink{
    int n,s,ansd;//列數 節點總數 
    int S[M],A[M],H[M];//S[]該列節點總數  A[]答案  H[]行首指標 
    int L[N],R[N],U[N],D[N]; //L[],R[],U[],D[] 上下左右 
    int X[N],C[N],vis[M];//X[] C[] 行列編號 
    void init(int n){//初始化 
        this->n=n;
        for(int i=0;i<=n;i++)
            U[i]=i,D[i]=i,L[i]=i-1,R[i]=i+1;
        R[n]=0,L[0]=n;s=n+1;
		memset(S,0,sizeof(S));
		memset(H,-1,sizeof(H));
    }
    void DelCol(int c){//刪除列 
        for(int i=D[c];i!=c;i=D[i])
        	L[R[i]]=L[i],R[L[i]]=R[i];
    }
    void ResCol(int c){//恢復列 
        for(int i=U[c];i!=c;i=U[i])
        	L[R[i]]=R[L[i]]=i;
    }
    void AddNode(int r,int c){//新增節點 
        ++S[c],C[++s]=c,X[s]=r;
        D[s]=D[c],U[D[c]]=s,U[s]=c,D[c]=s;
        if(H[r]<0) H[r]=L[s]=R[s]=s;//行首節點
        else  R[s]=R[H[r]],L[R[H[r]]]=s,L[s]=H[r],R[H[r]]=s;
    }
    int f(){
    	int ret=0;
        for(int i=R[0];i;i=R[i]) vis[i]=0;
        for(int i=R[0];i;i=R[i])
			if(!vis[i]){
                ret++;vis[i]=1;
                for(int j=D[i];j!=i;j=D[j])
                	for(int k=R[j];k!=j;k=R[k])
                		vis[C[k]]=1;	
            }
        return ret;
	}
    void dfs(int d){//深度，深搜遍歷 
    	if(d+f()>=ansd) return;
        if(!R[0]){
            ansd=min(d,ansd);
			return;
        }
        int c=R[0];
        for(int i=R[0];i;i=R[i]) if(S[i]<S[c]) c=i;
        for(int i=D[c];i!=c;i=D[i]){
        	DelCol(i);
            A[d]=X[i];
            for(int j=R[i];j!=i;j=R[j]) DelCol(j);
            dfs(d+1);
            for(int j=L[i];j!=i;j=L[j]) ResCol(j);
            ResCol(i);
        }
    }

} dlx; 

double eps = 1e-8;
struct point{
	double x, y;
} ci[60], sta[60];

double dist(point a, point b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int main(){
	int n, m, k, cas;
	scanf("%d", &cas);
	while(cas--) {
		scanf("%d%d%d", &n, &m, &k);
		for(int i=0; i<n; i++) 	scanf("%lf%lf", &ci[i].x, &ci[i].y);
		for(int i=0; i<m; i++) 	scanf("%lf%lf", &sta[i].x, &sta[i].y);
		double l=0, r=2000, mid;
		while(r-l>=eps){
			dlx.init(n);
			mid = (l+r)/2;
			for(int i=0; i<m; i++)
				for(int j=0; j<n; j++)
					if(dist(ci[j],sta[i])<=mid)
						dlx.AddNode(i+1,j+1);
			dlx.ansd = 1<<30;
			dlx.dfs(0);
			if(dlx.ansd<=k) r = mid-eps;
			else l = mid+eps;
		}
		printf("%.6f\n", l);
	}
}