Description
你是一艘战列巡洋舰的引擎操作人员,这艘船的船员在空间中侦测到了一些无法辨识的异常信号。你的指挥官给你下达了命令,让你制定航线,驾驶战列巡洋舰到达那里。
船上老旧的曲速引擎的速度是0.1AU/s。然而,在太空中分布着许多殖民星域,这些星域可以被看成一个球。在星域的内部,你可以在任何地方任意次跳跃到星域内部的任意一个点,不花费任何时间。
你希望算出到达终点的最短时间。
Input
输入包含多组测试数据。
对于每一组数据,第一行包含一个正整数n,表示殖民星域的数量。
接下来n 行,第i 行包含四个整数Xi,Yi,Zi,Ri,表示第i个星域的中心坐标为(Xi, Yi,Zi),星域的半径是Ri。
接下来两行,第一行包含值Xa,Ya,Za,告诉你当前坐标为(Xa, Ya,Za)。
第二行包含值Xo,Yo,Zo,告诉你目的地坐标为(Xo, Yo,Zo)。
输入以一行单独的-1 结尾。所有坐标的单位都是天文单位(AU)。
Output
对于每一组输入数据,输出一行表示从目前的位置到达指定目的地的最短时间,取整到最近整数。输入保证取整是明确的。
Sample Input
1
20 20 20 1
0 0 0
0 0 10
1
5 0 0 4
0 0 0
10 0 0
-1
Sample Output
100
20
Data Constraint
每个输入文件至多包含10 个测试数据。
对于10% 的数据,n = 0。
对于30% 的数据,0<=n<=10。
对于100% 的数据,0<=n<=100,所有坐标的绝对值<=10000 ,半径r<=10000。
你可以认为,你所在的星区的大小为无限大。
简单的spfa建边即可。
#pragma GCC optimize(3)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int n;
long double ans,map[102][102];
struct star
{
int x,y,z,r;
}a[101];
long double dis(star a,star b)
{
int t1=abs(a.x-b.x),t2=abs(a.y-b.y),t3=abs(a.z-b.z);
long double t4=sqrt(1.0*t1*t1+1.0*t2*t2),t5=sqrt(t4*t4+1.0*t3*t3);
return t5;
}
int main()
{
freopen("warp.in","r",stdin);
freopen("warp.out","w",stdout);
while(true)
{
scanf("%d",&n);
if(n==-1)break;
memset(a,0,sizeof(a));
memset(map,0,sizeof(map));
for(int i=1;i<=n;i++)
scanf("%d%d%d%d",&a[i].x,&a[i].y,&a[i].z,&a[i].r);
scanf("%d%d%d",&a[0].x,&a[0].y,&a[0].z);
scanf("%d%d%d",&a[n+1].x,&a[n+1].y,&a[n+1].z);
a[0].r=a[n+1].r=0;
for(int i=0;i<=n+1;i++)
for(int j=0;j<=n+1;j++)
if(i!=j)
{
map[i][j]=dis(a[i],a[j])-a[i].r-a[j].r;
if(map[i][j]<0)map[i][j]=0;
}
for(int k=0;k<=n+1;k++)
for(int i=0;i<=n+1;i++)
for(int j=0;j<=n+1;j++)
if(k!=i&&k!=j&&i!=j)
if(map[i][k]+map[k][j]<map[i][j])
map[i][j]=map[i][k]+map[k][j];
map[0][n+1]*=10;
if((int)(map[0][n+1]*10)%10>=5)printf("%d\n",(int)(map[0][n+1]+1));
else printf("%d\n",(int)map[0][n+1]);
}
return 0;
}