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

原题链接:

平面最近点对（加强版） - 洛谷https://www.luogu.com.cn/problem/P1429

题目背景

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

题目描述

给定平面上 n 个点，找出其中的一对点的距离，使得在这 n 个点的所有点对中，该距离为所有点对中最小的

输入格式

第一行：n ，保证 2≤n≤200000 。

接下来 n 行：每行两个实数：x y ，表示一个点的行坐标和列坐标，中间用一个空格隔开。

输出格式

仅一行，一个实数，表示最短距离，精确到小数点后面 4 位。

输入输出样例

输入 #1复制

3
1 1
1 2
2 2

输出 #1复制

1.0000

说明/提示

数据保证 0≤x,y≤10^9 

Solution:

        玄学题目get√,鬼想得到解法.先把所有点按x(或y)坐标排序,然后所有的点x坐标近了距离更容易接近了。然后每个点都和后五个点的距离和ans取min,再将随机旋转几个角度。旋转之后每个点的相对位置其实是不会改变的,于是就可以再排序,再跟后五个点的距离和ans取min最后就怎么玄学AC了呗。好像也只能这样了

        一个点(x2,y2)依照(x1,y1)点旋转的公式为x = (x2-x1)*cos(角度)-(y2-y1)*sin(角度)+x1,y = (x2-x1)*sin(角度)+(y2-y1)*cos(角度)+y1

        别问我为啥。。看题解的

Code:

#include<bits/stdc++.h>
#define pi acos(-1.0)
using namespace std;
struct node{
	double x,y;
}a[200005];
int n;
bool cmp(node a,node b)
{
	return a.x<b.x;
}
double dis(node a,node b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double ans = 1e9;
void around(double x)
{
	x = x/180.0*pi;
	for(int i = 1;i<=n;++i)
	{
		double ddx = a[i].x,ddy = a[i].y;
		double xn,yn;
		xn = ddx*cos(x)-ddy*sin(x);
		yn = ddx*sin(x)+ddy*cos(x);
		a[i].x = xn,a[i].y = yn;
	}
	sort(a+1,a+1+n,cmp);
	for(int i = 1;i<=n;++i)
		for(int j = i+1;j<=n&&j<=i+6;++j)
			ans = min(ans,dis(a[i],a[j]));
}
int main()
{
	cin>>n;
	for(int i = 1;i<=n;++i)
		cin>>a[i].x>>a[i].y;
	around(0.0);
	around(90.0);
	around(45.0); 
	printf("%.4lf",ans);
	return 0;
}