导弹拦截

前言
看前方……对面鲁班七号发射了一个炮弹，快跑…Oh，no！团灭。这个导弹真“捣蛋”。今天，嘟嘟老师给大家讲一道含金量也挺高，挺搞事情的一道题——导弹拦截。
题目描述
经过 11 年的韬光养晦，某国研发出了一种新的导弹拦截系统，凡是与它的距离不超过其工作半径的导弹都能够被它成功拦截。当工作半径为 0 时，则能够拦截与它位置恰好相同的导弹。但该导弹拦截系统也存在这样的缺陷：每套系统每天只能设定一次工作半径。而当天的使用代价，就是所有系统工作半径的平方和。
某天，雷达捕捉到敌国的导弹来袭。由于该系统尚处于试验阶段，所以只有两套系统投入工作。如果现在的要求是拦截所有的导弹，请计算这一天的最小使用代价。
输入
第一行包含 4 个整数x1、y1、x2、y2，每两个整数之间用一个空格隔开，表示这两套导弹拦截系统的坐标分别为(x1, y1)、(x2, y2)。
第二行包含 1 个整数 N，表示有 N颗导弹。接下来 N行，每行两个整数 x、y，中间用一个空格隔开，表示一颗导弹的坐标(x, y)。不同导弹的坐标可能相同。
输出
输出只有一行，包含一个整数，即当天的最小使用代价。
提示
两个点(x1, y1)、(x2, y2)之间距离的平方是(x1− x2)^2+(y1−y2)^2。
两套系统工作半径 r1、r2的平方和，是指 r1、r2 分别取平方后再求和，即 r1^2+r2^2。
样例输入1
0 0 10 0
2
-3 3
10 0
样例输出1
18
样例 1 说明
样例 1 中要拦截所有导弹，在满足最小使用代价的前提下，两套系统工作半径的平方分别为 18 和 0。
样例输入2
0 0 6 0
5
-4 -2
-2 3
4 0
6 -2
9 1
样例输出2
30
样例 2 说明
样例2中的导弹拦截系统和导弹所在的位置如下图所示。要拦截所有导弹，在满足最小使用代价的前提下，两套系统工作半径的平方分别为 20 和 10。

思路
首先嘟嘟老师想吐个槽，像这题一样，很多题的题目描述不清楚，导致看不懂题而爆蛋（0分）。最近的题目描述真难读懂啊。第一次读此题就能读懂的孩纸很从命（聪明），李姐（理解）能腻特别高。反正嘟嘟老师第一次没读懂此题。目前AC了，都对此题还有点疑问。
此题纯属搞事情。
首先用勾股定理算出每个点离系统的距离，用a【i，1】表示。a【i，2】来表示第几个导弹。这是第一个防导弹系统。第二个系统也是这样，用b【i，1】和b【i，2 】来表示。而还要用两个数组，因为我们后面要算距离和，用c【i】表示a【i，1】，用d【i】来表示b【i，1】。然后用快排排一下a数组和b数组。然后选出最大的a【i，1】和b【i，1】，取其中较小的，然后每个数都判断一下（除了a【n，1】和b【n，1】，因为在上一步已经判断过了）。算出最大中的最小的距离和，再跟之前的a【n，1】或b【n，1】来比较，输出那个较小值就行了。就是这么easy！！！
对面的童鞋看过来，看过来，看过来，下面的标程很精彩，给我的博客刷访问量。
……
代码

uses math;
var
        a,b:array[1..100000,1..2] of longint;
        c,d:array[1..100000] of longint;
        i,j,k,n,x,y,x1,x2,y1,y2,ans:longint;
procedure sb(l,r:longint);
var
        i,j,k,p:longint;
begin
        i:=l;
        j:=r;
        k:=a[(l+r) div 2,1];
        while (i<j) do
        begin
                while (a[i,1]<k) and (r>i) do inc(i);
                while (a[j,1]>k) and (l<j) do dec(j);
                if (i<=j) then
                begin
                        p:=a[i,1];
                        a[i,1]:=a[j,1];
                        a[j,1]:=p;
                        p:=a[i,2];
                        a[i,2]:=a[j,2];
                        a[j,2]:=p;
                        inc(i);
                        dec(j);
                end;
        end;
        if (i<r) then sb(i,r);
        if (l<j) then sb(l,j);
end;
procedure lsb(l,r:longint);
var
        i,j,k,p:longint;
begin
        i:=l;
        j:=r;
        k:=b[(l+r) div 2,1];
        while (i<j) do
        begin
                while (b[i,1]<k) and (r>i) do inc(i);
                while (b[j,1]>k) and (l<j) do dec(j);
                if (i<=j) then
                begin
                        p:=b[i,1];
                        b[i,1]:=b[j,1];
                        b[j,1]:=p;
                        p:=b[i,2];
                        b[i,2]:=b[j,2];
                        b[j,2]:=p;
                        inc(i);
                        dec(j);
                end;
        end;
        if (i<r) then lsb(i,r);
        if (l<j) then lsb(l,j);
end;
begin
        assign(input,'missile.in');reset(input);
        assign(output,'missile.out');rewrite(output);
        readln(x1,y1,x2,y2);
        readln(n);
        for i:=1 to n do
        begin
                read(x,y);
                a[i,1]:=sqr(abs(x1-x))+sqr(abs(y1-y));
                a[i,2]:=i;
                c[i]:=a[i,1];
                b[i,1]:=sqr(abs(x2-x))+sqr(abs(y2-y));
                b[i,2]:=i;
                d[i]:=b[i,1];
        end;
        sb(1,n);
        lsb(1,n);
        ans:=min(a[n,1],b[n,1]);
        j:=-1;
        k:=-1;
        for i:=n-1 downto 1 do
        begin
                j:=max(j,c[b[i+1,2]]);
                k:=max(k,d[a[i+1,2]]);
                ans:=min(ans,min(k+a[i,1],j+b[i,1]));
        end;
        writeln(ans);
        close(input);
        close(output);
end.