求平面里最近两个点的距离,推荐两种做法。
洛谷 P1429
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+100,inf=1e9;
int temp[N];
struct point{
double x,y;
}S[N];
bool cmp1(point a,point b){
if(a.x==b.x)return a.y<b.y;
return a.x<b.x;
}
bool cmp2(int a,int b){
return S[a].y<S[b].y;
}
double dist(int i,int j){
return sqrt((S[i].x-S[j].x)*(S[i].x-S[j].x)+(S[i].y-S[j].y)*(S[i].y-S[j].y));
}
double merge(int l,int r){
if(l==r) return inf;
if(l+1==r) return dist(l,r);
int mid=(l+r)/2;
double d=min(merge(l,mid),merge(mid+1,r));
int k=0;
for(int i=l;i<=r;i++)
if(fabs(S[i].x-S[mid].x)<d)
temp[k++]=i;
sort(temp,temp+k,cmp2);
for(int i=0;i<k;i++)
for(int j=i+1;j<k&&S[temp[j]].y-S[temp[i]].y<d;j++)
//写这里的时候忘记for里面的&&不满足会直接跳出循环,调了一段时间
d=min(d,dist(temp[i],temp[j]));
return d;
}
int main(){
int n;
scanf("%d",&n);
for(int i=0;i<n;i++)
scanf("%lf%lf",&S[i].x,&S[i].y);
sort(S,S+n,cmp1);
printf("%.4lf\n",merge(0,n-1));
}
- 洛谷旋转坐标题解
通过旋转坐标x,来得到_x,通过_x排序,最近点大概率在排序好的后面,不得已的时候大家可以用这种方法水一水。
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+100;
const double inf=1e7+0.1;
struct point{
double x,y,_x;
}S[N];
bool cmp(point a,point b){
return a._x<b._x;
}
int main(){
int n;
scanf("%d",&n);
double d,ans=inf,s=sin(1),c=cos(1);
for(int i=0;i<n;i++){
scanf("%lf%lf",&S[i].x,&S[i].y);
S[i]._x=S[i].x*c-S[i].y*s;
}
sort(S,S+n,cmp);
for(int i=0;i<n;i++)
for(int j=i+1;j<n&&j<i+5;j++){
d=sqrt((S[i].x-S[j].x)*(S[i].x-S[j].x)+(S[i].y-S[j].y)*(S[i].y-S[j].y));
ans=min(ans,d);
}
printf("%.4lf\n",ans);
}
/*旋转方式
x'=xcosθ-ysinθ
y'=xsinθ+ycosθ
*/