最近点对（蛮力法）

 

1. #include <stdio.h>
2. //rand() and srand()
3. #include <stdlib.h>
4. #include <time.h>
5. //sqrt()
6. #include <math.h>
7. //INT_MAX
8. #include <limits.h>
10. #define distance(x1,y1,x2,y2)   /
11.   sqrt((x1-x2)*(x1-x2) + (y1-y2)*(y1-y2))
12. #define MAX 10          //生成点的个数
14. struct point
15. {
16.   int x;
17.   int y;
18. }point[MAX];
20. struct pair
21. {
22.   int point1;
23.   int point2;
24.   int min_distance;
25. };
27. void init_point(struct point *p);
28. void show_point(struct point *p);
29. struct pair closest_pair(struct point *p);
31. int main(void)
32. {
33.   struct pair point_pair;
34.   init_point(point);
35.   show_point(point);
37.   point_pair = closest_pair(point);
38.   
39.   printf("/nThe closest points are point %d and point %d./n", point_pair.point1+1, point_pair.point2+1 );
40.   printf("The distance is %d/n", point_pair.min_distance);
41.   return 0;
42. }//main()
44. void init_point(struct point *p)
45. {
46.   int i;
47.   srand(time(NULL));
49.   for (i=0; i<MAX; i++)
50.   {
51.     p[i].x = rand()%100;
52.     p[i].y = rand()%100;
53.   }//for
55. }//init_point()
57. void show_point(struct point *p)
58. {
59.   int i;
61.   for (i=0; i<MAX; i++)
62.   {
63.     printf("point[%d]/t(%d,%d)/n", i+1, p[i].x,p[i].y);
64.   }//for
65. }//show_point()
67. struct pair closest_pair(struct point *p)
68. {
69.   int i,j;
70.   int temp;
71.   struct pair point_pair;
73.   point_pair.min_distance = INT_MAX;
75.   for (i=0; i<MAX-1; i++)
76.   {
77.     for (j=i+1; j<MAX; j++)
78.     {
79.       temp = distance(p[i].x,p[i].y,p[j].x,p[j].y);
80.       if (temp < point_pair.min_distance)
81.       {
82.         point_pair.min_distance = temp;
83.         point_pair.point1=i;
84.         point_pair.point2=j;
85.       }//if
86.     }//inner_for
87.   }//for()
89.   return point_pair;
90. }//closest_pair()