最接近点问题的分治法

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#include <float.h>



// 定义一个二维点的结构体

struct Point {

    int x, y;

};



// 比较函数，用于根据x坐标对点进行排序

int compareX(const void* a, const void* b) {

    return ((struct Point*)a)->x - ((struct Point*)b)->x;

}



// 比较函数，用于根据y坐标对点进行排序

int compareY(const void* a, const void* b) {

    return ((struct Point*)a)->y - ((struct Point*)b)->y;

}



// 计算两点之间的距离

float dist(struct Point p1, struct Point p2) {

    return sqrt((p1.x - p2.x)*(p1.x - p2.x) + (p1.y - p2.y)*(p1.y - p2.y));

}



// 计算最小距离的函数

float min(float x, float y) {

    return (x < y) ? x : y;

}



// 在给定点集中找到最小距离的函数

float bruteForce(struct Point P[], int n) {

    float min_distance = FLT_MAX;

    for (int i = 0; i < n; ++i) {

        for (int j = i+1; j < n; ++j) {

            if (dist(P[i], P[j]) < min_distance) {

                min_distance = dist(P[i], P[j]);

            }

        }

    }

    return min_distance;

}



// 使用分治法找到最小距离的函数

float closestUtil(struct Point Px[], struct Point Py[], int n) {

    if (n <= 3) {

        return bruteForce(Px, n);

    }



    // 找到中间点

    int mid = n / 2;

    struct Point midPoint = Px[mid];



    // 分别在左右子集中递归查找最小距离

    float dl = closestUtil(Px, Py, mid);

    float dr = closestUtil(Px + mid, Py + mid, n - mid);



    // 取两个子集中的最小距离

    float d = min(dl, dr);



    // 构建y坐标在(d - 2d, d + 2d)范围内的点集

    struct Point strip[n];

    int j = 0;

    for (int i = 0; i < n; i++) {

        if (abs(Py[i].x - midPoint.x) < d) {

            strip[j] = Py[i];

            j++;

        }

    }



    // 找到strip中的最小距离

    float minStrip = bruteForce(strip, j);



    // 返回最小距离

    return min(d, minStrip);

}



// 主函数，计算最小距离

float closest(struct Point P[], int n) {

    struct Point Px[n];

    struct Point Py[n];

    for (int i = 0; i < n; i++) {

        Px[i] = P[i];

        Py[i] = P[i];

    }



    // 按x坐标排序

    qsort(Px, n, sizeof(struct Point), compareX);

    // 按y坐标排序

    qsort(Py, n, sizeof(struct Point), compareY);



    // 使用分治法计算最小距离

    return closestUtil(Px, Py, n);

}



// 示例用法

int main() {

    struct Point P[] = {{2, 3}, {12, 30}, {40, 50}, {5, 1}, {12, 10}, {3, 4}};

    int n = sizeof(P) / sizeof(P[0]);

    printf("最接近点的距离是 %f\n", closest(P, n));

    return 0;

}