jzoj. 3838. 【NOIP2014模拟9.14】Super Big Stupid Cross

Description

“我是超级大沙茶”——Mato_No1
为了证明自己是一个超级大沙茶，Mato 神犇决定展示自己对叉（十字型）有多么的了解。
Mato 神犇有一个平面直角坐标系，上面有一些线段，保证这些线段至少与一条坐标轴平行。Mato 神犇需要指出，这些线段构成的最大的十字型有多大。
称一个图形为大小为R（R 为正整数）的十字型，当且仅当，这个图形具有一个中心点，它存在于某一条线段上，并且由该点向上下左右延伸出的长度为R 的线段都被已有的线段覆盖。
你可以假定：没有两条共线的线段具有公共点，没有重合的线段。

Input

第一行，一个正整数N，代表线段的数目。
以下N 行，每行四个整数x1,y1,x2,y2（x1=x2 或y1=y2），描述了一条线段。

Output

当不存在十字型时：输出一行“Human intelligence is really terrible”（不包括引号）。
否则：输出一行，一个整数，为最大的R。

Sample Input

输入1：
1
0 0 0 1
输入2：
3
-1 0 5 0
0 -1 0 1
2 -2 2 2

Sample Output

输出1：
Human intelligence is really terrible
输出2：
2

Data Constraint

对于50%的数据：N≤1000。
对于100%的数据：1≤N≤100000，所有坐标的范围在-10^9~10^9 中。
后50%内，所有数据均为随机生成。

背景：这题原题好像叫超流超导对撞机，好奇怪的名字。

分析：一开始我是暴力的，50分。然后我发现一个很惊讶的事实。先把长度大到小排序，如果当前长度的一半小于我们已经计算出最大的ans，就可以break掉。正解好像是用二分答案+set维护。

代码：

const
 maxn=100001;
type
 node=record
  x,y,u,v,l:longint;
 end;
arr=array [1..maxn] of node;

var
 a,b:arr;
 numa,numb,n,i,j,ans,x0,y0,x1,y1,p,l,r,d:longint;
 flag:boolean;

procedure qsort(l,r:longint;var b:arr);
  var
    i,j,key:longint;
    temp:node;
  begin
    if l>=r then exit;
    i:=l;j:=r;
    key:=b[l+random(r-l+1)].l;
    repeat
      while  (b[i].l>key) do inc(i);
      while  (b[j].l<key) do dec(j);
      if i<=j then
      begin
        temp:=b[i];b[i]:=b[j];b[j]:=temp;
        inc(i);dec(j);
      end;
    until i>j;
    qsort(l,j,b);
    qsort(i,r,b);
  end;

function min(x,y:longint):longint;
 begin
  if x<y then exit(x)
         else exit(y);
 end;

procedure swep(var x,y:longint);
 begin
  x:=x xor y;
  y:=x xor y;
  x:=x xor y;
 end;

begin
 readln(n);
 flag:=true;
 for i:=1 to n do
  begin
   readln(x0,y0,x1,y1);
   if x0=x1 then
    begin
     if y0>y1 then swep(y0,y1);
     inc(numb);
     b[numb].x:=x0;
     b[numb].y:=y0;
     b[numb].u:=x1;
     b[numb].v:=y1;
     b[numb].l:=y1-y0;
    end
   else
    begin
     if x0>x1 then swep(x0,x1);
     inc(numa);
     a[numa].x:=x0;
     a[numa].y:=y0;
     a[numa].u:=x1;
     a[numa].v:=y1;
     a[numa].l:=x1-x0;
    end;
  end;
 qsort(1,numb,b);
 qsort(1,numa,a);
 for i:=1 to numa do
  begin
   for j:=1 to numb do
    begin
     if (b[j].l+1) div 2<=ans then break;
     x0:=b[j].x;
     y0:=a[i].y;
     d:=min(min(x0-a[i].x,a[i].u-x0),min(y0-b[j].y,b[j].v-y0));
     if d>ans then ans:=d;
    end;
   if (a[i].l+1) div 2<=ans then break;
  end;
 if ans=0 then writeln('Human intelligence is really terrible')
          else writeln(ans);
end.