POJ - 2031

题意：给你n个球，如果两个球是接触的那么连接就不花费，如果不是接触的那么就是接触的距离，很容易看出是个最小生成树的。

最近对ACM有点迷茫，不知道为什么，感觉看到每个题都感觉很难，也不愿意去思考，有点脑子疼。。很无奈，不知怎么取舍，毕竟我的专业不是专攻这个的。。哎，有点无奈。

#include <iostream>
#include <math.h>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
#include <iomanip>
using namespace std;
const int maxn = 20000 + 20;
struct node
{
    int come, to;
    double cost;
}edge[maxn];
double x[maxn], y[maxn], z[maxn], r[maxn];
int  fa[maxn];
void init()
{
    for(int i = 0; i < maxn; i++)
        fa[i] = i;
}
int getf(int x)
{
    if(x == fa[x])
        return x;
    return fa[x] = getf(fa[x]);
}
double distance_(int i, int j)
{
    return sqrt((x[i] - x[j]) * (x[i] - x[j]) + (y[i] - y[j]) * (y[i] - y[j]) + (z[i] - z[j]) * (z[i] - z[j]));
}
int same(int x, int y)
{
    x = getf(x);
    y = getf(y);
    if(x == y)
        return 1;
    fa[x] = y;
    return 0;
}
bool cmp(node p1, node p2)
{
    return p1.cost < p2.cost;
}
double kursual(int sum)
{
    sort(edge, edge + sum, cmp);
    init();
    double ans = 0;
    for(int i = 0; i < sum; i++)
    {
        if(!same(edge[i].come, edge[i].to))
            ans += edge[i].cost;

    }
    return ans;
}
int main()
{
    int n;
    ios::sync_with_stdio(false);
    while(cin >> n && n)
    {
        for(int i = 0; i < n; i++)
            cin >> x[i] >> y[i] >> z[i] >> r[i];
        int sum = 0;
        for(int i = 0; i < n - 1; i++)
        {
            for(int j = i + 1; j < n; j++)
            {

                double dis = distance_(i ,j);
                if(dis <= r[i] + r[j])
                {
                    edge[sum].come = i;
                    edge[sum].to = j;
                    edge[sum++].cost = 0;
                }
                else
                {
                    edge[sum].come = i;
                    edge[sum].to = j;
                    edge[sum++].cost = dis - r[i] - r[j];
                }
            }
        }
        double ans = kursual(sum);
        cout <<setprecision(3) << fixed <<ans << endl;
    }
    return 0;
}