这个东西比较强。。。。。
当然,可以交智商税(开发KD—ree)
然后我并不会。。
于是乎只能看题解。。。。
就是把平面一直分割
然后每次根据左右得出的结论进行剪枝。。。
然后就可以优化过掉这题
大概是 O ( N l o g N ) O(NlogN) O(NlogN)上界挺松的
#include<bits/stdc++.h>
#define MAXN 3000001
using namespace std;
int n,q[MAXN],js;
struct node{
double x,y;
}t[MAXN];
bool cmp(node a , node b){
if(a.x != b.x)return a.x < b.x;
else return a.y < b.y;
}
bool cmp2(const int &a , const int &b){
return t[a].y < t[a].y;
}
double dis(int a , int b){
return sqrt((t[a].x - t[b].x) * (t[a].x - t[b].x) +
(t[a].y - t[b].y) * (t[a].y - t[b].y));
}
void init(){
cin>>n;
for(int i = 1 ; i <= n ; i++)cin>>t[i].x>>t[i].y;
sort(t + 1 , t + 1 + n , cmp);
}
double cal(int l , int r){
if(l == r)return 1 << 20;
if(l + 1 == r)return dis(l , r);
int mid = (l + r) >> 1;
double zz = min(cal(l , mid) , cal(mid + 1 , r));
js = 0;
for(int i = l ; i <= r ; i++)
if(fabs(t[mid].x - t[i].x) <= zz)
q[++js] = i;
sort(q + 1 , q + 1 + js , cmp2);
for(int i = 1 ; i <= js ; i++){
for(int j = i + 1 ; j <= js && t[q[j]].y - t[q[i]].y < zz ; j++){
zz = min(zz , dis(q[i] , q[j]));
}
}
return zz;
}
int main(){
init();
double jj = cal(1 , n);
printf("%.4lf\n" , jj);
}
启发的话。。。
想到暴力怎么做。。有些眉目后,要大胆去做。。。。
一些很像nlogn的数据,大多数都是分治结构。。。
考虑优化暴力。。