畅通工程再续
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的整数。
每组数据首先是一个整数C(C <= 100),代表小岛的个数,接下来是C组坐标,代表每个小岛的坐标,这些坐标都是 0 <= x, y <= 1000的整数。
Output
每组输入数据输出一行,代表建桥的最小花费,结果保留一位小数。如果无法实现工程以达到全部畅通,输出”oh!”.
Sample Input
2 2 10 10 20 20 3 1 1 2 2 1000 1000
Sample Output
1414.2oh!
小总结:最小生成树的典型应用。
AC代码:
# include <stdio.h>
# include <string.h>
# include <math.h>
# define INF 1000000
# define M 205
struct 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;
}