NOIP2010 普级组 导弹拦截

Description

　　经过11年的韬光养晦，某国研发出了一种新的导弹拦截系统，凡是与它的距离不超过其工作半径的导弹都能够被它成功拦截。当工作半径为0时，则能够拦截与它位置恰好相同的导弹。但该导弹拦截系统也存在这样的缺陷：每套系统每天只能设定一次工作半径。而当天的使用代价，就是所有系统工作半径的平方和。 
　　某天，雷达捕捉到敌国的导弹来袭。由于该系统尚处于试验阶段，所以只有两套系统投入工作。如果现在的要求是拦截所有的导弹，请计算这一天的最小使用代价。

Input

输入文件名missile.in。 
第一行包含4 个整数x1、y1、x2、y2，每两个整数之间用一个空格隔开，表示这两套导弹拦截系统的坐标分别为(x1, y1)、(x2,y2)。 
第二行包含1 个整数N，表示有N 颗导弹。接下来N 行，每行两个整数x、y，中间用一个空格隔开，表示一颗导弹的坐标(x,y)。不同导弹的坐标可能相同。

Output

输出文件名missile.out。 
输出只有一行，包含一个整数，即当天的最小使用代价。 
两个点(x1, y1)、(x2, y2)之间距离的平方是(x1−x2)^2+(y1−y2)^2。 
两套系统工作半径r1、r2 的平方和，是指r1、r2 分别取平方后再求和，即r1^2+r2^2。

Sample Input

 

样例10 0 10 02-3 310 0样例20 0 6 05-4 -2-2 34 06 -29 1

 

Sample Output

 

样例118【样例 1 说明】样例1 中要拦截所有导弹，在满足最小使用代价的前提下，两套系统工作半径的平方分别为18 和0。样例230【样例2 说明】样例中的导弹拦截系统和导弹所在的位置如下图所示。要拦截所有导弹，在满足最小使用代价的前提下，两套系统工作半径的平方分别为20 和10。

 

Hint

【数据范围】 
对于10%的数据，N = 1 
对于20%的数据，1 ≤ N ≤ 2 
对于40%的数据，1 ≤ N ≤ 100 
对于70%的数据，1 ≤ N ≤ 1000 
对于100%的数据，1 ≤ N ≤ 100000，且所有坐标分量的绝对值都不超过1000。

解题思路：先读入数据，先求出他们与两套系统距离的平方，储存在数组里，然后用快排对他们根据其中一个数组对两个数组进行排序，最后用循环求出当天的代价输出即可。

程序：

var

  a,b:array[0..100001] of longint;

  x1,y1,x2,y2,n:longint;

procedure init;

  var

    i,x,y:longint;

  begin

    readln(x1,y1,x2,y2);

    readln(n);

    for i:=1 to n do

     begin

       readln(x,y);

       a[i]:=sqr(x1-x)+sqr(y1-y);

       b[i]:=sqr(x2-x)+sqr(y2-y);

     end;

end;

procedure qsort(l,r:longint);

  var

    i,j,k,t:longint;

  begin

    if l>=r then exit;

    i:=l; j:=r;

    k:=a[(i+j) div 2];

    repeat

      whilea[i]>k do inc(i);

      whilea[j]

      ifi<=j then

       begin

        t:=a[i];a[i]:=a[j];a[j]:=t;

        t:=b[i];b[i]:=b[j];b[j]:=t;

         inc(i);dec(j);

       end;

    until i>j;

    qsort(i,r);

    qsort(l,j);

end;

procedure main;

  var

    i,max,ans:longint;

  begin

    init;

    qsort(1,n);

    max:=0;

    ans:=maxlongint;

    for i:=1 to n+1 do

     begin

       if a[i]+max

       if b[i]>max then max:=b[i];

     end;

    writeln(ans);

end;

begin

  main;

end.

版权属于: Chris

原文地址: http://blog.sina.com.cn/s/blog_83ac6af80102v66q.html

转载时必须以链接形式注明原始出处及本声明。