HDU 3585 maximum shortest distance （二分+最大团）

题意：有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;
}