教主的花园
(p1.pas/cpp/in/out)
【问题背景】
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),并且这些点均不在屏障上,询问从一个点走到另一个点最短距离是多少,每次只能从一个格子走到相邻的格子。
【输入格式】
输入的第1行为一个正整数N,为屏障上入口的个数。
第2行有N个整数,a1, a2, …, aN,之间用空格隔开,为这N个入口的横坐标。
第3行为一个正整数M,表示了M个询问。
接下来M行,每行4个整数x1, y1, x2, y2,有y1与y2均不等于0,表示了一个询问从(x1, y1)到(x2, y2)的最短路。
【输出格式】
输出共包含m行,第i行对于第i个询问输出从(x1, y1)到(x2, y2)的最短路距离是多少。
【样例输入】
2
2 -1
2
0 1 0 -1
1 1 2 2
【样例输出】
4
2
【数据规模】
对于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。
【考察点】
二分
【思路】
通过画图可知:所有可能的情况共3种:都在上面,都在下面,一上一下。
其中,目测可知:前两种情况可以直接横纵坐标分别相减,然后加起来即可(因为没有障碍)。后面一种稍微麻烦一点,通过画图分析:如果(x1,x2)区间上有通道,那么实际效果同前两种情况,如果没有,那么最短路一定是走最靠近(x1,x2)这个区间的某个通道。
从而问题转化成了找一个距离(x1,x2)最近的通道了,这里是坐标,很明显可以让那些通道有序,所以我们可以二分解决。
【提交情况】
1次AC
【经验】
用y1*y2<0来判断的总要中枪嘛……叫你秀高端……
【收获】
熟悉二分操作
ACCode:
Program p1;
Var
n,m:longint;
a:array[0..100010]oflongint;
Procedure terminate;
begin
close(input);
close(output);
halt;
end;
Procedure qsort; //排序,题目没说数据有序
procedure sort(l,r: longint);
var
i,j,x,y:longint;
begin
i:=l;
j:=r;
x:=a[(l+r)div 2];
repeat
whilea[i]<x do
inc(i);
whilex<a[j] do
dec(j);
ifnot(i>j) then
begin
y:=a[i];
a[i]:=a[j];
a[j]:=y;
inc(i);
j:=j-1;
end;
until i>j;
if l<j then sort(l,j);
if i<r then sort(i,r);
end;
begin
sort(1,n);
end;
Procedure init;
var
i:longint;
begin
assign(input,'p1.in');
assign(output,'p1.out');
reset(input);
rewrite(output);
readln(n);
for i:=1 to n do
begin
read(a[i]);
end;
a[0]:=99999999;
a[n+1]:=99999999; //端点要注意,不然后面计算会把这里算上去,答案就诡异了……当然,也可以在后面特判,这种方法比较懒= =
readln(m);
qsort;
//for i:=1 to n do
//write(a[i],'');
//writeln;
end;
Function divide(l,r,x:longint):longint; //二分查找,找一个能把X插进去的地方
var
m:longint;
begin
whilel<r do
begin
m:=(l+r)>>1;
if x<a[(l+r)>>1] then r:=m
else l:=m+1;
end;
exit(m);
end;
Procedure main;
var
wjmzbmr,x1,y1,x2,y2:longint;
function solve:longint;
var
id1,id2,min:longint;
begin
id1:=divide(1,n,x1);
//write('id1:',id1,'');
id2:=divide(1,n,x2); //记录下两个x的位置
//writeln('id2:',id2);
if id1<>id2 then exit(abs(x1-x2)+abs(y1-y2)) //如果不相同,则代表a中有数位于(x1,x2)区间内
else //否则就找附近的某个最优值,其实就是左右两个点就行了……不放心可以左右多枚举几个
begin
min:=$7f7f7f7f;
if abs(y1-y2)+abs(x1-a[id1])+abs(x2-a[id1])<min then min:=abs(y1-y2)+abs(x1-a[id1])+abs(x2-a[id1]);
//writeln('当前:',abs(y1-y2)+abs(x1-a[id1])+abs(x2-a[id1]));
if abs(y1-y2)+abs(x1-a[id1+1])+abs(x2-a[id1+1])<min then min:=abs(y1-y2)+abs(x1-a[id1+1])+abs(x2-a[id1+1]);
//writeln('右:',abs(y1-y2)+abs(x1-a[id1+1])+abs(x2-a[id1+1]));
if abs(y1-y2)+abs(x1-a[id1-1])+abs(x2-a[id1-1])<min then min:=abs(y1-y2)+abs(x1-a[id1-1])+abs(x2-a[id1-1]);
//writeln('左',abs(y1-y2)+abs(x1-a[id1-1])+abs(x2-a[id1-1]));
exit(min);
end;
end;
begin
for wjmzbmr:=1 to m do
begin
readln(x1,y1,x2,y2);
//坑爹
if(y1=0)or(y2=0) then begin writeln('欧教坑爹啊!!'); terminate; end; //输入数据不合法
//同侧
//if(y1*y2>0) then begin writeln(abs(x1-x2)+abs(y1-y2)); continue; end;
If((y1<0) and (y2<0))or((y1>0)and(y2>0)) then begin writeln(abs(x1-x2)+abs(y1-y2));continue; end; //千万注意不要y1*y2,会爆Longint
//不同侧
//if(y1*y2<0) then writeln(solve);
if((y1>0) and (y2<0))or((y1<0)and(y2>0)) then writeln(solve); //不同侧处理
end;
end;
Begin
init;
main;
terminate;
End.