【洛谷 2872】道路建设 Building Roads

题目描述

Farmer John had just acquired several new farms! He wants to connect the farms with roads so that he can travel from any farm to any other farm via a sequence of roads; roads already connect some of the farms.

Each of the N (1 ≤ N ≤ 1,000) farms (conveniently numbered 1..N) is represented by a position (Xi, Yi) on the plane (0 ≤ Xi ≤ 1,000,000; 0 ≤ Yi ≤ 1,000,000). Given the preexisting M roads (1 ≤ M ≤ 1,000) as pairs of connected farms, help Farmer John determine the smallest length of additional roads he must build to connect all his farms.

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)条路分别连接了哪两个农场，他希望你计算一下，为了使得所有农场连通，他所需建造道路的最小总长是多少。

输入格式

* Line 1: Two space-separated integers: N and M

* Lines 2..N+1: Two space-separated integers: Xi and Yi

* Lines N+2..N+M+2: Two space-separated integers: i and j, indicating that there is already a road connecting the farm i and farm j.

• 第1行: 2个用空格隔开的整数：N 和 M

• 第2..N+1行: 第i+1行为2个用空格隔开的整数：X_i、Y_i

• 第N+2..N+M+2行: 每行用2个以空格隔开的整数i、j描述了一条已有的道路， 这条道路连接了农场i和农场j

输出格式

* Line 1: Smallest length of additional roads required to connect all farms, printed without rounding to two decimal places. Be sure to calculate distances as 64-bit floating point numbers.

输出使所有农场连通所需建设道路的最小总长，保留2位小数，不必做 任何额外的取整操作。为了避免精度误差，计算农场间距离及答案时 请使用64位实型变量

输入输出样例

输入 #1复制

4 1
1 1
3 1
2 3
4 3
1 4

输出 #1复制

4.00

说明/提示

题目简述：给出n个点的坐标,其中一些点已经连通,现在要把所有点连通,求修路的最小长度.

感谢@睿屿青衫丶 提供翻译

题解 ：几乎是裸的生成树，不过每道题都有坑啊！这道题的坑在于sqrt，相乘容易爆int，

需要强制转换成double，再计算(double 和long long范围差不多，都是18位)

#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<cstdio>
using namespace std;
typedef double db;
const int N=1005;

struct node{
    int u,v;
    db w;
}e[1002019];
int n,m,x[N],y[N],cp,fa[N];

int find(int x){
    if(x!=fa[x]) 
       fa[x]=find(fa[x]);
    return fa[x];
} 

db how_long(int uu,int vv){
    int x5=x[uu]; int x6=x[vv];
    int y5=y[uu]; int y6=y[vv];
    return (db)(sqrt((db)(x5-x6)*(db)(x5-x6)+(db)(y5-y6)*(db)(y5-y6)));
}
void add_YCLL(int aa,int bb){
    cp++; e[cp].u=aa; e[cp].v=bb;
    e[cp].w=how_long(aa,bb);
}
bool cmp(node aa,node bb)
    { return aa.w<bb.w; }
int main(){
    freopen("2872.in","r",stdin);
    freopen("2872.out","w",stdout);
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d %d",&x[i],&y[i]);
    for(int i=1;i<=n;i++) fa[i]=i;
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
            add_YCLL(i,j);
    for(int i=1;i<=m;i++){
        int xx; int yy;
        scanf("%d %d",&xx,&yy);
        add_YCLL(xx,yy); e[cp].w=0.0;
    }
    sort(e+1,e+cp+1,cmp);
    db ans=0.0; 
    //int tot=0;
    for(int i=1;i<=cp;i++){
        int u=find(e[i].u);
        int v=find(e[i].v);
        if(u==v) continue;
        ans+=e[i].w; 
        fa[u]=v; //tot++;
        //if(tot==(n-1)) break;
    } 
    printf("%.2lf",ans);
    return 0;
}