HDU1007Quoit Design（分治）

真令人头大，按照这个博客的代码写的https://blog.csdn.net/wr_technology/article/details/51057937

一直超时。

最后发现不能用cin,原因如下：

scanf是格式化输入，printf是格式化输出。
cin是输入流，cout是输出流。效率稍低，但书写简便。

效率低的原因可以自行百度。

数据量大的题目最好一开始就用scanf，printf写，不要等超时了再回来改。

 

 

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <string>
#include <algorithm>
#include <stdlib.h>
#include <cmath>


    
using namespace std;

const int MAXN = 100005;
const double INF=1e20;
struct Point{
    double x;
    double y;
};
Point point[MAXN];
int index[MAXN];
double dis(Point a,Point b){
    
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmpx(Point x,Point y){
    return x.x<y.x;
}
bool cmpy(int a,int b){
    return point[a].y<point[b].y;
}

double closest_point(int s,int e){
    double d=INF;    if(s==e) return d; 
else if(s+1==e) return dis(point[s],point[e]);

    int mid = (s+e)>>1;//相当于(s+e)/2
    d = min(closest_point(s,mid),closest_point(mid+1,e));
    //递归分解问题，找到子区域中的最近点距离
     /***********找到后，以他们的最近距离来分割区间***********/
    int i,j,cut = 0;
    for(i=s;i<=e;i++){
        if(fabs(point[i].x-point[mid].x)<d)//abs（）函数主要是对int求绝对值
            index[cut++]=i;                  //fabs（）对float，double求绝对值
    } 
    //在区间里找x宽度小于最近距离的把他的下标存在index数组里面
    sort(index,index+cut,cmpy);
    for(i=0; i<cut; i++){
        for(j=i+1; j<cut; j++){
            if(fabs(point[index[i]].y-point[index[j]].y)>=d)
                break;
            //sort之后，只要当前超过了了，后面的数字一点不可能比这个小，所以不找，省时间
            d=min(dis(point[index[i]],point[index[j]]),d);
        }
    }
    return d;
}

int main(){
    int n;
    while(scanf("%d",&n)!=EOF&&n){
        int i;
        for(i=0;i<n;i++)
        //    cin>>point[i].x>>point[i].y;
          scanf("%lf %lf",&point[i].x,&point[i].y);
        sort(point,point+n,cmpx);
        printf("%.2lf\n",closest_point(0,n-1)/2);
        //cout<<setiosflags(ios::fixed)<<setprecision(2)<<getMin(0,n-1)/2<<endl;
    }
}