[SDOI2012] 最近最远点对

---

题目描述

给定平面直角坐标系上的N个点，分别求出距离最近的两个点的距离和距离最远的两个点的距离。注意，距离为直线距离。

---

输入格式

第一行一个整数，N。
接下来N行每行两个整数，xi,yi，表示第i个点的X坐标与Y坐标。

---

输出格式

总共一行，两个浮点数，为最短距离与最长距离，保留两位小数。

---

样例数据

样例输入

4
0.0 0.0
0.0 1.0
1.0 0.0
1.0 1.0

样例输出

1.00 1.41

---

数据范围

对于30%的数据，1<=N<=2000
对于70%的数据，1<=N<=20000
对于100%的数据，1<=N<100000，输入数据中所有数不超过10^9。

说明

本题采用自定义校验器。对于一个测验点，若只是答案1或答案2正确，即可得到一半的分数。但前提保证输出且只输出两个浮点数，或多或少都不给分。

---

题目分析

智障写法
最近点对（二分）+最远点对（旋转卡壳）

---

源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
inline const int Get_Int() {
    int num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
const double eps=1e-10;
int DoubleCompare(double x) { //精度三出口判断与0关系
    if(fabs(x)<eps)return 0; //=0
    else if(x<0)return -1; //<0
    else return 1; //>0
}
struct Point {
    double x,y;
};
bool cmp(Point a,Point b) {
    return a.x<b.x||(a.x==b.x&&a.y<b.y);
}
double Dist(Point a,Point b) {
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
struct Vector {
    double x,y;
};
Vector operator - (Point a,Point b) {
    Vector c;
    c.x=b.x-a.x;
    c.y=b.y-a.y;
    return c;
}
double Cross(Vector a,Vector b) { //叉积
    return a.x*b.y-b.x*a.y;
}
double Area(Point a,Point b,Point c) { //三点的平行四边形有向面积
    Vector u=b-a;
    Vector v=c-a;
    return Cross(u,v);
}
double Area(int n,Point* P) { //计算多边形有向面积（剖分法）
    double ans=0;
    for(int i=2; i<n; i++)ans+=Area(P[1],P[i],P[i+1]);
    return ans/2;
}
int ConvexHull(int n,Point* p,Point* ans) { //Andrew算法求凸包（求上链与下链）：p是点数组，ch是凸包顶点，返回顶点数
    //输入不能有重复点，若要凸包边上没有输入点，将两个<=改为<
    sort(p+1,p+n+1,cmp);
    int top=0;
    for(int i=1; i<=n; i++) {
        while(top>1&&DoubleCompare(Cross(ans[top-1]-ans[top-2],p[i]-ans[top-2]))<=0)top--;
        ans[top++]=p[i];
    }
    int k=top;
    for(int i=n-1; i>=1; i--) {
        while(top>k&&DoubleCompare(Cross(ans[top-1]-ans[top-2],p[i]-ans[top-2]))<=0)top--;
        ans[top++]=p[i];
    }
    if(n>1)top--;
    return top;
}
double Rotating_Calipers(int n,Point* p) { //旋转卡壳求最远点对距离
    int b=2; //旋转顶点b
    double ans=0;
    p[n+1]=p[1];
    for(int a=1; a<=n; a++) { //旋转顶点a
        while(abs(Area(p[a],p[a+1],p[b]))<abs(Area(p[a],p[a+1],p[b+1])))b=(b+1)%n;
        ans=max(ans,max(Dist(p[a],p[b]),Dist(p[a+1],p[b])));
    }
    return ans;
}
//
Point a[200005],b[200005];
bool cmpx(Point a,Point b) {
    return a.x<b.x;
}
bool cmpy(Point a,Point b) {
    return a.y>b.y;
}
double Find(int Left,int Right) {
    if(Left==Right)return 1e20; //只有一个点
    if(Left+1==Right)return (a[Left].x-a[Right].x)*(a[Left].x-a[Right].x)+(a[Left].y-a[Right].y)*(a[Left].y-a[Right].y);
    int mid=(Left+Right)/2;
    double pos=(a[mid].x+a[mid+1].x)/2;
    double minn=min(Find(Left,mid),Find(mid+1,Right));
    int cnt=0;
    Point tmp[Right-Left+5];
    for(int i=Left; i<=Right; i++)
        if(abs(a[i].x-pos)<minn) {
            cnt++;
            tmp[cnt].x=a[i].x;
            tmp[cnt].y=a[i].y;
        }
    sort(tmp+1,tmp+cnt+1,cmpy);
    for(int i=1; i<cnt; i++)
        for(int j=i+1; j<=i+8; j++) {
            if(j>cnt)break;
            minn=min(minn,double((tmp[i].x-tmp[j].x)*(tmp[i].x-tmp[j].x)+(tmp[i].y-tmp[j].y)*(tmp[i].y-tmp[j].y)));
        }
    return minn;
}
int n;
int main() {
    scanf("%d",&n);
    for(int i=1; i<=n; i++)scanf("%lf%lf",&a[i].x,&a[i].y);
    int cnt=ConvexHull(n,a,b);
    sort(a+1,a+n+1,cmpx);
    printf("%0.2lf %0.2lf\n",sqrt(Find(1,n)),Rotating_Calipers(n,b));
    return 0;
}

---