BZOJ1626 修建道路【例题精讲】
Description
Farmer John最近得到了一些新的农场,他想新修一些道路使得他的所有农场可以经过原有的或是新修的道路互达(也就是说,从任一个农场都可以经过一些首尾相连道路到达剩下的所有农场)。有些农场之间原本就有道路相连。 所有N(1 <= N <= 1,000)个农场(用1..N顺次编号)在地图上都表示为坐标为(X_i, Y_i)的点(0 <= X_i <= 1,000,000;0 <= Y_i <= 1,000,000),两个农场间道路的长度自然就是代表它们的点之间的距离。现在Farmer John也告诉了你农场间原有的M(1 <= M <= 1,000)条路分别连接了哪两个农场,他希望你计算一下,为了使得所有农场连通,他所需建造道路的最小总长是多少。
Input
* 第1行: 2个用空格隔开的整数:N 和 M
* 第2..N+1行: 第i+1行为2个用空格隔开的整数:X_i、Y_i * 第N+2..N+M+2行: 每行用2个以空格隔开的整数i、j描述了一条已有的道路, 这条道路连接了农场i和农场j
Output
* 第1行: 输出使所有农场连通所需建设道路的最小总长,保留2位小数,不必做 任何额外的取整操作。为了避免精度误差,计算农场间距离及答案时 请使用64位实型变量
Sample Input
4 1
1 1
3 1
2 3
4 3
1 4
输入说明:
FJ一共有4个坐标分别为(1,1),(3,1),(2,3),(4,3)的农场。农场1和农场
4之间原本就有道路相连。
1 1
3 1
2 3
4 3
1 4
输入说明:
FJ一共有4个坐标分别为(1,1),(3,1),(2,3),(4,3)的农场。农场1和农场
4之间原本就有道路相连。
Sample Output
4.00
输出说明:
FJ选择在农场1和农场2间建一条长度为2.00的道路,在农场3和农场4间建一
条长度为2.00的道路。这样,所建道路的总长为4.00,并且这是所有方案中道路
总长最小的一种。
输出说明:
FJ选择在农场1和农场2间建一条长度为2.00的道路,在农场3和农场4间建一
条长度为2.00的道路。这样,所建道路的总长为4.00,并且这是所有方案中道路
总长最小的一种。
此题是一道比较裸的最小生成树,用kruskal就好,只不过题目中有一个干扰项,就是可能会给出一部分已经修建好的道路,然后呢我们可以把所有给了边的两个点加到kruskal的并查集里合并一下,这样在扫描的时候就可以不管这两个点;
但是呢,我一开始就是这样写的,各种地方都修改过了,题解也是这么写的,但是我死活wa五组数据!
迫不得已换招!如果两个点之间给了边,那么就在这两点之间添加一条权值为零的边,然后直接跑最小生成树就好,因为权值为零嘛,一定会加入最小生成树的!
下面上代码:
1 //苍山负雪,明烛天南; 2 #include<iostream> 3 #include<cstdio> 4 #include<cstdlib> 5 #include<cstring> 6 #include<cmath> 7 #include<string> 8 #include<queue> 9 #include<algorithm> 10 using namespace std; 11 const int MAXM=1000005; 12 const int MAXN=1005; 13 struct node{ 14 int x,y; 15 double z; 16 }edge[MAXM]; 17 bool cmp(node a,node b){ 18 return a.z<b.z; 19 } 20 int father[MAXN],n,m,k; 21 int tot=0; 22 long long xx[MAXN],yy[MAXN]; 23 double dist(int i,int j) 24 { 25 long long fuck1=xx[i]-xx[j]; 26 long long fuck2=yy[i]-yy[j]; 27 return sqrt(fuck1*fuck1*1.0+fuck2*fuck2*1.0); 28 } 29 void add(int x,int y,double z) 30 { 31 edge[++tot].x=x; 32 edge[tot].y=y; 33 edge[tot].z=z; 34 } 35 int get(int x) 36 { 37 if(father[x]==x)return x; 38 return father[x]=get(father[x]); 39 } 40 int main() 41 { 42 //freopen("2.in","r",stdin); 43 //freopen("2.out","w",stdout); 44 double ans=0.0; 45 scanf("%d%d",&n,&m); 46 for(int i=1;i<=n;i++){ 47 scanf("%lld%lld",&xx[i],&yy[i]); 48 father[i]=i; 49 } 50 for(int i=1;i<=m;i++){ 51 int x,y; 52 scanf("%d%d",&x,&y); 53 add(x,y,0); 54 } 55 for(int i=1;i<=n;i++){ 56 for(int j=i+1;j<=n;j++){ 57 if(i!=j){ 58 int x=get(i),y=get(j); 59 if(x!=y){ 60 add(i,j,dist(i,j)); 61 } 62 } 63 } 64 } 65 sort(edge+1,edge+1+tot,cmp); 66 ans=0.0; 67 for(int i=1;i<=tot;i++){ 68 int x=get(edge[i].x); 69 int y=get(edge[i].y); 70 if(x==y)continue; 71 ans=ans+edge[i].z; 72 father[x]=y; 73 } 74 printf("%0.2lf",ans); 75 puts(""); 76 //fclose(stdin); 77 //fclose(stdout); 78 return 0; 79 }

浙公网安备 33010602011771号