1、算法:
1) 总体的思路:
A、首先根据快排算法将数组按照x轴对数组a进行排序;
B、总体采用递归的表现形式,处理方法为conquerPair,当点的数目为1时返回无穷大,为2时,直接返回两点的距离;
C、对start到mid范围的点使用conquerPair方法得到l_min, 对mid+1到end范围的点使用conquerPair方法得到r_min
D、取d_min=Min(l_min,r_min);
E、在x[mid]-d_min到x[mid]+d_min进行操作findMidMin方法操作,找出横跨left区和right区的最小点,与d_min比较大小,得到的最小值即为最小值;
2) findMidMin方法的处理:
A、遍历在x[mid]-d_min范围的每一个点i;
B、对于每一个点i,计算它与点j的距离,j的限制条件是:x在x[mid]到x[mid]+d_min,y在y[i]+d_min到y[i]-d_min;
C、同时,满足要求的个数k<6,这一点后面证明;
D、将计算结果与d_min比较,返回小的一个即可;
3)为什么在限定范围中只需要找小于6个点,具体查看:
由下图知,假设在图中范围存在6个点满足小于d_min的点,分别位于不同的矩形中,在一个1/6的矩形中,最长为(5/6)*d_min,即一个矩形只能有一个点,再增加一个点就会与d_min的大小产生矛盾。
2、关键代码:
1)代码1,对应1中1)“总体的思路”:
private static double conquerPair(double[] x, double[] y, int start, int end) {
int mid=(start+end)/2,i=start,j=end;
int num=end-start+1;
if(num==1)
return Integer.MAX_VALUE;
else if(num==2)
return getDistance(x[i],y[i],x[j],y[j]);
double l_min=conquerPair(x,y,start,mid);
double r_min=conquerPair(x,y,mid+1,end);
double d_min=l_min<r_min?l_min:r_min;
return findMidMin(x,y,mid,d_min,start,end);
}
2)代码2,对应1的2)中“findMidMin方法的处理”:
private static double findMidMin(double[] x, double[] y, int mid, double d_min,
int start, int end) {
double min=d_min,l_border=x[mid]-d_min,r_border=x[mid]+d_min,t_border,b_border,temp;
int k;
for(int i=mid;i>=start&&x[i]>l_border;i--){
t_border=y[i]+d_min;
b_border=y[i]-d_min;
k=0;//满足要求的6个点的计数
for(int j=mid+1;j<=end&&k<6&&x[j]<r_border;j++){
if(y[j]<t_border&&y[j]>b_border){
k++;
temp=getDistance(x[i], y[i], x[j], y[j]);
if(min>temp)
min=temp;
}
}
}
return min;
}