quoit Design

使用二分法进行加速

#include<iostream>
#include<cmath>
#include<algorithm>
#include<iomanip>
using namespace std;
const int mm=100008;
const double oo= 1e18;
class Point
{
public:
    double x,y;
}poin[mm];
int m;
int tmp[mm];
bool cmpxy(Point a,Point b)
{
    if(a.x==b.x)
        return a.y<b.y;
    else
        return a.x<b.x;
}
bool cmpy(int a,int b)
{
    return poin[a].y<poin[b].y;
}
double dist(Point a,Point b)
{
    double dis = (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
    return sqrt(dis);
}
double min_dist(int left,int right)
{
    double d=oo;
    if(left == right)                                                           //left=right，只有一个点，没有意义
        return d;
    else if(left+1==right)                                                      //这只有两个点
    return dist(poin[left],poin[right]);
    int mid = (left+right)/2;                                                   //取中间点，在主函数中已经进行了排序
    d= min(min_dist(left,mid),min_dist(mid+1,right));                           //递归，取出两边的最小的距离
    int pos=0;
    for(int i=left;i<=right;i++)                                                //判断在mid之间的2d范围内是否还有更小的d的值
    {
        if(fabs(poin[i].x-poin[mid].x)<=d)                                      //这里可以看出，如果双x的距离小于d，，才有可能使得两点的距离小于d。这里使用x进行判断的原因在于先按x进行排序的
            tmp[pos++]=i;
    }
    sort(tmp,tmp+pos,cmpy);
    for(int i=0;i<pos;i++)
        for(int j=i+1;j<pos&&poin[tmp[j]].y-poin[tmp[i]].y<=d;j++)             //将x之间的距离小于d的所有的点取出来以后开始再判断y的距离小于d的点，只有这样都满足，才进行计算。
    {
        double dd = dist(poin[tmp[i]],poin[tmp[j]]);
        d=min(d,dd);
    }
    return d;
}
int main()
{
    while(cin>>m&&m)
    {
        for(int i=0;i<m;i++)
            cin>>poin[i].x>>poin[i].y;
        sort(poin,poin+m,cmpxy);
        cout.setf(ios::fixed);
        double ans = min_dist(0,m-1)/2.0;
        cout<<setprecision(2)<<ans<<endl;
    }
}