poj2031

#include<iostream>
#include<algorithm>
#include<cstring>
#include<math.h>
#define INF 0x3f3f3f3f
using namespace std;

const int N=110;
double g[N][N];
double dist[N];
int visited[N];
int n;

struct Range{
    double x,y,z,r;
}range[N];

double calculate( Range &a, Range &b){
    double v;
    v=sqrt((pow((a.x-b.x),2)+pow((a.y-b.y),2)+pow((a.z-b.z),2)))-a.r-b.r;
    
    if(v>0){
    return v;
    }
    else
    return 0;    
}
void init()
{
    for(int l=0;l<n;l++){
        visited[l]=0;
        dist[l]=INF;
}
}
double prim(){
    
    
    double res=0;
    for(int i=0;i<n;i++){
        int t=-1;
        dist[0]=0;
        for(int j=0;j<n;j++)
        if(!visited[j]&&(t==-1||dist[j]<dist[t]))
        t=j;
        if(t==-1)
    break;
        if(i) 
        res+=dist[t];
        visited[t]=1;
        for(int k=0;k<n;k++){
        
        dist[k]=min(dist[k],g[t][k]);

    }
        
        
    
    }
    return res;
    
    
} 

int main(){
    while(scanf("%d",&n),n){
        memset(g, 0x3f, sizeof g);
        init();
        for(int i=0;i<n;i++){
            
            scanf("%lf%lf%lf%lf",&range[i].x,&range[i].y,&range[i].z,&range[i].r);
        }
        
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++){
                g[i][j]=calculate(range[i],range[j]);
            
            }
        }
        double t=prim();
        printf("%.3f\n",t);
        
    }
    return 0;
}