HDU 1007 Quoit Design

644人阅读 评论(0) 收藏 举报
分类:

题目地址:点击打开链接

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


以下是代码实现,耗时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;
}


查看评论

杭电ACM第1007题——Quoit Design

WA #include #include #include using namespace std; const int SIZE = 100005; typedef struct { d...
  • miomiomiomio
  • miomiomiomio
  • 2017-03-23 14:30:27
  • 212

HDU-1007 Quoit Design

题目大意:给你n个点,代表toy,求一个使得可以套住一个toy且只能套住一个toy的环的最大半径。解法:只需要求出所有点的最近点对的长度就可以解决了。但是如何求最近点对? 最容易想到的就是暴力枚举所...
  • Snow_Me
  • Snow_Me
  • 2016-05-20 10:14:30
  • 459

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

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) ...
  • sunrainchy
  • sunrainchy
  • 2013-07-19 11:12:34
  • 786

HDU 1007

  • sky_miange
  • sky_miange
  • 2014-12-31 13:25:11
  • 703

杭电 ACM 1007 Quoit Design

Quoit DesignTime Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Tota...
  • hsm159
  • hsm159
  • 2017-04-26 12:36:00
  • 147

hdu1007最近点对-分治策略-具体优化

1)AC总是让人感到兴奋。 上一篇求最大连续子数组之和了解到分治策略之后,本题先是找到中间界线,依然分成三部分,左边集合、右边集合、两个点分别在左边和右边的集合的情况。前两种情况不断递归,直到该集合内...
  • a272846945
  • a272846945
  • 2016-03-08 21:23:52
  • 436

HDU 1007 Quoit Design (分治)

HDU 1007 Quoit Design (分治) 题目大意: 给你很多点,问你最近点的距离一半是多少。 解题思路: 最近点对,参照:http://blog.csdn.net/hellobaby...
  • a1061747415
  • a1061747415
  • 2014-06-04 19:14:08
  • 705

hdu 1007 && zoj 2107 Quoit Design [最近点对]

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1007                    http://acm.zju.edu.cn/on...
  • u011394362
  • u011394362
  • 2014-08-11 17:48:07
  • 640

hdu1007 Quoit Design 分治

Quoit DesignTime Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Tota...
  • liu_tian_wei
  • liu_tian_wei
  • 2017-07-10 17:26:22
  • 418

hdu 1007 zoj 2107 Quoit Design 求平面最近点对 分治法

本题题目很裸,只要求出平面中最近点对的距离,然后除以2就是答案,方法就是分治,具体的可以参见算法导论 简单分析:对于该问题的一个子问题,只需把当前点的集合再分成两部分,分的方法是:按照x(或者y)排序...
  • cwb_where
  • cwb_where
  • 2012-08-13 23:50:28
  • 1715
    个人资料
    专栏达人 持之以恒
    等级:
    访问量: 6万+
    积分: 2252
    排名: 2万+
    联系方式
    欢迎谈论交流:1245985209
    博客专栏