最小生成树———— 畅通工程再续

N - 畅通工程再续

 

相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其他的小岛时都要通过划小船来实现。现在政府决定大力发展百岛湖，发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！经过考察小组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!

思路：按照最小生成树的思想，找到一棵树。本题我用的是prim算法，欢迎指正，代码中有注释。

代码：

#include<string.h>
#include<stdio.h>
#include<math.h>
#define inf 0x3f3f3f3f
double dis[105],e[105][105];
int x[105],y[105];
int vis[105],n;
double prim()
{
    memset(vis,0,sizeof(vis));
    double ans=0;
    for(int i=1;i<=n;i++)
        e[i][i]=0;
    for(int i=2;i<=n;i++)
       dis[i]=e[1][i];
    vis[1]=1;
    for(int i=2;i<=n;i++)
    {
        int u=-1;
        double minn=inf;
        for(int j=2;j<=n;j++)
        {
            if(!vis[j]&&dis[j]<minn)//找到没有被标记的点作为试探点（最近的那个），每次试探得到一个最小值，赋给dis
            {
                u=j;
                minn=dis[j];
            }
        }
        if(u==-1) return -1;
        vis[u]=1;
        ans+=minn;
        for(int v=1;v<=n;v++)
        {
            if(!vis[v]&&dis[v]>e[u][v])//找到最佳试探点之后，以该点为中心，再找下一个点（注意：与dij算法不同之处在于，dij是if(dis[j]>dis[u]+e[u][j]),而这个是if(dis[j]>e[u][j]))
                dis[v]=e[u][v];
        }
    }
    return ans;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        double tmp,tmp_X,tmp_y;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d%d",&x[i],&y[i]);
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(i==j) e[i][j]=0;
                else
                {
                    tmp_X=(x[i]-x[j])*(x[i]-x[j]);
                    tmp_y=(y[i]-y[j])*(y[i]-y[j]);
                    tmp=sqrt(tmp_X+tmp_y);
                    if(tmp<10||tmp>1000) e[i][j]=e[j][i]=inf;
                    else e[i][j]=e[j][i]=tmp;
                }
            }
        }
            double ans=prim();
            if(ans==-1) printf("oh!\n");
            else printf("%.1lf\n",ans*100);
    }
}