百度之星 / 初赛第二场 B题

怎么说呢，只能说自己还不够熟练，能力还不够，细心成都还不够吧。这样的二分题目在POJ的训练计划里面有类似的题目，自己也是都刷了，可是在调这道题目的时候废了老大的尽了，比赛结束后听polla的一句<t才属于同一类，才知道自己怎么犯了这么二逼的错误呢！唉。。。。自己给弄成>t了。改过后终于出来了。

思路：找出最大的t，如果按t分类，可分出1组，按0分类分出n组，二分枚举t然后利用并查集判断分出多少种类。

#include <cstdio>

#include <cstring>

#include <iostream>

#define maxn 1007



using namespace std;



const double eps = 1e-8;

int n,m,k;

double t;

int f[maxn];

struct node

{

    double x,y,z;

}p[maxn];





int cmp(double x)

{

    if (x > eps) return 1;

    else if (x < -eps) return -1;

    else return 0;

}

int find(int x)

{

    if (x != f[x])

    f[x] = find(f[x]);

    return f[x];

}



void Union(int x,int y)

{

    x = find(x);

    y = find(y);

    if (x != y)

    f[y] = x;

}



double getse(int a,int b)

{

    double x = p[b].x - p[a].x;

    double y = p[b].y - p[a].y;

    double z = p[b].z - p[a].z;

    return (x*x + y*y + z*z);

}

//并查集判断是否属于同一类

int getnum(double tx)

{

    int i,j;

    for (i = 0; i < n; ++i) f[i] = i;

    for (i = 0; i < n; ++i)

    {

        for (j = i + 1; j < n; ++j)

        {

            double tmp = getse(i,j);

            if (cmp(tmp - tx) < 0) //就是这里的二逼错误让我调了很长时间

            Union(i,j);

        }

    }

    int ct = 0;

    for (i = 0; i < n; ++i)

    {

        if (i == find(i)) ct++;

    }

    return ct;

}

int main()

{

    int i,j;

    t = 0;

    scanf("%d%d",&n,&k);

    for (i = 0; i < n; ++i) scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].z);

    double pt;

    for (i = 0; i < n; ++i)

    {

        for (j = 0; j < n; ++j)

        {

            if (i != j)

            {

                pt = getse(i,j);

                if (pt > t) t = pt;

            }

        }

    }

    //printf("%lf\n",t);

    double l = 0;

    double r = t;

    double mid = 0;

    //二分枚举

    while (cmp(l - r) < 0)

    {

        mid = (l + r)/2.0;

        int tt = getnum(mid);

        //printf(">> %lf %lf %lf %d\n",l,r,mid,tt);

        if (tt >= k)  l = mid;

        else r = mid;

    }

    printf("%.6lf\n",l);

}