题意
- 给你平面上n个点,和2个圆心点,问你最短的两个半径平方和,使得所有点都能至少在其中一个圆内或圆上
思路
- 这题暴力就可解了,对一个圆试它的半径,距离可以到第i个点,这样会覆盖一些点,然后算另一个圆半径,即另一个圆要能覆盖到其它所有点
- 这样O(n ^ 2)即可解
- 但是单这么做实现会麻烦一些,我们可以先把点按照到第一个圆心点的距离排序,这样我们遍历到第i个点时可知,此时前i个点会被第一个圆包住,后n-i个点则需要被第二个圆包住,所以求圆心2到i+1~ n点距离的最大值即可得到r2
- 但是这样,我们可以想到,其实我们只要预处理出来这些排好序的点到圆心2的最大距离的后缀即可,这样我们就可以在O(nlog n)的时间内解决问题了
实现的是n ^ 2算法,nlogn的差不多
实现
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
const ll inf = 1e16;
const int maxn = 2005;
pii p[maxn];
pii pos[2];
int mark[maxn];
inline ll dist(pii& x,pii& y){
return (ll)(x.fi - y.fi) * (ll)(x.fi - y.fi) + (ll)(x.se - y.se) * (ll)(x.se - y.se);
}
bool cmp(pii x,pii y){
return dist(x,pos[0]) < dist(y,pos[0]);
}
int main(){
int n;
cin>>n>>pos[0].fi>>pos[0].se>>pos[1].fi>>pos[1].se;
for (int i=1;i<=n;i++){
scanf("%d%d",&p[i].fi,&p[i].se);
}
sort(p+1,p+n+1,cmp);
ll ans = inf;
p[0] = pos[0];
for (int i=0;i<=n;i++){
ll r1 = dist(pos[0],p[i]);
ll r2 = 0;
for (int j=i+1;j<=n;j++){
r2 = max(r2,dist(pos[1],p[j]));
}
ans = min(ans,r1+r2);
}
cout << ans << endl;
}