导弹拦截--upc寒假训练赛

题目

题目背景

经过11年的韬光养晦，某国研发出了一种新的导弹拦截系统，凡是与它的距离不超过其工作半径的导弹都能够被它成功拦截。当工作半径为0时，则能够拦截与它位置恰好相同的导弹。但该导弹拦截系统也存在这样的缺陷：每套系统每天只能设定一次工作半径。而当天的使用代价，就是所有系统工作半径的平方和。

某天，雷达捕捉到敌国的导弹来袭。由于该系统尚处于试验阶段，所以只有两套系统投入工作。如果现在的要求是拦截所有的导弹，请计算这一天的最小使用代价。

输入

第一行包含4个整数x1、y1、x2、y2，每两个整数之间用一个空格隔开，表示这两套导弹拦截系统的坐标分别为(x1,y1)、(x2,y2)。
第二行包含1个整数N(1≤N≤100000)，表示有N颗导弹。接下来N行，每行两个整数x、y，中间用一个空格隔开，表示一颗导弹的坐标(x,y)。不同导弹的坐标可能相同。

输出

输出只有一行，包含一个整数，即当天的最小使用代价。

样例输入

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。 

---

思路分析

错误思路：这个题我一开始想的是这个导弹由离他最近的拦截系统来拦截，但是发现结果是WA。哎，还是太单纯了。如果按照我的这个想法，后来想到一个很好解释的例子，在一个数轴上，0和5是拦截系统的位置，现在如果在-5和4处有一个导弹，很显然我们之前的思路就不对了，这个就只需要用在0位置的拦截系统就行。

正确思路：正确的做法是按照到1的距离进行排序，然后从找到一个合适的分界点 i 使得一部分由第一个拦截系统拦截，另一部分由第二个拦截系统拦截。这样就可以找出最小的平方和

废话不多说，上代码 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=1e5+10;
ll xa,xb,ya,yb;
struct node{
	ll x,y;
	ll disa,disb;
}a[N];
ll da,db;
ll n;
bool cmp(node c,node d){
	return c.disa>d.disa;
}
int main()
{
	cin>>xa>>ya>>xb>>yb;//输入两套拦截系统所在的位置坐标 
	cin>>n;//输入导弹的个数 
	for(ll i=1;i<=n;i++){
		cin>>a[i].x>>a[i].y;//输入第i个导弹的位置坐标 
		a[i].disa=(a[i].x-xa)*(a[i].x-xa)+(a[i].y-ya)*(a[i].y-ya);//计算到a的距离的平方 
		a[i].disb=(a[i].x-xb)*(a[i].x-xb)+(a[i].y-yb)*(a[i].y-yb);//计算到b的距离的平方 
	}
	sort(a+1,a+1+n,cmp);//根据距离第一个拦截系统的位置进行排序 
	ll ans=INT_MAX;//定义一个变量来储存答案，赋初值为无穷 
	ll r2=0;//定义一个变量来记录第二个拦截系统设定的半径 
	for(ll i=1;i<=n+1;i++){
		ll r1=a[i].disa;//第一个拦截系统的半径 
		r2=max(r2,a[i-1].disb);//求出第二个拦截系统的拦截半径 
		ans=min(r1+r2,ans);//求出最合适的分解点是的最小距离的平方和 
	}
	cout<<ans;
	return 0;
}

 操作过程举例