畅通project再续


Problem Description
相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其它的小岛时都要通过划小船来实现。

如今政府决定大力发展百岛湖。发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！

经过考察小组RPRush对百岛湖的情况充分了解后。决定在符合条件的小岛间建上桥，所谓符合条件，就是2个小岛之间的距离不能小于10米，也不能大于1000米。

当然，为了节省资金。仅仅要求实现随意2个小岛之间有路通就可以。当中桥的价格为 100元/米。

 

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

Output
每组输入数据输出一行，代表建桥的最小花费。结果保留一位小数。假设无法实现project以达到所有畅通，输出”oh!”.
 

Sample Input
2
2
10 10
20 20
3
1 1
2 2
1000 1000
 

Sample Output
1414.2
oh!
 

Author
8600
 

Source
2008浙大研究生复试热身赛（2）——全真模拟
 



AC代码：

kruskal：

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<algorithm>
using namespace std;
int per[111];
double x[111],y[111];
struct node{
	int start,end;
	double distance;
};
node p[10000];
int cmp(node a,node b)	//距离从小到大 
{
	return a.distance < b.distance;
}
void init()
{
	for(int i = 1;i < 111;i++ )
		per[i] = i;
}
int find(int x)
{
	if(x == per[x])
		return x;
	return per[x] = find(per[x]);	
}
bool join(int x,int y)
{
	int fx = find(x);
	int fy = find(y);
	if(fx != fy)  //推断是否成环 
	{
		per[fx] = fy;
		return true;  //没成环 
	}
	return false;	
}
int main()
{
	int i,j,k;
	int t,n,c;
	double d,cost;  //一定要注意题目所要求的数据类型 
	scanf("%d",&t);
	while(t--)
	{
		init();
		scanf("%d",&c);
		for(i = 1;i <= c; i++)
		{
			scanf("%lf%lf",&x[i], &y[i]);
		}
		k = 0;
		for(i = 1;i <= c;i++)
		{
			for(j = i+1;j <= c;j++)
			{
				d = sqrt( (x[j]-x[i])*(x[j]-x[i])+(y[j]-y[i])*(y[j]-y[i]) );//距离公式 
				if(d >=10.0 && d <= 1000.0)	//一定要筛选完再存进结构体 要不然会WA.
				{							
					p[k].start = i;   //起点 
					p[k].end = j;	  //终点 
					p[k].distance = d;//起点到到终点的距离 
					k++;
				}
			}
		}
		sort(p,p+k,cmp);   //按距离从小到大排序 
		int num = 0;
		cost = 0.0;
		for( i = 0;i < k; i++)
		{
			if(join(p[i].start, p[i].end))
			cost += p[i].distance;
		}
		for( i = 1;i <= c;i++ )
		{
			if(per[i]==i)
			num++;
			if(num > 1)		// 及时跳出循环节省时间 
			break;
		}
		if(num>1)  //无法连接全部小岛 
			printf("oh!\n");
		else
			printf("%.1lf\n",100*cost);
	}
	return 0;
}



===============================切割线======================================

prim：

#include<cstdio>
#include<cstring>
#include<cmath>
#define mem(a, b) memset(a, (b), sizeof(a))
#define Wi(a) while(a--)
#define Si(a) scanf("%d", &a)
#define Pi(a) printf("%d\n", (a))
#define Pf(a) printf("%.1lf\n", (a))
#define INF 0x3f3f3f
double map[150][150];
int x[150],y[150];
int n;
double d[150], vis[150];
void prim()
{
	mem(vis, 0);
	int i, j, k;
	double ans = 0, minn;
	for(i = 1; i <= n; i++)
		d[i] = map[1][i];
	vis[1] = 1;
	for(i = 1; i < n; i++)
	{
		k = 1;
		minn = INF;
		for(j = 1; j <= n; j++)
		{
			if(!vis[j] && d[j] < minn)
			{
				minn = d[j];
				k = j;
			}
		}
		if(minn == INF){
			puts("oh!");return;
		}
		vis[k] = 1;
		ans += minn;
		for(j = 1; j <= n; j++)
		{
			if(!vis[j] && d[j] > map[j][k])
				d[j] = map[j][k];
		}
	}
	Pf(ans*100);
} 
int main(){
	int t; Si(t);
	Wi(t){
		Si(n);
		int i , j, k;
		for(i = 1; i <= n; i++)
		{
			scanf("%d%d", &x[i], &y[i]);
		}
		for(i = 1; i <= n; i++)
		{
			for(j = 1; j <= n; j++)
			{
				double d = sqrt( (x[i]-x[j])*(x[i]-x[j]) + (y[i]-y[j])*(y[i]-y[j]));
				if(d >= 10.0 && d <= 1000.0)
					map[i][j] = map[j][i] =  d;
				else
					map[i][j] = map[j][i] = INF;
			}
			map[i][i] = 0;
		}
		prim();
	}
	return 0;
}