题目描述
经过11 年的韬光养晦,某国研发出了一种新的导弹拦截系统,凡是与它的距离不超过其工作半径的导弹都能够被它成功拦截。当工作半径为0 时,则能够拦截与它位置恰好相同的导弹。
但该导弹拦截系统也存在这样的缺陷:每套系统每天只能设定一次工作半径。而当天的使用代价,就是所有系统工作半径的平方和。
某天,雷达捕捉到敌国的导弹来袭。由于该系统尚处于试验阶段,所以只有两套系统投入工作。如果现在的要求是拦截所有的导弹,请计算这一天的最小使用代价。
输入
第一行包含4 个整数x1、y1、x2、y2,每两个整数之间用一个空格隔开,表示这两套导弹拦截系统的坐标分别为(x1, y1)、(x2, y2)。
第二行包含1 个整数N,表示有N 颗导弹。接下来N 行,每行两个整数x、y,中间用一个空格隔开,表示一颗导弹的坐标(x, y)。不同导弹的坐标可能相同。
输出
输出只有一行,包含一个整数,即当天的最小使用代价。
样例
输入:
样例一:
0 0 10 0
2
-3 3
10 0
样例二:
0 0 6 0
5
-4 -2
-2 3
4 0
6 -2
9 1
输出:
样例一:
18
样例二:
30
提示
样例解释 1
样例1 中要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分别为18 和0。
样例解释 2
样例中的导弹拦截系统和导弹所在的位置如下图所示。要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分别为20 和10。
数据范围与提示
对于10%的数据,N = 1
对于20%的数据,1 ≤ N ≤ 2
对于40%的数据,1 ≤ N ≤ 100
对于70%的数据,1 ≤ N ≤ 1000
对于100%的数据,1 ≤ N ≤ 100000,且所有坐标分量的绝对值都不超过1000。
参考代码
#include<bits/stdc++.h>
using namespace std;
int p1,q1,p2,q2;
int n, x, y;
struct Node{
int dis1;//导弹到一系统的距离平方
int dis2;//导弹到二系统的距离平方
};
Node a[100005];
bool cmp(Node x, Node y){
return x.dis1<y.dis1;
}
int main(){
cin >> p1 >> q1 >> p2 >> q2 >> n;
for(int i=1;i<=n;i++){
cin >> x >> y;
a[i].dis1=(x-p1)*(x-p1)+(y-q1)*(y-q1);
a[i].dis2=(x-p2)*(x-p2)+(y-q2)*(y-q2);
}
sort(a+1,a+n+1, cmp);//按照一系统的拦截的半径平方从小到大排序
int r=0; //二系统拦截的半径平方的最大值
int ans=INT_MAX; //答案一系统和二系统 截的半径平方和
for(int i=n;i>=0;i--){ //注意这里循环n+1次 (i=0所有导弹都被二系统拦截)
r=max(r, a[i+1].dis2);
ans=min(ans, a[i].dis1+r);
}
cout << ans;
return 0;
}