题目描述每一个要上天堂的人都要经历一番考验,当然包括小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.
时间: 2024-08-04 14:23:53