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
2 2 10 10 20 20 3 1 1 2 2 1000 1000
 

Sample Output
1414.2 oh!
 

Author
8600
 

Source

 

先保存所有坐标然后对于任意两个不同的坐标求距离,如果距离不符合要求则将它们距离置于无限大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!\n");
       else printf("%.1lf\n",re*100);
    }
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

posted @ 2015-08-27 09:05  Thereisnospon  阅读(141)  评论(0编辑  收藏  举报