简单的两圆关系

$题目地址$

有一定总结性的题目，所以我记下来了。

题意就是对于三维中的n个球找到一条最小生成树，相碰或者想包含的小球视为可以直接到达。所以这样的情况直接让边长为$0$即可。那么怎么判断这样的情况呢。

$d$为圆心距，$R$为大圆半径，$r$为小圆半径

$d=R-r$时，两圆内切

$R-r<d<R+r$时，两圆相交

$d=R+r$时，两圆外切

$d>R+r$时，两圆不相交

那么直接求就可以了。唯一需要注意的是：$G++$提交的时候double输出应该用$\%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);
    }
}