P1257 平面上的最接近点对&P1429 平面最近点对（加强版）

题目

题目&题目

思路

利用二分思想，划分子问题，显然这一步非常简单。
本题难点在于合并，显然左右2边的点可以合并，但是我们还需要在左右各选一个点合并，怎么办？
暴力
n方过百万，dp得零蛋
以上纯属搞笑
这里我们已经有了一个候选答案d（来自2个子问题），我们离选择的中心点最多不能超过d，显然如果超过一定不是更优的，使用这条性质可以筛走一部分数。同时我们发现，若2点的纵坐标>=d，也不会更优，这样再暴力两两匹配即可。
code：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<map> 
using namespace std;
int n;
struct f{
	double x,y;
} a[100001],b[100001];
bool cmp2(f a,f b)
{
	if (a.x==b.x) return a.y<b.y;
	return a.x<b.x;
}
bool cmp(f a,f b)
{
	if (a.y==b.y) return a.x<b.x;
	return a.y<b.y;
}
double dis(f a,f b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double ff(int l,int r)
{
	if (l==r) return 0x7fffffff;
	if (l+1==r)
	{
		return dis(a[l],a[r]);
	}
	double d1=ff(l,(l+r)>>1),d2=ff(((l+r)>>1)+1,r);
	d1=min(d1,d2);
	int tot=0,mid=(l+r)>>1;
	for (int i=l;i<=r;i++) if (fabs(a[i].x-a[mid].x)<=d1) b[tot++]=a[i];
	sort(b,b+tot,cmp);
	for (int i=0;i<tot;i++)
	{
		for (int j=i+1;j<tot;j++)
		{
			if (b[j].y-b[i].y<d1)
			{
				d1=min(d1,dis(b[i],b[j]));
			}
		}
	}
	return d1;
}
int main()
{
	double x,y;
	scanf("%lf%lf",&x,&y);
	while (x!=0||y!=0)
	{
		yy--;
		if (yy<=0) break;
		a[++n].x=x,a[n].y=y;
		scanf("%lf%lf",&x,&y); 
	}
	a[++n].x=x,a[n].y=y;
	sort(a+1,a+n+1,cmp2);
	for (int i=1;i<=n;i++)
	{
		if (i!=1&&a[i].x==a[i-1].x&&a[i].y==a[i-1].y)
		{
			cout<<"0.0000";
			return 0;
		}
	}
	printf("%.4f",ff(1,n));
    return 0;
}

同时放一下学校OJ的代码：
code：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<map> 
using namespace std;
int n;
struct f{
	double x,y;
} a[100001],b[100001];
bool cmp2(f a,f b)
{
	if (a.x==b.x) return a.y<b.y;
	return a.x<b.x;
}
bool cmp(f a,f b)
{
	if (a.y==b.y) return a.x<b.x;
	return a.y<b.y;
}
double dis(f a,f b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double ff(int l,int r)
{
	if (l==r) return 0x7fffffff;
	if (l+1==r)
	{
		return dis(a[l],a[r]);
	}
	double d1=ff(l,(l+r)>>1),d2=ff(((l+r)>>1)+1,r);
	d1=min(d1,d2);
	int tot=0,mid=(l+r)>>1;
	for (int i=l;i<=r;i++) if (fabs(a[i].x-a[mid].x)<=d1) b[tot++]=a[i];
	sort(b,b+tot,cmp);
	for (int i=0;i<tot;i++)
	{
		for (int j=i+1;j<tot;j++)
		{
			if (b[j].y-b[i].y<d1)
			{
				d1=min(d1,dis(b[i],b[j]));
			}
		}
	}
	return d1;
}
int main()
{
	double x,y;
	scanf("%lf%lf",&x,&y);
	while (x!=0||y!=0)
	{
		a[++n].x=x,a[n].y=y;
		scanf("%lf%lf",&x,&y); 
	}
	sort(a+1,a+n+1,cmp2);
	for (int i=1;i<=n;i++)
	{
		if (i!=1&&a[i].x==a[i-1].x&&a[i].y==a[i-1].y)
		{
			cout<<"0.0000";
			return 0;
		}
	}
	printf("%.4f",ff(1,n));
    return 0;
}