题目描述
每一个要上天堂的人都要经历一番考验，当然包括小X，小X开始了他进入天堂的奇异之旅。地狱有18层，天堂竟然和地狱一样，也有很多很多层，天堂共有N层。从下到上依次是第1，2，3，…，N层，天堂的每一层都是一个延伸无限远的地板，在地板上人可以任意走动，层与层之间是平行关系，每一层的地板都是由人不能穿过的物质构成，幸好每一层地板上有且仅有1个人可以通过的洞口。
我们可以把小X和洞口，还有下面提到的气球店都看成点，坐标是二维的。小X开始在第
3
1层的（0，0）.
小X的重量为M，第i层与第i+1层之间的特殊气体能浮起的重量为Wi ,每一层的地面上散落了若干个气球店，多个气球店可以在同一点，每个气球可以浮起的重量是1，去一个气球店一次只能领取一个气球，不能连续在一个气球店领取气球，当然你可以在两个气球店之间来回跑，每个气球店供应的气球都是无限多的。第i层的气球只能在第i层进入第i+1层时使用，当小X在第i层，只有站到了第i+1层洞口的位置（在其它位置不会浮起），并且自身的重量小于等于气球和特殊气体浮起重量的总和，才可以进入第i+1层。小X想知道他要到达第N层走过的长度最少是多少？题目保证有解。
输入文件
第1行： 三个正整数N，M，Q（Q表示气球店）
第2行： 共2*（N-1）个整数，每两个数描述1个洞口坐标，第i对xi,yi表示第i+1层的洞口位置(xi,yi)。
第3行： 共N-1个整数，第i个数为Wi。
往后Q行，每行三个整数x,y,z , 表示第Z层有一个气球店，坐标为（x,y）
输出文件
1个实数L，保留两位小数，表示小X最少要走的长度。
样例输入
3 10 4
0 0 1 2
9 0
0 1 1
2 3 1
0 1 2
1 1 2
样例输出
13.00
注释
【样例解释】
在第一层从（0，0）出发到（0，1）取得1个气球并返回（0，0）即可到达第二层。长度：2.00
在第二层，从（0，0）到（0，1）领取气球，再到（1，1）领取气球，两个点来回跑，第5次到达（1，1）时恰好气球数达到10，走到（1，2）即可到达第3层终点。长度：11.00
总长度：13.00
【数据范围】
2<=N<=100
每层的气球店数目不超过50。
0<=M<=100, 0<=Wi<=100
坐标-3000<=x,y<=3000
可以看出，每一层与其它层是互不影响的。
对于第x层，共有number[x]个气球店，气体浮力为air[x]，即需领m-air[x]个气球
f[i,j]表示在第i个气球点领取第j个气球
f[i,j]:=min(f[k,j-1]+dist[i,j]) {i<>k}
f[i,1]:=起点到每一个气球店的距离
第x层最小的移动距离为min{f[i,m-air[x]]+终点到i的距离}
and 不能再把数据范围看错了
code：
var n,m,q:longint;
pole:array[1..200]of record
x,y:longint;
end;
air:array[1..200]of longint;
number:array[1..200]of longint;
position:array[1..200,1..60]of record
x,y:longint;
end;
f:array[1..1010,1..1000]of double;
dist1,dist2:array[1..100]of double;
dist:array[1..100,1..100]of double;
ii,i,j,k:longint;
sx,sy,fx,fy:longint;
x,y,z:longint;
min,ans:double;
temp:double;
begin
readln(n,m,q);
for i:=2 to n do
read(pole[i].x,pole[i].y);
for i:=1 to n-1 do
read(air[i]);
fillchar(position,sizeof(position),0);
fillchar(number,sizeof(number),0);
for i:=1 to q do
begin readln(x,y,z);
inc(number[z]);
position[z,number[z]].x:=x;
position[z,number[z]].y:=y;
end;
ans:=0;
sx:=pole[1].x;
sy:=pole[1].y;
for ii:=1 to n-1 do
begin if air[ii]<m
then begin fillchar(dist1,sizeof(dist1),0);
fillchar(dist2,sizeof(dist2),0);
fillchar(dist,sizeof(dist),0);
for i:=1 to number[ii] do
dist1[i]:=sqrt(sqr(position[ii,i].x-sx)+
sqr(position[ii,i].y-sy));
for i:=1 to number[ii] do
for j:=1 to number[ii] do
dist[i,j]:=sqrt(sqr(position[ii,i].x-position[ii,j].x)+
sqr(position[ii,i].y-position[ii,j].y));
for i:=1 to number[ii] do
for j:=1 to m-air[ii] do
f[i,j]:=maxint*5000;
for i:=1 to number[ii] do
f[i,1]:=dist1[i];
fx:=pole[ii+1].x;
fy:=pole[ii+1].y;
for i:=1 to number[ii] do
dist2[i]:=sqrt(sqr(position[ii,i].x-fx)+
sqr(position[ii,i].y-fy));
for j:=2 to m-air[ii] do
for i:=1 to number[ii] do
for k:=1 to number[ii] do
if (f[i,j]>f[k,j-1]+dist[i,k])and(i<>k)
then f[i,j]:=f[k,j-1]+dist[i,k];
min:=maxlongint;
for i:=1 to number[ii] do
if min>f[i,m-air[ii]]+dist2[i]
then min:=f[i,m-air[ii]]+dist2[i];
ans:=ans+min;
sx:=pole[ii+1].x;
sy:=pole[ii+1].y;
end
else begin fx:=pole[ii+1].x;
fy:=pole[ii+1].y;
ans:=ans+sqrt(sqr(sx-fx)+sqr(sy-fy));
sx:=pole[ii+1].x;
sy:=pole[ii+1].y;
end;
end;
writeln(ans:0:2);
end.
