Quoit Design---hdu1007(最近点对问题 分治法)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1007

题意:给你n(2<=n<=10^6)个点的坐标,然后找到两个点使得他们之间的距离最小,然后输出最小距离的一半;

先把n个点按x坐标排序,然后求左边n/2个和右边n/2个的最近距离,最后合并。

首先,假设点是n个,编号为1到n。我们要分治求,则找一个中间的编号mid,先求出1到mid点的最近距离设为d1,还有mid+1到n的最近距离设为d2。这里的点需要按x坐标的顺序排好,并且假设这些点中,没有2点在同一个位置。(若有,则直接最小距离为0了)。

然后,令d为d1, d2中较小的那个点。如果说最近点对中的两点都在1-mid集合中,或者mid+1到n集合中,则d就是最小距离了。但是还有可能的是最近点对中的两点分属这两个集合,所以我们必须先检测一下这种情况是否会存在,若存在,则把这个最近点对的距离记录下来,去更新d。这样我们就可以得道最小的距离d了。

关键是要去检测最近点对,理论上每个点都要和对面集合的点匹配一次,那效率还是不能满足我们的要求。所以这里要优化。怎么优化呢?考虑一下,假如以我们所选的分割点mid为界,如果某一点的横坐标到点mid的横坐标的绝对值超过d1并且超过d2,那么这个点到mid点的距离必然超过d1和d2中的小者,所以这个点到对方集合的任意点的距离必然不是所有点中最小的。

    所以我们先把在mid为界左右一个范围内的点全部筛选出来,放到一个集合里。筛选好以后,当然可以把这些点两两求距离去更新d了,不过这样还是很慢,万一满足条件的点很多呢。这里还得继续优化。首先把这些点按y坐标排序。假设排序好以后有cnt个点,编号为0到cnt-1。那么我们用0号去和1到cnt-1号的点求一下距离,然后1号和2到cnt-1号的点求一下距离。。。如果某两个点y轴距离已经超过了d,这次循环就可以直接break了,开始从下一个点查找了.

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <bitset>
#include <iostream>
#include <time.h>
#include <vector>
#include <queue>

typedef long long LL;

using namespace std;

const int N = 1e6+1;
const double eps = 1e-10;
const int INF = 0x3f3f3f3f;
const int mod = 1000000007;
const double PI = 4*atan(1.0);

struct point
{
    double x, y;
}p[N], t[N];

double dist(point p1, point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y));
}

bool cmpX(point p1, point p2)
{
    return p1.x < p2.x;
}
bool cmpY(point p1, point p2)
{
    return p1.y < p2.y;
}

double Find(int L, int R)
{
    if(L+1 == R)///当只有两个点的时候;
        return dist(p[L], p[R]);
    if(L+2 == R)///三个点的时候;
        return min(min(dist(p[L], p[L+1]), dist(p[L], p[R])), dist(p[L+1], p[R]));

    int Mid = (L+R)/2;

    double Min_d  = min(Find(L, Mid), Find(Mid+1, R));///找到两边的最小值, 下面更新中间部分的;

    int cnt = 0;

    for(int i=L; i<=R; i++)
    {
        if(fabs(p[i].x-p[Mid].x) <= Min_d)
            t[cnt++] = p[i];///把可能是最近点对中的点加入t集合中去;
    }

    sort(t, t+cnt, cmpY);///排序,下面找到,t集合中最近的两点间的距离,跟新Min_d;

    for(int i=0; i<cnt; i++)
    {
        for(int j=i+1; j<cnt; j++)
        {
            if( t[j].y - t[i].y > Min_d )
                break;
            Min_d = min(Min_d, dist(t[i], t[j]));
        }
    }
    return Min_d;
}

int main()
{
    int n;

    while(scanf("%d", &n), n)
    {
        for(int i=0; i<n; i++)
            scanf("%lf %lf", &p[i].x, &p[i].y);

        sort(p, p+n, cmpX);///按横坐标x排序;

        double ans = Find(0, n-1);///递归 求解;

        printf("%.2f\n", ans/2);
    }
    return 0;
}

时间: 2024-10-05 04:45:05

Quoit Design---hdu1007(最近点对问题 分治法)的相关文章

Quoit Design(hdu1007)最近点对问题。模版哦!

Quoit Design Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 30919 Accepted Submission(s): 8120 Problem Description Have you ever played quoit in a playground? Quoit is a game in which flat rings

Quoit Design(最近点对+分治)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1007 Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 42865    Accepted Submission(s): 11128 Problem Description Have you ever played

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

Have you ever played quoit in a playground? Quoit is a game in which flat rings are pitched at some toys, with all the toys encircled awarded. In the field of Cyberground, the position of each toy is fixed, and the ring is carefully designed so it ca

【HDOJ】P1007 Quoit Design (最近点对)

题目意思很简单,意思就是求一个图上最近点对. 具体思想就是二分法,这里就不做介绍,相信大家都会明白的,在这里我说明一下如何进行拼合. 具体证明一下为什么只需要检查6个点 首先,假设当前左侧和右侧的最小值为d,那么对于一个点,如果有个最小值小于d,那么一定存在于上d下d左d右d的一块区域内,又因为是从左到右,从上到下,所以左侧的那部分匹配的区域会重叠,也就是对于左侧的区域,完全没有必要去进行匹配.所以只需要对右侧d,上下d的区域进行匹配,而假设这个区域内的所以点的距离为d那么最多为6个,那么如果长

HDOJ-1007 Quoit Design(最近点对问题)

http://acm.hdu.edu.cn/showproblem.php?pid=1007 给出n个玩具(抽象为点)的坐标 求套圈的半径 要求最多只能套到一个玩具 实际就是要求最近的两个坐标的距离 典型的最近点对问题 最近点对详解 http://blog.csdn.net/lonelycatcher/article/details/7973046 //最近点对 # include <stdio.h> # include <algorithm> # include <math

hdu 1007 Quoit Design, 平面最近点对

平面最近点对的经典做法就是分治.所有点先按x再按y排序,然后取中间位置的点,利用其x坐标值将点划分成左右两部分,分别求出两侧的最近点对(设其距离为δ),然后合并.参看下面的图片(来自guyulongcs的专栏,http://blog.csdn.net/guyulongcs/article/details/6841550)合并的时候需要注意,大致有两种做法,一种是把分界线两侧δ范围内的点分别加进一个数组,然后对右侧数组按y坐标值排序,之后依次为每个左侧数组中的点,在右侧数组中寻找与其y坐标值最近的

HDU ACM 1007 Quoit Design 分治法求最近点对

题意:给n个点的坐标,求距离最近的一对点之间距离的一半. 分析:分治法求最近点对. #include<iostream> #include<algorithm> #include<cmath> using namespace std; #define N 100005 double min(double a,double b) { return a<b?a:b; } struct POINT { double x,y; }; POINT point[N],*px[

hdu 1007 Quoit Design 分治求最近点对

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 29344    Accepted Submission(s): 7688 Problem Description Have you ever played quoit in a playground? Quoit is a game in which flat

HDU1007 Quoit Design 【分治】

Quoit Design Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 30505    Accepted Submission(s): 8017 Problem Description Have you ever played quoit in a playground? Quoit is a game in which flat