poj2349——Arctic Network

题目大意：给出P个哨所的坐标，他们相互可以在D半径的范围内通信，其中有S个哨所还配备了卫星，有卫星的哨所之间通信不受范围限制，问D的最小值

输入：（可以有很多case）

            case个数N

            S  P

            第i个哨所的坐标xi  yi（共P行）

输出：D的最小值

分析：最小生成树的应用，最小生成树中最大的S-1条边可以分配给卫星哨所，剩下的边中最大的就是D的长度。

代码：转载自http://blog.youkuaiyun.com/lyhvoyage/article/details/19930551

1. #include<stdio.h>  
2. #include<math.h>  
3. #include<algorithm>  
4. using namespace std;  
5. int father[550], m, k;  
6. double d[550];  
7. struct post  
8. {  
9.     double x, y;  
10. }p[550];  
11. struct edge  
12. {  
13.     int u, v;  
14.     double w;  
15. }e[500005];  
16. bool comp(edge e1, edge e2)  
17. {  
18.     return  e1.w < e2.w;  
19. }  
20. double get_dis(double x1, double y1, double x2, double y2)  
21. {  
22.     return sqrt((x1-x2)*(x1-x2) + (y1-y2)*(y1-y2));  
23. }  
24. void Init(int n)  
25. {  
26.     for(int i = 1; i <= n; i++)  
27.         father[i] = i;  
28. }  
29. int Find(int x)  
30. {  
31.     if(x != father[x])  
32.         father[x] = Find(father[x]);  
33.     return father[x];  
34. }  
35. void Merge(int a, int b)  
36. {  
37.     int p = Find(a);  
38.     int q = Find(b);  
39.     if(p > q)  
40.         father[p] = q;  
41.     else  
42.         father[q] = p;  
43. }  
44. void Kruskal(int n)  
45. {  
46.     k = 0;  
47.     double Max = 0;  
48.     for(int i = 0; i < m; i++)  
49.         if(Find(e[i].u) != Find(e[i].v))  
50.         {  
51.             Merge(e[i].u, e[i].v);  
52.             d[k++] = e[i].w;  
53.             n--;  
54.             if(n == 1)  
55.                 return;  
56.         }  
57. }  
58. int main()  
59. {  
60.     int t, S, P, i, j;  
61.     double x, y;  
62.     scanf("%d",&t);  
63.     while(t--)  
64.     {  
65.         m = 0;  
66.         scanf("%d%d",&S,&P);  
67.         Init(P);  
68.         for(i = 1; i <= P; i++)  
69.             scanf("%lf%lf",&p[i].x, &p[i].y);  
70.         for(i = 1; i <= P; i++)  
71.             for(j = i + 1; j <= P; j++)  
72.             {  
73.                 e[m].u = i;  //感觉边存一次就够了，可以试一试
74.                 e[m].v = j;  
75.                 e[m++].w = get_dis(p[i].x, p[i].y, p[j].x, p[j].y);  
76.                 e[m].u = j;  
77.                 e[m].v = i;  
78.                 e[m++].w = get_dis(p[i].x, p[i].y, p[j].x, p[j].y);  
79.             }  
80.         sort(e, e+m, comp);  
81.         Kruskal(P);  
82.         printf("%.2lf\n",d[P-S-1]);  
83.     }  
84.     return 0;  
85. }