【POJ】[2031]Building a Space Station
给出球心的位置（x,y,z）和半径(r)，如果两球相离，那么两球之间的距离是两球心间的距离减去两半径，否则两球之间的距离为0
要求所有的点都相连并输出最小的距离和
感觉就是3D版的
【杭电】[1875]畅通工程再续
不过好像比那个还简单些
思路还是根据点求出来边长
注意的是如果两圆不相离 那么距离为0
相离了是两表面的最短距离
记得两点距离减去半径
英语还是硬伤啊……
还有就是不知道为什么C++能过……
G++好多都WA了……
看算法……
也暂时不考虑这个问题了
#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
struct xyz {
double x,y,z,r;
} a[120];
int main() {
double dis[120];
double e[120][120];
bool flag[120];
double inf=999999;
int N;
while(scanf("%d",&N),N) {
for(int i=1; i<=N; i++) {
dis[i]=inf;
flag[i]=false;
}
for(int i=1; i<=N; i++) {
for(int j=1; j<=N; j++) {
e[i][j]=inf;
}
}
for(int i=1; i<=N; i++) {
scanf("%lf %lf %lf %lf",&a[i].x,&a[i].y,&a[i].z,&a[i].r);
for(int j=1; j<=N; j++) {
if(i==j)
continue;
else {
double l=sqrt(pow(a[i].x-a[j].x,2)+pow(a[i].y-a[j].y,2)+pow(a[i].z-a[j].z,2));
if(l>a[i].r+a[j].r) {
e[j][i]=e[i][j]=l-a[i].r-a[j].r;
} else {
e[j][i]=e[i][j]=0;
}
}
}
}
double sum=0;
dis[1]=0;
while(true) {
int v=-1;
for(int u=1; u<=N; u++) {
if(!flag[u]&&(v==-1||dis[u]<dis[v]))
v=u;
}
if(v==-1)
break;
flag[v]=true;
sum+=dis[v];
for(int u=1; u<=N; u++) {
dis[u]=min(dis[u],e[v][u]);
}
}
printf("%.3lf\n",sum);
}
return 0;
}
