HDU 1875 prim算法复习

思想一开始就对了，不过在置无穷大的时候犯了个错误，这个题目的路程距离dis是double型的，以前写的都是int 类型的，置无穷大都是用memset(map,0x3f,sizeof(map));怪我乱来，double也这么置无穷大，结果就出错了。当然题目还是典型的最小生成树问题，练练手，不然忘记了；。。。

畅通工程再续

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

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

  
  

   
   2
2
10 10
20 20
3
1 1
2 2
1000 1000
  
  

 

Sample Output

  
  

   
   1414.2
oh!
  
  

 

代码：

#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstdio>
#define maxn 102
#define INF 999999999.0
using namespace std;

double map[maxn][maxn],dis[maxn];
int visit[maxn],n;
struct pointer
{
	double x;
	double y;
}point[maxn];
double prim()
{
	for(int i = 2; i <= n; i++)
	{
		dis[i] = map[1][i];
	}
	dis[1] = 0;
	visit[1] = 1;
	double sum = 0.0;
	for(int i = 1; i <= n-1; i++)
	{
		double temp = INF;
		int pos = -1;
		for(int j = 1; j <= n; j++)
		{
			if(!visit[j] && temp > dis[j])
			{
				temp = dis[j];
				pos = j;
			}
		}
		if(pos == -1)return 0.0;
		visit[pos] = 1;
		sum += dis[pos];
		for(int j = 1; j <= n; j++)
		{
			if(!visit[j] && map[pos][j] < dis[j])
				dis[j] = map[pos][j];
		}
	}
	return sum;
}
int main()
{
	int T;
	cin>>T;
	while(T--)
	{
		cin>>n;
		double x,y;
		for(int i = 1; i <= n; i++)
		{
			scanf("%lf%lf",&point[i].x,&point[i].y);
		}
		memset(visit,0,sizeof(visit));
		for(int i = 1; i <= n; i++)
		{
			double xx,yy,distan;
			map[i][i] = 0.0;
			for(int j = i+1; j <= n; j++)
			{
				map[i][j] = map[j][i] = INF;
				xx = (point[i].x-point[j].x)*(point[i].x-point[j].x);
				yy = (point[i].y-point[j].y)*(point[i].y-point[j].y);
			 	distan  = sqrt(xx+yy);
				if((distan>=10.0) && (distan<=1000.0))
					map[i][j] = map[j][i] = distan;
			}
		}
	
		double ans = prim()*100;
		if(ans == 0.0)puts("oh!");
		else
			printf("%.1lf\n",ans);
	}
}