最近点对

    #include<iostream>  
    #include<cstdio>  
    #include<string>  
    #include<cstring>  
    #include<vector>  
    #include<cmath>  
    #include<queue>  
    #include<stack>  
    #include<map>  
    #include<set>  
    #include<algorithm>  
    using namespace std;  
    const int maxn=100010;  
    const double INF=10000000000000.0;  
    const int eps=1e-6;  
    int N;  
    double mind;  
    struct node  
    {  
        double x,y;  
        bool operator<(const node &a)const  
        {  
            if(fabs(x-a.x)<eps)return y<a.y;  
            return x<a.x;  
        }  
    }p[maxn],tmp[maxn];  
    bool cmp(node a,node b)  
    {  
        return a.y<b.y;  
    }  
    double cal(node a,node b)  
    {  
        double x=a.x-b.x;  
        double y=a.y-b.y;  
        return sqrt(x*x+y*y);  
    }  
    double solve(int l,int r)  
    {  
        if(l==r)return INF;  
        if(l+1==r)return cal(p[l],p[r]);  
        int mid=(l+r)>>1;  
        double d=INF;  
        d=min(d,solve(l,mid-1));  
        d=min(d,solve(mid+1,r));  
        int k=0;  
        for(int i=l;i<=r;i++)  
            if(fabs(p[mid].x-p[i].x)<=d)tmp[k++]=p[i];  
        sort(tmp,tmp+k,cmp);  
        for(int i=0;i<k;i++)  
        {  
            for(int j=i+1;j<k&&fabs(tmp[j].y-tmp[i].y)<=d;j++)  
            {  
                double dis=cal(tmp[i],tmp[j]);  
                if(dis<d)  
                    d=dis;  
            }  
        }  
        return d;  
    }  
    int main()  
    {//freopen("in.txt","r",stdin);  
        while(scanf("%d",&N)!=EOF)  
        {  
            for(int i=1;i<=N;i++)scanf("%lf%lf",&p[i].x,&p[i].y);  
            sort(p+1,p+1+N);  
            double ans=solve(1,N);  
            printf("%.2lf\n",ans/2);  
        }  
        return 0;  
    }  

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值