HDU 3585

这道题说实话。。如果说跟最大团找关系。。我看不出来。。。看来对最大团理解还是太浅了。。。

 

题意：

有 n 个点 现在给出这 n 个点的坐标 x y ，选择 k 个点做一个集合，使得集合内部最近的两个点 距离最大。

 

n范围不超过 50 证明 最大团的 dfs 加剪枝 可以跑的过。。

首先 我们把所有点之间的距离存在一起，然后排序，既然答案是存在的，那么就代表 可以二分

我们对排序好的距离进行二分，之后开始跑最大团，一旦集合内存在 k 个点使得 dis [ i ] [ j ]  >  long [ mid ] 向上二分，反之向下二分，最终找出结果。

 

既然是选择 k 个点，那么这 k 个点肯定能组成一集合，我们只把距离 大于二分结果的点进行联通，之后找出整个图里的最大团，再计算一下最大团里节点数目是不是恰好等于 k 来决定二分方向。。最终得出结果

 

 

 

以下是 AC 代码

 

 

 

 

 

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn = 100;
int n,k;
bool map[maxn][maxn];//确认联通图
int mx,dp[maxn];// dp 优化
int stk[maxn][maxn];//求最大团
double dis[maxn][maxn],l[maxn*maxn];
double x[maxn],y[maxn];
double dist(double a,double b,double c,double d)
{
    double x=a-c;
    double y=b-d;
    return sqrt(x*x+y*y);
}
bool dfs(int now,int dis)
{
	if(now == 0)
	{
		mx = max(mx, dis);
		return true;
	}
	for(int i=0; i<now; i++)
	{
		int u = stk[dis][i];
		if(dis + n - u + 1 <= mx || dis + dp[u] <= mx)//剪枝
		return false;
		int cnt = 0;
		for(int j=i+1; j<now; j++)
		{
			int v = stk[dis][j];
			if(map[u][v])
			{
				stk[dis+1][cnt++] = v;
			}
		}
		dfs(cnt, dis+1);
	}
	return false;
}
bool judge(int mid)
{
    memset(map,false,sizeof(map));
    for(int i=1;i<=n;i++)
    {
        for(int j=i+1;j<=n;j++)
        {
            if(dis[i][j]>=l[mid])
            {
                map[i][j]=map[j][i]=true;
            }
        }
    }
    mx=0;
    int ns;
    for(int i=n;i>=1;i--)
    {
        ns=0;
        for(int j=i+1;j<=n;j++)
        {
            if(map[i][j])
            {
                stk[1][ns++]=j;
            }
        }
        dfs(ns,1);
        dp[i]=mx;
        if(mx>=k)
			return true;
    }
    return false;
}
int main()
{
    while(~scanf("%d%d",&n,&k))
    {
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf",&x[i],&y[i]);
        }
        int cnt=0;
        for(int i=1;i<=n;i++)
        {
            for(int j=i+1;j<=n;j++)
            {
                dis[j][i] = dis[i][j] = dist(x[i],y[i],x[j],y[j]);
                l[++cnt] = dis[i][j];
            }
            dis[i][i] = 0;
        }
        sort(l+1,l+cnt+1);
        int lt=1,rt=cnt,ans;
        while(lt <= rt)
        {
            int mid=(lt+rt)>>1;
            if(judge(mid))
            {
            	ans = mid;
                lt=mid + 1;
            }
            else
            {
		    	rt=mid - 1;
		    }
        }
        printf("%.2lf\n",l[ans]);
    }
    return 0;
}