有一定总结性的题目,所以我记下来了。
题意就是对于三维中的n个球找到一条最小生成树,相碰或者想包含的小球视为可以直接到达。所以这样的情况直接让边长为
0
0
0即可。那么怎么判断这样的情况呢。
d
d
d为圆心距,
R
R
R为大圆半径,
r
r
r为小圆半径
d
=
R
−
r
d=R-r
d=R−r时,两圆内切
R
−
r
<
d
<
R
+
r
R-r<d<R+r
R−r<d<R+r时,两圆相交
d
=
R
+
r
d=R+r
d=R+r时,两圆外切
d
>
R
+
r
d>R+r
d>R+r时,两圆不相交
那么直接求就可以了。唯一需要注意的是:
G
+
+
G++
G++提交的时候double输出应该用
%
f
\%f
%f
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<map>
#include<vector>
#include<set>
#include<queue>
#include<cctype>
#include<string>
#include<cmath>
#define FOR(i,a,b) for(int i=a;i<=b;i++)
#define eps 1e-8
using namespace std;
int n;
const int maxn=500050;
const int maxm=1000100;
struct circle{
double cx,cy,cz;
double r;
circle(){}
circle(double _cx,double _cy,double _cz,double _r):cx(_cx),cy(_cy),cz(_cz),r(_r){}
}A[1010];
struct Edge{
int from,to;double dist;
Edge(){}
Edge(int _from,int _to,double _dist):from(_from),to(_to),dist(_dist){}
};
Edge ed[maxm];int he[maxn],ne[maxm],etop=1;
void insert(int u,int v,double w){ed[etop]=Edge(u,v,w);ne[etop]=he[u];he[u]=etop++;}
bool vis[maxn];double minc[maxn];double ans=0;
priority_queue<pair<double,int>,vector<pair<double,int> >,greater<pair<double,int> > >q;
void prim(){
memset(vis,0,sizeof(vis));
memset(minc,0x7f,sizeof(minc));
q.push(make_pair(0.0,1));
while(!q.empty()){
double c=q.top().first;
int now=ed[q.top().second].to;
q.pop();
if(!vis[now]){
vis[now]=true;
ans+=c;
for(int i=he[now];i;i=ne[i]){
Edge &e = ed[i];
if(minc[e.to]>e.dist&&!vis[e.to]){
minc[e.to]=e.dist;
q.push(make_pair(e.dist,i));
}
}
}
}
}
double cal(circle a,circle b){
double d=(a.cx-b.cx)*(a.cx-b.cx)+(a.cy-b.cy)*(a.cy-b.cy)+(a.cz-b.cz)*(a.cz-b.cz);
d=sqrt(d);
if(d-(a.r+b.r)>eps)return d-(a.r+b.r);
else return 0.0;
}
int main(){
while(cin>>n){
if(n==0)break;
memset(he,0,sizeof(he));
memset(ne,0,sizeof(ne));
etop=1,ans=0;
double x,y,z,r;
FOR(i,1,n){
scanf("%lf%lf%lf%lf",&x,&y,&z,&r);
A[i]=circle(x,y,z,r);
}
FOR(i,1,n)FOR(j,i+1,n)if(i!=j){
double tmp=cal(A[i],A[j]);
insert(i,j,tmp);
insert(j,i,tmp);
}
prim();
printf("%.3lf\n",ans);
}
}