【题目描述】
经过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。
【输出】
只有一行,包含一个整数,即当天的最小使用代价。
【输入样例】
0 0 10 0
2
-3 3
10 0
【输出样例】
18
【提示】
【样例 1 说明】
样例1 中要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分
别为18 和0。
【输入输出样例2】
输入:
0 0 6 0
5
-4 -2
-2 3
4 0
6 -2
9 1
输出:
30
【样例2 说明】
样例中的导弹拦截系统和导弹所在的位置如下图所示。要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分别为20 和10。
【数据范围】
对于10%的数据,N = 1
对于20%的数据,1 ≤ N ≤ 2
对于40%的数据,1 ≤ N ≤ 100
对于70%的数据,1 ≤ N ≤ 1000
对于100%的数据,1 ≤ N ≤ 100000,且所有坐标分量的绝对值都不超过1000。
【解题思路】
考虑了每个导弹系统的拦截范围,并且通过枚举不同的拦截半径来找到能够拦截所有导弹的最小总代价。
具体来说:
-
计算导弹到导弹系统的距离:首先,计算每个导弹到两个导弹系统的距离。这个距离代表了导弹离导弹系统的远近,越远意味着需要更大的拦截半径才能拦截到导弹。
-
按照距离排序:将导弹按照到第一个导弹系统的距离进行排序,这样可以让距离较远的导弹排在前面,方便后续的枚举拦截半径。
-
枚举拦截半径:从最近的导弹开始,逐步增大拦截半径,尝试拦截所有导弹。在枚举过程中,记录每个拦截半径下的总代价,并选择最小的总代价。
-
输出结果:最终输出能够拦截所有导弹的最小总代价。
通过这种思路,能够遍历所有可能的拦截方案,并且选择最优的方案来拦截所有导弹,
【代码实现】
#include <cstdio>
#include <algorithm>
using namespace std;
int x1, y1, x2, y2, n, ans, r1, r2;
struct node {
int x, y, d1;
}mis[100005]; // 定义结构体表示导弹的横纵坐标和到第一个导弹系统的距离平方
int dis(int x, int y, int a, int b) {
return (x - a) * (x - a) + (y - b) * (y - b);
} // 计算两点之间的距离的平方
bool cmp(node a, node b) {
return a.d1 > b.d1;
} // 按照距离平方由大到小排序
int main() {
int x, y;
scanf("%d%d%d%d%d", &x1, &y1, &x2, &y2, &n);
// 读取第一个导弹系统和第二个导弹系统的坐标以及导弹的数量
for (int i = 1; i <= n; i++) {
scanf("%d%d", &x, &y);
// 读取每颗导弹的坐标
mis[i].x = x;
mis[i].y = y;
mis[i].d1 = dis(x1, y1, mis[i].x, mis[i].y);
// 计算每颗导弹到第一个导弹系统的距离平方
}
sort(mis + 1, mis + n + 1, cmp);
// 按照距离平方由大到小排序
r1 = mis[1].d1; // 初始化最小使用代价为最近一颗导弹到第一个系统的距离
ans = r1;
for (int i = 2; i <= n + 1; i++) {
// 一定要注意这里是n+1,因为可能一系统一个导弹都没拦截
r1 = mis[i].d1;
r2 = max(r2, dis(x2, y2, mis[i - 1].x, mis[i - 1].y));
// 更新第二个系统拦截半径的最大值
ans = min(ans, r1 + r2);
// 更新总代价的最小值
}
printf("%d", ans);
// 输出最小总代价
return 0;
}