题意:有n个点。选取k个点(k> = 2),并在这k点中最近的两个点距离最远。
题解:二分+最大团
最近的最远,往二分上面靠。
由于最近的两个点所连的边是k点构成完全图中边权最小的,那么其他的边必然大于等于该值,那么我们二分这个答案即可。
对于二分的mid,每次将大于等于mid的边连上,然后跑最大团(点点直接相连才能保证最短边权,若不是完全图则必然存在小于mid的边),若最大团的结点数大于等于k,则最大团存在结点数等于k,且最小边权大于等于mid的子图。
#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<set>
#include<map>
#include<sstream>
#include<iomanip>
#define ll long long
#define pii pair<int, int>
using namespace std;
struct node {
double x, y;
}p[55];
double dis(node a, node b) {
return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
int n, k;
struct node1 {
int x, y;
double dis;
bool operator<(const node1 a)const {
return dis < a.dis;
}
}edge[5555];
/*
最大团 = 补图G的最大独立集数
———>最大独立集数 = 补图G'最大团
*/
#define N 55
int mx;//最大团数(要初始化为0)
int x[N], tuan[N];
int can[N][N];//can[i]表示在已经确定了经选定的i个点必须在最大团内的前提下还有可能被加进最大团的结点集合
int num[N];//num[i]表示由结点i到结点n构成的最大团的结点数
bool g[N][N];//邻接矩阵(从1开始)
bool dfs(int tot, int cnt) {
int i, j, k;
if (tot == 0) {
if (cnt > mx) {
mx = cnt;
for (i = 0; i < mx; i++) tuan[i] = x[i];
return true;
}
return false;
}
for (i = 0; i < tot; i++) {
if (cnt + (tot - i) <= mx)return false;
if (cnt + num[can[cnt][i]] <= mx)return false;
k = 0;
x[cnt] = can[cnt][i];
for (j = i + 1; j < tot; j++) {
if (g[can[cnt][i]][can[cnt][j]]) can[cnt + 1][k++] = can[cnt][j];
}
if (dfs(k, cnt + 1))return false;
}
return false;
}
void MaxTuan() {
int i, j, k;
mx = 1;
for (i = n; i >= 1; i--) {
k = 0;
x[0] = i;
for (j = i + 1; j <= n; j++) {
if (g[i][j]) can[1][k++] = j;
}
dfs(k, 1);
num[i] = mx;
}
}
int nu;
bool check(double mid) {
int pos = lower_bound(edge + 1, edge + nu + 1, node1{ 0, 0, mid }) - edge;
memset(g, 0, sizeof(g));
for (int i = pos; i <= nu; i++) {
g[edge[i].x][edge[i].y] = g[edge[i].y][edge[i].x] = 1;
}
mx = 0;
MaxTuan();
if (num[1] >= k) return true;
else return false;
}
int main() {
while (~scanf("%d%d", &n, &k)) {
for (int i = 1; i <= n; i++) scanf("%lf%lf", &p[i].x, &p[i].y);
nu = 0;
double l = 0, r = 0;
for (int i = 1; i <= n; i++) {
for (int j = i + 1; j <= n; j++) {
double di = dis(p[i], p[j]);
edge[++nu] = { i, j, di };
r = max(r, di);
}
}
sort(edge + 1, edge + nu + 1);
double mid;
while (r - l > 1e-6) {
mid = (l + r) / 2;
if (check(mid))l = mid;
else r = mid;
}
printf("%.2f\n", mid);
}
return 0;
}