最近点对

1. #include <iostream>  
2. #include <cstdio>  
3. #include <cstring>  
4. #include <cmath>  
5. #include <algorithm>  
6. using namespace std;  
7. const double INF = 1e20;  
8. const int N = 100005;  
9.   
10. struct Point  
11. {  
12.     double x;  
13.     double y;  
14. }point[N];  
15. int n;  
16. int tmpt[N];  
17.   
18. bool cmpxy(const Point& a, const Point& b)  
19. {  
20.     if(a.x != b.x)  
21.         return a.x < b.x;  
22.     return a.y < b.y;  
23. }  
24.   
25. bool cmpy(const int& a, const int& b)  
26. {  
27.     return point[a].y < point[b].y;  
28. }  
29.   
30. double min(double a, double b)  
31. {  
32.     return a < b ? a : b;  
33. }  
34.   
35. double dis(int i, int j)  
36. {  
37.     return sqrt((point[i].x-point[j].x)*(point[i].x-point[j].x)  
38.                 + (point[i].y-point[j].y)*(point[i].y-point[j].y));  
39. }  
40.   
41. double Closest_Pair(int left, int right)  
42. {  
43.     double d = INF;  
44.     if(left==right)  
45.         return d;  
46.     if(left + 1 == right)  
47.         return dis(left, right);  
48.     int mid = (left+right)>>1;  
49.     double d1 = Closest_Pair(left,mid);  
50.     double d2 = Closest_Pair(mid+1,right);  
51.     d = min(d1,d2);  
52.     int i,j,k=0;  
53.     //分离出宽度为d的区间  
54.     for(i = left; i <= right; i++)  
55.     {  
56.         if(fabs(point[mid].x-point[i].x) <= d)  
57.             tmpt[k++] = i;  
58.     }  
59.     sort(tmpt,tmpt+k,cmpy);  
60.     //线性扫描  
61.     for(i = 0; i < k; i++)  
62.     {  
63.         for(j = i+1; j < k && point[tmpt[j]].y-point[tmpt[i]].y<d; j++)  
64.         {  
65.             double d3 = dis(tmpt[i],tmpt[j]);  
66.             if(d > d3)  
67.                 d = d3;  
68.         }  
69.     }  
70.     return d;  
71. }  
72.   
73.   
74. int main()  
75. {  
76.     while(true)  
77.     {  
78.         scanf("%d",&n);  
79.         if(n==0)  
80.             break;  
81.         for(int i = 0; i < n; i++)  
82.             scanf("%lf %lf",&point[i].x,&point[i].y);  
83.         sort(point,point+n,cmpxy);  
84.         printf("%.2lf\n",Closest_Pair(0,n-1)/2);  
85.     }  
86.     return 0;  
87. }