HDOJ---1875 畅通工程再续[Prim()算法]

畅通工程再续

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


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

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

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

[align=left]Sample Input[/align]

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

[align=left]Sample Output[/align]

1414.2 oh!

[align=left]Author[/align]
8600

[align=left]Source[/align]
2008浙大研究生复试热身赛（2）——全真模拟

[align=left]Recommend[/align]
lcy

code:

#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>
#include <algorithm>
#include <string>
#include <set>
#include <utility>
#include <queue>
#include <stack>
#include <list>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <ctype.h>
using namespace std;

#define MAXN 1005

int t,c;
double map[MAXN][MAXN];
double dis[MAXN];
int vst[MAXN];

typedef struct island
{
double x,y;
}Island;
Island island[MAXN];

double fun_dis(double x1,double y1,double x2,double y2)
{
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}

double prim()
{
int i,j,k;
double mindis,sum=0;
for(i=2;i<=c;i++)
{
vst[i]=0;
dis[i]=map[1][i];
}
vst[1]=1;
for(i=2;i<=c;i++)
{
k=1;
mindis=INT_MAX;
for(j=2;j<=c;j++)
if(!vst[j]&&dis[j]<mindis)
{
k=j;
mindis=dis[j];
}
if(mindis<INT_MAX)
{
sum+=mindis;
vst[k]=1;
}
else
return INT_MAX;
for(j=2;j<=c;j++)
if(map[k][j]<dis[j]&&!vst[j])
dis[j]=map[k][j];
}
return sum;
}

int main()
{
int i,j;
double temp;
while(~scanf("%d",&t))
{
while(t--)
{
scanf("%d",&c);
for(i=1;i<=c;i++)
scanf("%lf%lf",&island[i].x,&island[i].y);
for(i=1;i<=c;i++)
for(j=i;j<=c;j++)
{
temp=fun_dis(island[i].x,island[i].y,island[j].x,island[j].y);
if(temp>1000||temp<10)
map[i][j]=map[j][i]=INT_MAX;
else
map[i][j]=map[j][i]=temp;
}
temp=prim();
if(temp<INT_MAX)
printf("%.1lf\n",temp*100);
else
printf("oh!\n");
}
}
return 0;
}