畅通工程再续
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2267 Accepted Submission(s): 630
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
Author
8600
解题:
最小生成树,用的克鲁斯卡尔算法。过程中,一个错误改了好久。100我是最后才乘的,造成最后的精度不正确。后来发现,改在中间就乘以100,才AC了。队友很多用普里姆算法,都是在最后乘以100。我的代码如下:
#include <iostream>
#include <iomanip>
#include <stdlib.h>
#include <math.h>
using namespace std;
#define MaxSize 101
#define INF 0x7fffffff
void toEdge(int a[][2],int n,double b[][MaxSize])
{
int i,j;
for (i=0;i<n;i++)
for(j=i+1;j<n;j++)
{
b[i][j]=sqrt((double)((a[j][0]-a[i][0])*(a[j][0]-a[i][0])+(a[j][1]-a[i][1])*(a[j][1]-a[i][1])));
if (b[i][j]<10 || b[i][j]>1000)
b[i][j]=INF;
}
}
typedef struct node
{
int u; //边的起始顶点
int v; //边的终止顶点
double w; //边的权值
}Edge;
int cmp( const void *a , const void *b )
{
struct node *c = (node *)a;
struct node *d = (node *)b;
if(c->w != d->w) return c->w - d->w;
else return d->u - c->u;
}
void kruskal(double edges[][MaxSize],int n,double& sum)
{
int i,j,u1,v1,sn1,sn2,k;
int vset[MaxSize]; //辅助数组,判定两个顶点是否连通
Edge E[MaxSize*MaxSize]; //存放所有的边
k=0; //E数组的下标从0开始
for (i=0;i<n;i++)
{
for (j=0;j<n;j++)
{
if (edges[i][j]!=0 && edges[i][j]!=INF)
{
E[k].u=i;
E[k].v=j;
E[k].w=edges[i][j]*100; //得在这边就乘以100
k++;
}
}
}
if (k<n-1)
return ;
qsort(E,k,sizeof(E[0]),cmp); //快速排序,按权值从小到大排列
for (i=0;i<n;i++) //初始化辅助数组
{
vset[i]=i;
}
k=1; //生成的边数,最后要刚好为总边数
j=0; //E中的下标
while (k<n)
{
u1=E[j].u;
v1=E[j].v;
sn1=vset[u1];
sn2=vset[v1]; //得到两顶点属于的集合编号
if (sn1!=sn2) //不在同一集合编号内的话,把边加入最小生成树
{
//cout<<"( "<<G.vertex[u1]<<"--->"<<G.vertex[v1]<<" ) "<<E[j].w<<endl;
if(E[j].w>INF-1)
{
sum=0;
break;
}
sum+=E[j].w;
k++;
for (i=0;i<n;i++)
{
if (vset[i]==sn2)
{
vset[i]=sn1;
}
}
}
j++;
}
}
int main()
{
int i,j,t,c,island[MaxSize][2];
double sum=0.0,edges[MaxSize][MaxSize];
//freopen("d://1875.txt","r",stdin);
cin>>t;
while (t--)
{
cin>>c;
sum=0.0;
for(i=0;i<c;i++)
for(j=0;j<c;j++)
edges[i][j]=INF;
for(i=0;i<c;i++)
{
cin>>island[i][0]>>island[i][1];
}
toEdge(island,c,edges);
kruskal(edges,c,sum);
if (sum<0.000001)
cout<<"oh!"<<endl;
else
cout<<fixed<<setprecision(1)<<sum<<endl;
}
return 0;
}