HDU 1007 Quoit Design

原创 2017年05月11日 18:43:37

题目地址:点击打开链接

一道典型的分治法求最小点对的问题,很难搞懂,请参考大神的博客:点击打开链接


以下是代码实现,耗时889ms,已经相当快了,注意如果用cin输入的话会超时。

/*分治法求最小点对*/ 
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iomanip>

using namespace std;
struct point {
	double x,y;
}p[100005],*px[100005],*py[100005];

double cmpx(point *a,point *b)
{
	if(a->x==b->x)
	return a->y<b->y;
	
	return a->x<b->x;
}

double cmpy(point *a,point *b)
{
	return a->y<b->y;
}

double dis(point *a,point *b)
{
	return sqrt((a->x-b->x)*(a->x-b->x)+(a->y-b->y)*(a->y-b->y));
}

double min(double a,double b)
{
	return a<b?a:b;
} 

double find(int s,int e)
{
	if(s+1==e)
	return dis(px[s],px[e]);
	if(s+2==e)  
    return min(dis(px[s],px[s+1]),min(dis(px[s+1],px[e]),dis(px[s],px[e])));
    
    int mid=(s+e)>>1;
    double ans=min(find(s,mid),find(mid+1,e));
    int i,j,cnt=0;  
    for(i=s;i<=e;i++)
    {  
        if(px[i]->x>=px[mid]->x-ans&&px[i]->x<=px[mid]->x+ans)  
        py[cnt++]=px[i];  
    }  
    sort(py,py+cnt,cmpy);//按y坐标排序  
    for(i=0;i<cnt;i++)  
    {  
        for(j=i+1;j<cnt;j++) 
        {  
            if(py[j]->y-py[i]->y>=ans)  
            break;  
            ans=min(ans,dis(py[i],py[j]));  
        }  
    }  
    return ans;  
	
}

int main()
{
	long long n,i;
	while(scanf("%d",&n)!=EOF&&n)
	{
		for(i=0;i<n;i++)
		{
			scanf("%lf%lf",&p[i].x,&p[i].y);
			px[i]=&p[i];
		}
		
		sort(px,px+n,cmpx);
		double dis=find(0,n-1)/2;
		
		cout<<fixed<<setprecision(2)<<dis<<endl;
	}
	return 0;
}


版权声明:欢迎转载,请注明此博客地址。

相关文章推荐

HDU 1007 Quoit Design

给一系列点的坐标,求距离最小的两个点。基本是照着别人的代码写的……我以为我看懂了的……怎么还是错了辣么多次……#include #include #include #include using ...
  • wen_ruo
  • wen_ruo
  • 2015年03月12日 22:07
  • 261

hdu 1007 Quoit Design 最近点对

一道题目,最近点对的模版题。

hdu1007 Quoit Design(最近点对模板题)

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) ...

HDU 1007 Quoit Design

题意: 给你n个全在一个平面上点的坐标,
  • wwx0848
  • wwx0848
  • 2014年09月27日 01:50
  • 277

hdu 1007 Quoit Design(分治法求最近点对)

大致题意:给N个点,求最近点对的距离 d ;输出:r = d/2。 // Time 2093 ms; Memory 1812 K #include #include #include #includ...

HDU 1007 Quoit Design 分治法

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) ...

HDU 1007 Quoit Design(分治)

题目链接:点击打开链接 思路: 经典的分治法, 网上讲解很多我就不多说了, 这是nlognlogn复杂度, 大多数情况是够用的。。优化了一下排序函数, 跑了780ms 细节参见代码: #inc...

HDU 1007 Quoit Design(计算几何 平面最近点对)

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) ...

hdu 1007 Quoit Design 最近点对(分治)

题意:给定一个整数n,以及n个x

HDU1007 Quoit Design 【分治】

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) ...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 1007 Quoit Design
举报原因:
原因补充:

(最多只允许输入30个字)