题目描述每一个要上天堂的人都要经历一番考验，当然包括小X，小X开始了他进入天堂的奇异之旅。地狱有18层，天堂竟然和地狱一样，也有很多很多层，天堂共有N层。从下到上依次是第1，2，3，…，N层，天堂的每一层都是一个延伸无限远的地板，在地板上人可以任意走动，层与层之间是平行关系，每一层的地板都是由人不能穿过的物质构成，幸好每一层地板上有且仅有1个人可以通过的洞口。我们可以把小X和洞口，还有下面提到的气球店都看成点，坐标是二维的。小X开始在第31层的（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 40 0 1 29 00 1 12 3 10 1 21 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.