最近點對，采用分治方法。過程：
1對原數組依據x左標從小到大排序。
2二分數組，左邊求出最小值，右邊求出最小值，我們求最小的。
3找出對于左右兩邊的可能小于當前最小值的最近點對，更新最小值。
這題目需要區分一下點，讓我們求的是閃兵到任意一個核電站的最短距離，加一個標志就可以了。
代碼：

#include <iostream> #include <stdio.h> #include <iomanip> #include <algorithm> #include <cmath> using namespace std; const int N=200000+1; const double INF=1e100; struct Point{ double x,y; bool flag; }; Point m[N+1]; int tmp[N+1]; double dis(Point p1,Point p2) { return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)); } double minL(double l1,double l2) { return l1<=l2?l1:l2; } bool cmpy(int a,int b) { return m[a].y<m[b].y; } bool cmpx(Point a,Point b) { return a.x<b.x; } double getMinLen(Point *s,int left,int right) { double rs=INF; if(left == right) return rs; if(left+1 == right) { if(s[left].flag == s[right].flag) return rs; return dis(s[left],s[right]); } int mid=(left+right)>>1; rs=getMinLen(s,left,mid); rs=minL(rs,getMinLen(s,mid+1,right));//rs現在是兩邊中最點值了 int i,j,num=0; for(i=left;i<=right;++i) { if(fabs(s[i].x-s[mid].x)<=rs)//過濾掉不可能是最近點對的點 tmp[num++]=i; } sort(tmp,tmp+num,cmpy);//按y坐標排序，減少比較次數 double d=INF; for(i=0;i<num;++i) { for(j=i+1;j<num;++j) { if(fabs(s[tmp[i]].y-s[tmp[j]].y)>=rs)//i，j兩點的y坐標已經大于rs,不需要在遍歷了 break; if(s[tmp[i]].flag!=s[tmp[j]].flag && (d=dis(s[tmp[i]],s[tmp[j]]))<rs) rs=d; } } // cout<<left<<" "<<right<<" "<<rs<<endl; return rs; } int main() { int t,n,i; scanf("%d",&t); while(t--) { scanf("%d",&n); for(i=0;i<n;++i) { scanf("%lf%lf",&m[i].x,&m[i].y); m[i].flag=0; } for(i=0;i<n;++i) { scanf("%lf%lf",&m[i+n].x,&m[i+n].y); m[i+n].flag=1; } n<<=1; sort(m,m+n,cmpx); double ans=getMinLen(m,0,n-1); cout<<setiosflags(ios::fixed)<<setprecision(3)<<ans<<endl; } return 0; }

總結

以上是生活随笔為你收集整理的poj3714 最近点对的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。