hdu 1875 畅通工程再续 【最小生成树-prim】

畅通工程再续

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 26547    Accepted Submission(s): 8613

Problem Description

相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其他的小岛时都要通过划小船来实现。现在政府决定大力发展百岛湖，发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！经过考察小组RPRush对百岛湖的情况充分了解后，决定在符合条件的小岛间建上桥，所谓符合条件，就是2个小岛之间的距离不能小于10米，也不能大于1000米。当然，为了节省资金，只要求实现任意2个小岛之间有路通即可。其中桥的价格为 100元/米。

 

Input

输入包括多组数据。输入首先包括一个整数T(T <= 200)，代表有T组数据。
每组数据首先是一个整数C(C <= 100),代表小岛的个数，接下来是C组坐标，代表每个小岛的坐标，这些坐标都是 0 <= x, y <= 1000的整数。

 

Output

每组输入数据输出一行，代表建桥的最小花费，结果保留一位小数。如果无法实现工程以达到全部畅通，输出”oh!”.

 

Sample Input

2210 1020 2031 12 21000 1000

 

Sample Output

1414.2oh!

小总结：最小生成树的典型应用。

AC代码：

# include <stdio.h># include <string.h># include <math.h># define INF 1000000# define M 205struct Dao  // 用来存放小岛的坐标{int x;int y;};Dao d[M];double map[M][M];double lowcost[M];double sum;int closevertex[M];int n;double Distance(int x1, int y1, int x2, int y2)  //计算两点之间的距离{return sqrt((1.0 * x1-x2) * (1.0 * x1-x2) + (1.0 * (y1 - y2)) * (1.0 * (y1 - y2)));}void Change()  //把坐标数据转换成带权值的有向图{int i, j;for (i = 0; i < n; i++){for (j = 0; j < n; j++){double t = Distance(d[i].x, d[i].y, d[j].x, d[j].y);if (t <= 1000.0 && t >= 10.0){map[i][j] = map[j][i] = t;}}}}void Prim()  //最小生成树的prim算法{int i, j, k;double mincost;int count = 1;for (i = 0; i < n; i++){lowcost[i] = map[0][i];closevertex[i] = 0;}lowcost[0] = 0;for (i = 1; i < n; i++){mincost = INF;j = 1;k = -1;while (j < n){if (lowcost[j] < mincost && lowcost[j] != 0){mincost = lowcost[j];k = j;}j++;}if (k != -1){sum += mincost;lowcost[k] = 0;count ++;}for (j = 1; j < n; j++){if (map[k][j] < lowcost[j]){lowcost[j] = map[k][j];closevertex[j] = k;}}}if (count == n) //判断是否能把所有小岛都相连{printf("%.1lf\n", sum * 100);}else{printf("oh!\n");}}int main(void){int t;scanf("%d", &t);while (t--){int i, j;sum = 0.0;scanf("%d", &n);for (i = 0; i < n; i++) //初始化地图{for (j = 0; j < n; j++){if (i == j){map[i][j] = 0.0;}else{map[i][j] = map[j][i] = INF;}}}for (i = 0; i < n; i++){scanf("%d %d", &d[i].x, &d[i].y);}Change();Prim();}return 0;}