POJ 2031 Building a Space Station(最小生成树)

POJ 2031 Building a Space Station(最小生成树)

http://poj.org/problem?id=2031

题意:

       有n个三维的球体,对于每个球体给出了它的圆心坐标(x,y,z)以及它的半径r. 现在要你用最短的线使得所有球体属于同一个连通分量.

如果两个球体在初始的时候已经相切或重叠了,就认为它们已经连接上了,不用另外连线.否则需要连接长度= 圆心距-r1-r2 的线.

分析:

       其实就是给了你n个点,你要算他们初始是否连通以及不连通的情况下的初始长度. 然后你需要求出该图的最小生成树的权值.

       注意,如果初始时两个球心已经连通了,那么我们把这两个点的权值设为0即可.

AC代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=100+10;
const int maxm=100*100+10;

struct Edge
{
    int u,v;
    double dist;
    Edge(){}
    Edge(int u,int v,double d):u(u),v(v),dist(d){}
    bool operator<(const Edge&rhs)const
    {
        return dist <rhs.dist;
    }
};

struct Kruskal
{
   int n,m;     //点数,边数
   int fa[maxn];//并查集
   Edge edges[maxm];

   void init(int n)
   {
       this->n=n;
       m=0;
       memset(fa,-1,sizeof(fa));
   }

   int findset(int x){ return fa[x]==-1?x: fa[x]=findset(fa[x]); }

   void AddEdge(int u,int v,double dist)
   {
       edges[m++]=Edge(u,v,dist);
   }

   double kruskal()
   {
       double sum=0;//最小生成树权值
       int cnt=0;   //最小生成树边数目
       sort(edges,edges+m);
       for(int i=0;i<m;i++)
       {
           int u=edges[i].u, v=edges[i].v;
           if(findset(u) != findset(v))
           {
               sum +=edges[i].dist;
               fa[findset(u)] = findset(v);
               if(++cnt>=n-1) break;
           }
       }
       return sum;
   }
}KK;

struct Point
{
    double x,y,z,r;
}p[maxn];

double get_dist(int i,int j)//得到第i个球心与第j个秋心需要连的线长度
{
    double len = sqrt((p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y)+(p[i].z-p[j].z)*(p[i].z-p[j].z));
    double ans = len-p[i].r-p[j].r;
    if(ans<0) return 0;
    return ans;
}

int main()
{
    int n;
    while(scanf("%d",&n)==1&&n)
    {
        for(int i=0;i<n;i++)
            scanf("%lf%lf%lf%lf",&p[i].x,&p[i].y,&p[i].z,&p[i].r);
        KK.init(n);
        for(int i=0;i<n;i++)
        for(int j=i+1;j<n;j++)
            KK.AddEdge(i,j,get_dist(i,j));
        printf("%.3lf\n",KK.kruskal());
    }
    return 0;
}