（最小生成树） 畅通工程再续 -- HDU --1875

链接：

http://acm.hdu.edu.cn/showproblem.php?pid=1875

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=82831#problem/N

畅通工程再续

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

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<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<math.h>

#define N 210
#define INF 0xfffffff
double  G[N][N],dis[N];
int vis[N];

struct node
{
    int x,y;
}p[N];

void Inn()
{
    int i,j;
    for(i=;i<=N;i++)
    {
        dis[i]=INF;
        for(j=;j<=N;j++)
        {
            G[i][j]=INF;
        }
    }
    memset(vis,,sizeof(vis));
}

double prim(int s,int n)
{
    double ans=;
    int i,j;

    for(i=;i<=n;i++)
        dis[i]=G[s][i];

    memset(vis, , sizeof(vis));
    vis[s]=;

    for(i=;i<n;i++)
    {
        int dist=;
        double Min=INF;
        for(j=;j<=n;j++)
        {
            if(!vis[j] && dis[j]<Min)
            {
                dist=j;
                Min=dis[j];
            }
        }
        vis[dist]=;
        ans+=Min;
        for(j=;j<=n;j++)
        {
            if(!vis[j] && dis[j]>G[dist][j])
            {
                dis[j]=G[dist][j];
            }
        }
    }
    return ans;
}
int main()
{
    int T,n,i,j;
    scanf("%d",&T);

    while(T--)
    {
        Inn();

        scanf("%d",&n);
        memset(p,,sizeof(p));
        for(i=;i<=n;i++)
            scanf("%d %d",&p[i].x,&p[i].y);

        for(i=;i<=n;i++)
        {
            for(j=;j<=i;j++)
            {
                double d=sqrt(1.0*(p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y));
                if(d>= && d<=)
                    G[i][j]=G[j][i]=d;
            }
        }
        double aa=prim(,n);
        if(aa>INF)
            printf("oh!\n");
        else
            printf("%.1f\n",aa*);
    }
    return ;
}