//这道题并不难,注意细节就能一次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; }