NOIP2010导弹拦截
题目描述
经过11 年的韬光养晦,某国研发出了一种新的导弹拦截系统,凡是与它的距离不超过其工作半径的导弹都能够被它成功拦截。当工作半径为0 时,则能够拦截与它位置恰好相同的导弹。但该导弹拦截系统也存在这样的缺陷:每套系统每天只能设定一次工作半径。而当天的使用代价,就是所有系统工作半径的平方和。
某天,雷达捕捉到敌国的导弹来袭。由于该系统尚处于试验阶段,所以只有两套系统投入工作。如果现在的要求是拦截所有的导弹,请计算这一天的最小使用代价。
「提示」
两个点(x1, y1)、(x2, y2)之间距离的平方是(x1− x2)2+(y1−y2)2。两套系统工作半径r1、r2 的平方和,是指r1、r2 分别取平方后再求和,即r12+r22。
「输入输出样例1」
missile.in
|
missile.out
|
0 0 10 0
2
-3 3
10 0
|
18
|
「样例 1 说明」
样例1 中要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分别为18 和0。
「输入输出样例2」
missile.in
|
missile.out
|
0 0 6 0
5
-4 -2
-2 3
4 0
6 -2
9 1
|
30
|
「样例2 说明」
样例中的导弹拦截系统和导弹所在的位置如下图所示。要拦截所有导弹,在满足最小使用代价的前提下,两套系统工作半径的平方分别为20 和10。
输入
输入文件名missile.in。
第一行包含4 个整数x1、y1、x2、y2,每两个整数之间用一个空格隔开,表示这两套导弹拦截系统的坐标分别为(x1, y1)、(x2, y2)。
第二行包含1 个整数N,表示有N 颗导弹。接下来N 行,每行两个整数x、y,中间用一个空格隔开,表示一颗导弹的坐标(x, y)。不同导弹的坐标可能相同。
输出
输出文件名missile.out。
输出只有一行,包含一个整数,即当天的最小使用代价。
提示
对于10%的数据,N = 1
对于20%的数据,1 ≤ N ≤ 2
对于40%的数据,1 ≤ N ≤ 100
对于70%的数据,1 ≤ N ≤ 1000
对于100%的数据,1 ≤ N ≤ 100000,且所有坐标分量的绝对值都不超过1000。
题解
设拦截系统为 a , b
按照导弹到其中一个拦截系统 a 的距离排序,将离 a 最近的 i 个导弹都交给 a ,其余给 b
倒序枚举断点,每次更新答案
需要用到结构体及其排序
结构体的定义
struct data{
//data 是结构体的名称
int x,y,s1,s2; // 定义元素
}
data a; //声明一个结构体a
data b[10]; //声明结构体组
每个结构体中有几个元素,这几个元素打包在一起,比如本题,一个导弹有其坐标和到两个系统的距离,我们按照一个元素排序时,整个结构体是一起移动位置的
b[0].x = 1; //赋值b[0]的元素 x 为1
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 |
#include<map> #include<set> #include<stack> #include<queue> #include<cmath> #include<cstdio> #include<vector> #include<cstring> #include<cstdlib> #include<iostream> #include<algorithm> #define inf 1000000000 #define ll long long using namespace std; int n; int x1,x2,y1,y2; int mn=inf; struct data{ //data 是结构体的名称 int x,y,s1,s2; }a[100005]; bool cmp(data a,data b) { return a.s1<b.s1; //结构体比较函数,可以理解为定义小于号 //即a.s1<b.s1时return 1 否则return 0 } int main() { scanf("%d%d%d%d",&x1,&y1,&x2,&y2); scanf("%d",&n); for(int i=1;i<=n;i++) { scanf("%d%d",&a[i].x,&a[i].y); a[i].s1=(a[i].x-x1)*(a[i].x-x1)+(a[i].y-y1)*(a[i].y-y1); a[i].s2=(a[i].x-x2)*(a[i].x-x2)+(a[i].y-y2)*(a[i].y-y2); //计算距离 } sort(a+1,a+n+1,cmp); int rb=0; for(int i=n;i;i--)//从离a最远的导弹开始枚举 { rb=max(a[i+1].s2,rb);//将i+1号导弹交给系统b,更新系统b的半径 mn=min(mn,a[i].s1+rb);//更新答案 } printf("%d",mn); return 0; } |
Subscribe