//这道题并不难,注意细节就能一次AC//
AC代码:
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
using namespace std;
#define N 100000
int p[N],n;
int L;
struct node
{
double r,x,y,z,c,d;
int k,k1,k2;
}t[N],T[N];
double cmp(node p,node q)
{
return p.c<q.c;
}
int find(int k)
{
if(p[k]==k)
{
return p[k];
}
else
{
return find(p[k]);
}
}
double cl()
{
double sum=0;
int cnt=0;
int i;
for(i=0;i<L;i++)
{
int x=find(T[i].k1);
int y=find(T[i].k2);
if(x!=y)
{
p[y]=x;
cnt++;
sum+=T[i].c;
}
if(cnt==n-1)
return sum;
}
return sum;
}
int main()
{
while(scanf("%d",&n)!=EOF&&n)
{
int i,j;
L=0;
for(i=0;i<n;i++)
{
scanf("%lf%lf%lf%lf",&t[i].x,&t[i].y,&t[i].z,&t[i].r);
t[i].k=i+1;
}
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
double l=sqrt((t[i].x-t[j].x)*(t[i].x-t[j].x)+(t[i].y-t[j].y)*(t[i].y-t[j].y)+(t[i].z-t[j].z)*(t[i].z-t[j].z));
double D=t[i].r+t[j].r;
if(l>D)
{
T[L].c=l-D;
}
else
{
T[L].c=0;
}
T[L].k1=t[i].k;
T[L].k2=t[j].k;
L++;
}
p[i]=i;
}
p[n]=n;
sort(T,T+L,cmp);
/* for(i=0;i<L;i++)
{
printf("%d %d %d\n",T[i].k1,T[i].k2,T[i].c);
}*/
printf("%0.3lf\n",cl());
}
return 0;
}