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;
}