标签：ch int 分治 db mid 最短 rg 平面 include

#include<iostream>
#include<cstdio>
#include<iomanip>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>

#define ll long long
#define db double
#define rg register int

using namespace std;

const db inf=1e16;

struct su{
    db x,y;
}a[200005],b[200005];

db xx,yy;
int n;

inline int qr(){
    char ch; int sign=1;
    while((ch=getchar())<'0'||ch>'9')
        if(ch=='-')sign=-1;
    int res=ch^48;
    while((ch=getchar())>='0'&&ch<='9')
        res=res*10+(ch^48);
    return res*sign;
}

inline bool cmp_x(su x,su y){return x.x<y.x;}
inline bool cmp_y(su x,su y){return x.y<y.y;}

inline db dis2(su x,su y){
    xx=(y.x-x.x),yy=(y.y-x.y);
    return sqrt(xx*xx+yy*yy);
}

inline db dis3(su x,su y,su z){
    return dis2(x,y)+dis2(x,z)+dis2(y,z);
}

inline db find(int l,int r){
    if(l+1>=r)return inf;
    int mid=(l+r)>>1;
    db d=min(find(l,mid),find(mid+1,r));
    while(a[l].x+d<a[mid].x)++l;
    while(a[r].x-d>a[mid].x)--r;
    int t=0;
    for(rg i=l;i<=r;++i)b[++t]=a[i];
    sort(b+1,b+t+1,cmp_y);
    for(rg i=1;i<=t;++i)
        for(rg j=i+1;j<=t;++j)
            if(b[j].y-b[i].y>=d)break;
            else for(rg k=j+1;k<=t;++k)
                     if(b[k].y-b[i].y>=d)break;
                     else d=min(d,dis3(b[i],b[j],b[k])/2);
    return d;
}

int main(){
    freopen("math.in","r",stdin);
    freopen("math.out","w",stdout);
    n=qr();
    for(rg i=1;i<=n;++i)
        a[i].x=qr(),a[i].y=qr();
    sort(a+1,a+n+1,cmp_x);
    printf("%.6lf\n",find(1,n)*2);
    return 0;
}

标签：ch,int,分治,db,mid,最短,rg,平面,include
来源： https://www.cnblogs.com/812-xiao-wen/p/10363129.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。