P1429 平面最近点对（加强版）

求平面里最近两个点的距离，推荐两种做法。
洛谷 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θ
*/