寻找平面上最近的点对

最新推荐文章于 2022-05-07 16:47:44 发布

mingxiaoshan

最新推荐文章于 2022-05-07 16:47:44 发布

阅读量295

点赞数

分类专栏： c、c++ 数据结构

本文链接：https://blog.csdn.net/u012979457/article/details/102567599

版权

这篇博客介绍了如何使用预排序和分治策略解决寻找平面上最近点对的问题，时间复杂度达到O(nlogn)。作者在实现过程中针对横坐标重复的情况调整了排序规则，修复了原有BUG，并欢迎读者指正。

摘要由CSDN通过智能技术生成

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define MAX_DIS 65536

struct point{
    int x;
    int y;
};

char result[64]="\0";

int cmp_x(const void * p1, const void * p2)
{
    if(((struct point *)p1)->x==((struct point *)p2)->x)
        return ((struct point *)p1)->y-((struct point *)p2)->y;
    else
        return ((struct point *)p1)->x-((struct point *)p2)->x;
}

int cmp_y(const void * p1, const void * p2)
{
    if(((struct point *)p1)->y==((struct point *)p2)->y)
        return ((struct point *)p1)->x-((struct point *)p2)->x;
    else
        return ((struct point *)p1)->y-((struct point *)p2)->y;
}

double min(double x, double y)
{
    if(x<y)
        return x;
    else
        return y;
}

double cal_dis(struct point *p1, struct point *p2)
{
    return sqrt((p1->x-p2->x)*(p1->x-p2->x)+(p1->y-p2->y)*(p1->y-p2->y));
}


void generate_points(struct point * P, int n)
{
    int i;
    for(i=0;i<n;i++)
    {
        P[i].x=rand()%100-50;
        P[i].y=rand()%100-50;
    }
}

void print_point(struct point * P, int n)
{
    int i;
    for(i=0;i<n;i++)
        printf("(%d,%d)",P[i].x,P[i].y);
    printf("\n\n");
}

void copy_points(struct point * P1, struct point * P2, int n)
{
    int i;
    for(i=0;i<n;i++