「LuoguP1429」 平面最近点对（加强版）

题目描述

给定平面上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

题解

考场清晰的记得以前听过，并且记错做法还觉得自己是天才......

考场思路：按横轴排序，然后从左往右枚举每个点，把和该点的横坐标之差小于$ans$的点暴力算。

然后怕构造数据被卡（很多个点的横坐标差不多，竖轴差很多的话可以卡成$N$方），于是考虑用权值线段树优化竖轴，也就是从权值线段树中取竖轴在当前点加减ans范围的点暴力算。（考试$x,y$范围1e6）

这样应该就不会T了......

然后毒瘤评测人只给64M啊！哪里够开权值线段树啊！

于是把扫描线转了45°（也可以理解为把纸转45°），再做类暴力，这样就不好卡了。

然后就过了2333

 /*
     qwerta
     P1429 平面最近点对（加强版）
     Accepted
     100
     代码 C++，1.16KB
     提交时间 2018-10-19 15:55:35
     耗时/内存
     348ms, 2064KB
 */
 #include<algorithm>
 #include<iostream>
 #include<cstdio>
 #include<cmath>
 using namespace std;
 #define R register
 inline int read()
 {
     char ch=getchar();
     int x=;
     while(!isdigit(ch))ch=getchar();
     while(isdigit(ch)){x=x*+ch-'';ch=getchar();}
     return x;
 }
 const int MAXN=+,MAXX=+;
 struct emm{
     int x,y;
 }a[MAXN];
 const double k=-0.5;
 bool cmp(emm qaq,emm qwq)
 {
     return qaq.y-k*qaq.x<qwq.y-k*qwq.x;
 }
 inline double dis(int u,int v)
 {
     return sqrt(1LL*(a[u].x-a[v].x)*(a[u].x-a[v].x)+1LL*(a[u].y-a[v].y)*(a[u].y-a[v].y));
 }
 double gen2=sqrt();
 inline double je(int u,int v)
 {
     return abs(((a[u].y-k*a[u].x)-(a[v].y-k*a[v].x)))/gen2;
 }
 int main()
 {
     //freopen("dark.in","r",stdin);
     //freopen("dark.out","w",stdout);
     int n=read();
     for(R int i=;i<=n;++i)
       a[i].x=read(),a[i].y=read();
     sort(a+,a+n+,cmp);
     int l=,r=;
     double ans=dis(,);
     for(R int u=;u<=n;++u)
     {
         while(je(u,l)>ans)++l;
         while(je(r+,u)<ans&&r<n)++r;
         //cout<<u<<" "<<l<<" "<<r<<endl;
         for(R int v=l;v<=r;++v)
         if(u!=v)
         ans=min(ans,dis(u,v));
     }
     printf("%.4f",ans);
     return ;
 }