warp

3517. 【NOIP2013模拟11.6B组】空间航行(warp) 
(File IO): input:warp.in output:warp.out

Time Limits:  1000 ms  Memory Limits:  262144 KB  Detailed Limits  

Goto ProblemSet

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。

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

uses math;
var
        a,b,c,d,e,i,j,k:longint;
        f:array[0..101,0..101]of real;
        v:array[0..101,1..4]of longint;
begin
        assign(input,'warp.in');reset(input);
        assign(output,'warp.out');rewrite(output);
        while not eof do
        begin
                readln(a);
                if a=-1 then
                        halt;
                fillchar(v,sizeof(v),0);
                for b:=1 to a do
                        readln(v[b,1],v[b,2],v[b,3],v[b,4]);
                readln(v[0,1],v[0,2],v[0,3]);
                readln(v[a+1,1],v[a+1,2],v[a+1,3]);
                fillchar(f,sizeof(f),0);
                for b:=0 to a+1 do
                begin
                        for c:=0 to a+1 do
                        begin
                                f[b,c]:=sqrt(sqr(v[b,1]-v[c,1])+sqr(v[b,2]-v[c,2])+sqr(v[b,3]-v[c,3]));
                                f[b,c]:=f[b,c]-v[b,4]-v[c,4];
                                if f[b,c]<0 then f[b,c]:=0;
                        end;
                end;
                for k:=0 to a+1 do
                begin
                        for i:=0 to a+1 do
                        begin
                                for j:=0 to a+1 do
                                begin
                                        f[i,j]:=min(f[i,k]+f[k,j],f[i,j]);
                                end;
                        end;
                end;
                writeln((f[0,a+1]*10):0:0);
        end;
        close(input);
        close(output);
end.