NOIP2010 普级组 导弹拦截

该博客介绍了一个数学问题,涉及导弹拦截系统的工作原理。系统每次设定工作半径的代价是半径平方和,目标是拦截所有导弹,求最小代价。通过读入数据、计算与系统的距离平方、排序和循环计算,可以找到最小代价并输出。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

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

 

样例1
0 0 10 0
2
-3 3
10 0

样例2
0 0 6 0
5
-4 -2
-2 3
4 0
6 -2
9 1

 

Sample Output

 

样例1
18

【样例 1 说明】
样例1 中要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分
别为18 和0。

样例2
30

【样例2 说明】
样例中的导弹拦截系统和导弹所在的位置如下图所示。要拦截所有导弹,在满足最小使
用代价的前提下,两套系统工作半径的平方分别为20 和10。
NOIP2010 <wbr>普级组 <wbr>导弹拦截

 

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

      while a[i]>k do inc(i);

      while a[j]

      if i<=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

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值