jzoj 3517. 【NOIP2013模拟11.6B组】空间航行(warp)

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。

你可以认为，你所在的星区的大小为无限大。

Solution

易知，在两个球的之间航行的时间等于两个球心之间的距离减去两个球的半径。 将每个球的球心当作点建图，然后运行最短路算法即可。 无论是 Floyd 算法还是 Dijkstra 算法都可以在限定时间内解决。需要注意的是不要弄出负 的航行时间。

Code

#include<cstdio>
#include<cmath>
#include<iostream>
#include<algorithm>
int n;
double map[110][110];
struct cj
{
	int x,y,z,r;
}st[21000];
using namespace std;
int main()

{
	freopen("warp.in","r",stdin);
	freopen("warp.out","w",stdout);
	while(cin>>n)
	{
		if(n==-1)
		return 0;
		for(int i=1;i<=n;i++)
		scanf("%d%d%d%d",&st[i].x,&st[i].y,&st[i].z,&st[i].r);
		scanf("%d%d%d",&st[0].x,&st[0].y,&st[0].z);
		scanf("%d%d%d",&st[n+1].x,&st[n+1].y,&st[n+1].z);
		st[0].r=0;st[n+1].r=0;
		for(int i=0;i<=n+1;i++)
		{
			for(int j=0;j<=n+1;j++)
			{
				if(i!=j)
				{
					double k=sqrt((st[i].x-st[j].x)*(st[i].x-st[j].x)+(st[i].y-st[j].y)*(st[i].y-st[j].y)+(st[i].z-st[j].z)*(st[i].z-st[j].z))-st[i].r-st[j].r;
					map[i][j]=max(0.0,k);
				}
			}	
		}
		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][j]>map[i][k]+map[k][j])
					{
						map[i][j]=map[i][k]+map[k][j];
					}
				}
			}
		}	
		double zz=map[0][n+1]*10;
		int kk=map[0][n+1]*10;
		if(zz-kk>=0.5)
		kk=kk+1;
		printf("%d\n",kk);
	}
}