HDOJ 1875 畅通工程再续 (最小生成树)

畅通工程再续

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

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

Sample Output

1414.2
oh!

Author

8600

Source

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

先保存所有坐标然后对于任意两个不同的坐标求距离，如果距离不符合要求则将它们距离置于无限大inf，因为边数目多，用prim算法

</pre><pre name="code" class="cpp">#include<iostream>
#include<vector>
#include<cstdio>
#define MAXN 1005
#include<cmath>
using namespace std;
int m,n;
const double inf=10000;
double point[MAXN][2];
vector<int>mmap[MAXN];
double val[MAXN][MAXN];
int vis[MAXN],t;
double  prim()
{
    double  re=0;
    int k=0;
    double d[MAXN];
    fill(d,d+n+1,inf);
    fill(vis,vis+n+1,0);
    d[1]=0;
    while(true)
    {
        int v=-1;
        for(int i=1;i<=n;i++){
            if(!vis[i]&&(v==-1||d[i]<d[v]))
                v=i;
        }
        if(v==-1||d[v]==inf)
                return -1;
        re+=d[v];
        k++;
        if(k==n)return re;
        vis[v]=1;
        for(int i=0;i<mmap[v].size();i++)
        {
            int x=mmap[v][i];
            if(!vis[x]&&val[v][x]<d[x])
                d[x]=val[v][x];
        }

    }
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        m=0;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
         {
             mmap[i].clear();
             scanf("%lf%lf",&point[i][0],&point[i][1]);
         }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<i;j++)
            {
                m++;
                double x1=point[i][0],y1=point[i][1];
                double x2=point[j][0],y2=point[j][1];
                mmap[i].push_back(j);
                mmap[j].push_back(i);
                double len=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
                if(len<10||len>1000)
                    len=inf;
                val[i][j]=val[j][i]=len;

            }
        }
       double re=prim();
       if(re==-1)
        printf("on!
");
       else printf("%.1lf
",re*100);
    }
}