Description
【问题背景】
LHX教主最近总困扰于前来膜拜他的人太多了,所以他给他的花园加上了一道屏障。
【问题描述】
可以把教主的花园附近区域抽像成一个正方形网格组成的网络,每个网格都对应了一个坐标(均为整数,有可能为负),若两个网格(x1, y1),(x2, y2)有|x1 – x2| + |y1 – y2| = 1,则说这两个网格是相邻的,否则不是相邻的。
教主在y = 0处整条直线上的网格设置了一道屏障,即所有坐标为(x, 0)的网格。当然,他还要解决他自己与内部人员的进出问题,这样教主设置了N个入口a1, a2, …, aN可供进出,即对于y = 0上的所有网格,只有 (a1, 0),(a2, 0), ……, (aN, 0) 可以通过,之外的所有纵坐标为0的网格均不能通过,而对于(x, y)有y不为0的网格可以认为是随意通过的。
现在教主想知道,给定M个点对(x1, y1),(x2, y2),并且这些点均不在屏障上,询问从一个点走到另一个点最短距离是多少,每次只能从一个格子走到相邻的格子。
Input
输入的第1行为一个正整数N,为屏障上入口的个数。
第2行有N个整数,a1, a2, …, aN,之间用空格隔开,为这N个入口的横坐标。
第3行为一个正整数M,表示了M个询问。
接下来M行,每行4个整数x1, y1, x2, y2,有y1与y2均不等于0,表示了一个询问从(x1, y1)到(x2, y2)的最短路。
Output
输出共包含m行,第i行对于第i个询问输出从(x1, y1)到(x2, y2)的最短路距离是多少。
Sample Input
2
2 -1
2
0 1 0 -1
1 1 2 2
Sample Output
4
2
Data Constraint
Hint
【数据规模】
对于20%的数据,有n,m≤10,ai,xi,yi绝对值不超过100;
对于40%的数据,有n,m≤100,ai,xi,yi绝对值不超过1000;
对于60%的数据,有n,m≤1000,ai,xi,yi绝对值不超过100000;
对于100%的数据,有n,m≤100000,ai,xi,yi绝对值不超过100000000。
.
.
.
.
.
.
.
.
分析
题目大意就是在一个平面直角坐标系中x轴上是一道屏障,但上面有n个开头,要求两个点的距离
显然,这题可以分类讨论来做
①两个点的y坐标同号,直接求曼哈顿距离输出
如果不行,则二分一个在中间的开口
②如果这个开口位于两个点中间,求出两个点分别到开口的曼哈顿距离相加并输出
③否则判断一下它距离哪一个端点近,就由哪里绕过去
虽然此时的pos值不一定准确,但十分接近,真实值一定出现在pos-1,pos,pos+1中,求最小值输出就好了
.
.
.
.
.
.
.
程序:
uses math;
var
x1,y1,x2,y2,i,n,m,pos:longint;
a:array[-1..100002]of longint;
t1,t2,t3:int64;
procedure swap(x,y:longint);
var
xc:longint;
begin
xc:=x1;x1:=x2;x2:=xc;
end;
function find(x:longint):longint;
var
l,r,mid,k:longint;
begin
l:=1;r:=n;k:=n;
while l<=r do
begin
mid:=(l+r) div 2;
if a[mid]<x then l:=mid+1 else
begin
r:=mid-1;
k:=mid;
end;
end;
exit(k);
end;
procedure kp(l,r:longint);
var
i,j,mid:longint;
begin
if l>=r then exit;
i:=l;j:=r;mid:=a[(i+j) div 2];
repeat
while a[i]<mid do inc(i);
while a[j]>mid do dec(j);
if i<=j then
begin
a[-1]:=a[i];a[i]:=a[j];a[j]:=a[-1];
inc(i);dec(j);
end;
until i>j;
kp(l,j);
kp(i,r);
end;
begin
readln(n);
for i:=1 to n do
read(a[i]);
kp(1,n);
readln;
readln(m);
for i:=1 to m do
begin
readln(x1,y1,x2,y2);
if (y1<0)and(y2<0)or(y1>0)and(y2>0) then
begin
writeln(abs(x1-x2)+abs(y1-y2));
continue;
end;
if x1>x2 then swap(x1,x2);
pos:=find((x1+x2) div 2);
if (a[pos]>=x1)and(x2>=a[pos]) then writeln(abs(x1-x2)+abs(y1-y2)) else
begin
t1:=abs(a[pos]-x1)+abs(a[pos]-x2)+abs(y1-y2);
t2:=200000000;
t3:=200000000;
if pos>1 then t2:=abs(a[pos-1]-x1)+abs(a[pos-1]-x2)+abs(y1-y2);
if pos<n then t3:=abs(a[pos+1]-x1)+abs(a[pos+1]-x2)+abs(y1-y2);
t1:=min(t1,t2);
t1:=min(t1,t3);
writeln(t1);
end;
end;
end.