最近点对问题
hdu 1007的AC代码,对于于《编程之美》一书,HDU上的测试集有漏洞,按X+Y排序也可AC,此代码使用使用动态规划算法,使用C++超时,G++ AC,原因是sort函数的算法实现不同。
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
double MAX=1e10; //定义的最大距离,以在只有一个点的时返回无穷大
struct Node{
double x,y;
};
Node ar[100005],br[50005];
inline bool cmpx(Node a,Node b){return a.x<b.x;} //x坐标升序
inline bool cmpy(Node a,Node b){return a.y<b.y;} //y坐标升序
inline double min(double a,double b){return a<b?a:b;} //返回最小值
inline double dis(Node a,Node b){return sqrt(pow(a.x-b.x,2)+pow(a.y-b.y,2));} //返回点与点之间的距离
double cal(int s,int e) //s、e用来表示当前处理的数组中的下标位置
{
int mid,i,j,tail=0; //mid表示数组中间的位置下标 ,tail作为计数变量,是用来br数组存储标号的
double d; //d表示点对之间的距离
if(s==e) return MAX; //如果只有一个点
mid=(s+e)/2;
d=min(cal(s,mid),cal(mid+1,e)); //递归求出左右两边的最小距离
for(i=mid; i>=s&&ar[mid].x-ar[i].x<d; i--){
br[tail++]=ar[i];
}
for(i=mid+1;i<e&&ar[i].x-ar[mid].x<d;i++) //同上,筛选右边的点
{
br[tail++]=ar[i];
}
sort(br,br+tail,cmpy); //对矩形内的点按照y坐标进行排序
for(i=0;i<tail;i++){
for(j=i+1;j<tail&&j<i+8&&br[j].y-br[i].y<d;j++){
double td = dis(br[i],br[j]);
if(d>td){ //更新点的值
d=td;
}
}
}
return d; //返回最小的点对之间的距离
}
int main()
{
int n;
while(cin>>n&&n){
for(int i=0;i<n;i++){
scanf("%lf %lf",&ar[i].x,&ar[i].y);
}
sort(ar,ar+n,cmpx); //按x对ar排序
double d=cal(0,n);
printf("%.2lf\n",d/2.0);
}
return 0;
}