HDOJ 1007 (分治法)

/** 用分治法求最近点对问题，时间复杂度为O(n*logn*logn) */
/** 最近点对问题，时间复杂度为O(n*logn*logn) */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const double INF = 1e20; //1e20 = 1^20 = 0.00000000000000000000000000000000000001 指最小的浮点数
const int N = 100005;
struct Point
{ double x; double y; }point[N];
int n;
int tmpt[N];
bool cmpxy(const Point& a, const Point& b)
{
if(a.x != b.x)
return a.x < b.x;
return a.y < b.y;
}

bool cmpy(const int& a, const int& b)
{
return point[a].y < point[b].y;
}

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

double dis(int i, int j) //求两点距离
{
return sqrt((point[i].x-point[j].x)*(point[i].x-point[j].x)
+ (point[i].y-point[j].y)*(point[i].y-point[j].y));
}

double Closest_Pair(int left, int right)
{
double d = INF;
if(left==right) //如果只有一个点，则最近点对的距离为0
return d;
if(left + 1 == right) //如果只有两个点，则返回两点之间的距离
return dis(left, right);
int mid = (left+right)>>1;
double d1 = Closest_Pair(left,mid);
double d2 = Closest_Pair(mid+1,right);
d = min(d1,d2);
int i,j,k=0;
//分离出宽度为d的区间
for(i = left; i <= right; i++)
{
if(fabs(point[mid].x-point[i].x) <= d) //fabs(float x) 求浮点数的绝对值
tmpt[k++] = i;
}
sort(tmpt,tmpt+k,cmpy);
//线性扫描
for(i = 0; i < k; i++)
{
for(j = i+1; j < k && point[tmpt[j]].y-point[tmpt[i]].y<d; j++)
{
double d3 = dis(tmpt[i],tmpt[j]);
if(d > d3)
d = d3;
}
}
return d;
}

int main()
{
while(true)
{
scanf("%d",&n);
if(n==0)
break;
for(int i = 0; i < n; i++)
scanf("%lf %lf",&point[i].x,&point[i].y);
sort(point,point+n,cmpxy);
printf("%.2lf\n",Closest_Pair(0,n-1)/2);
}
return 0;
}

//代码来源连接http://yzmduncan.iteye.com/blog/1432880

求坐标上最近两点的距离

HDOJ 1007 (分治法)相关推荐

详解分治法（divide-and-conquer）及其典型应用
什么是分治法在昨天的文章<漫谈数据库中的join>的最后,提到Grace hash join和Sort-merge join都是基于分治思想的.分治法(divide-and-conque ...
算法设计思想（4）— 分治法
1. 分治法概念分治,顾名思义,分而治之. 具体来说,它先将一个难以直接解决的大问题,分割成一些可以直接解决的小问题.如果分割后的问题仍然无法直接解决,那么就继续递归地分割,直到每个小问题都可解. ...
动态规划和分治法，贪心算法以及递归的再一次深刻理解和体会
每次体会算法都有新的感觉,刷题越多,对算法的理解感觉也就越深刻. 下面我们来重新体会下分治法,动态规划,贪心法,递归的理解. 1.分治法: 将问题分成单独的阶段,每个阶段互相不干扰很独立,如10米长的 ...
LeetCode 23. Merge k Sorted Lists--Python解法--优先队列，分治法
题目地址:Merge k Sorted Lists - LeetCode Merge k sorted linked lists and return it as one sorted list. A ...
Hanoi塔(分治法的应用)
1.分治法分治法的设计思想是将一个难以直接解决的大问题分解成一些规模较小的相同问题,以便各个击破,分而治之. 一般来说,分治算法在每一层递归上都有3个步骤: (1)分解:将问题分解成一系列子问题. ...
从合并排序算法看“分治法”
本文内容分治策略分治步骤从合并排序看"分治策略" 分治策略分治法(divide-and-conquer),"分治法策略"是一种很重要的算法.顾名思义,& ...
2_2 递归与分治策略（分治法的基本思想）
基本思想:将复杂问题简单化,大事化小,小事化了,很符合中国人的思想. 打官腔:分治算法的基本思想是将一个规模为N的问题分解为K个规模较小的子问题,这些子问题相互独立且与原问题性质相同.求出子问题的解, ...
python分治算法_python算法实现-分治法
分治法概念将一个复杂的问题分成两个或更多的相同或相似的子问题,再把子问题分成更小的子问题----"分" 将最后子问题可以简单的直接求解----"治" 将所有子问 ...
分治法解决最大子数组问题
分治法解决最大子数组问题参考文章: (1)分治法解决最大子数组问题 (2)https://www.cnblogs.com/Christal-R/p/Christal_R.html (3)https: ...

HDOJ 1007 (分治法)

HDOJ 1007 (分治法)相关推荐

最新文章

热门文章