经过11年的韬光养晦,某国研发出了一种新的导弹拦截系统,凡是与它的距离不超过其工作半径的导弹都能够被它成功拦截。当工作半径为0时,则能够拦截与它位置恰好相同的导弹。但该导弹拦截系统也存在这样的缺陷:每套系统每天只能设定一次工作半径。而当天的使用代价,就是所有系统工作半径的平方和。
某天,雷达捕捉到敌国的导弹来袭。由于该系统尚处于试验阶段,所以只有两套系统投入工作。如果现在的要求是拦截所有的导弹,请计算这一天的最小使用代价。
0 0 10 0
2
-3 3
10 0
18
两个点(x1,y1)、(x2,y2)之间距离的平方是(x1−x2)2+(y1−y2)2。
两套系统工作半径r1、r2的平方和,是指r1、r2分别取平方后再求和,即r12+r22。
对于100%的数据,1≤N≤100000,且所有坐标分量的绝对值都不超过1000。
思路:按照到a的距离排序,对b取最大可能到达的距离,在记录总的最小距离。
#include<queue> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define N 300010 #define maxn 100050 #define REP(i, a, b) for(int i = (a); i <= (b); ++ i) #define REP(j, a, b) for(int j = (a); j <= (b); ++ j) #define PER(i, a, b) for(int i = (a); i >= (b); -- i) using namespace std; template <class T> inline void rd(T &ret){ char c; ret = 0; while ((c = getchar()) < ‘0‘ || c > ‘9‘); while (c >= ‘0‘ && c <= ‘9‘){ ret = ret * 10 + (c - ‘0‘), c = getchar(); } } struct node{ int x,y,d,dis; }p[maxn]; int n,v,s,u,t,ans; bool cmp(node a,node b){return a.dis<b.dis;} int main(){ cin>>u>>v>>s>>t; rd(n); REP(i, 1, n){ cin>>p[i].x>>p[i].y; p[i].dis=(p[i].x-u)*(p[i].x-u)+(p[i].y-v)*(p[i].y-v); p[i].d=(p[i].x-s)*(p[i].x-s)+(p[i].y-t)*(p[i].y-t); } sort(p+1,p+1+n,cmp); int j=0,k=0x3f3f3f3f; PER(i, n, 0)j=max(j,p[i+1].d),k=min(k,j+p[i].dis); cout<<k<<endl; }
原文:https://www.cnblogs.com/czy-power/p/10369585.html