151. [USACO Dec07] 建造路径

★★   输入文件：roads.in   输出文件：roads.out   简单对比
时间限制：1 s  
内存限制：128 MB

译 by CmYkRgB123

描述

Farmer John 刚刚得到了几个新农场！他想把这几个农场用路连接起来，这样他就可以通过笔直的公路从一个农场到另一个农场了。现在已经有了几条连接着的农场。

N (1 ≤ N ≤ 1,000) 个农场中，每个农场的位置在坐标平面的 (Xi, Yi) (0 ≤ Xi ≤ 1,000,000; 0 ≤
Yi ≤ 1,000,000)。已经有 M (1 ≤ M ≤ 1,000) 条路以前就被建好了。请你帮助 Farmer John
考虑建设尽量少长度的额外的路，使他的农场连在一起。

输入

* 第 1 行: 两个整数: N , M
* 第 2..N+1 行: 两个整数 Xi , Yi
* 第 N+2..N+M+2 行: 两个整数: i , j, 表示已经存在从农场i到农场j的路。

输出

* 第 1 行: 额外的路的最少长度，保留2小数。 请使用 64 位的浮点数。

样例输入

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

样例输出

4.00

好久没写最小生成树了结果爆了一堆bug。。

对于已经建造的道路，我们可以把它的权值设置成0

然后跑裸地kruskal

注意内存限制！！！！！！！！！！！

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<algorithm>
 #define INF 0x7ffff
 using namespace std;
 const int MAXN=;
 int vis[MAXN][MAXN];// 记录两个城市之间是否已经有道路相连
 double dis[MAXN][MAXN];
 struct node
 {
     int u,v,nxt;
     double w;
 }edge[];
 int f[MAXN*];
 int num=;
 struct pos
 {
     long long int  x,y;
 }where[MAXN*];
 int n,m;
 double ans=;
 int read(int & n)
 {
     int flag=,x=;char c='/';
     while(c<''||c>''){c=getchar();if(c=='-')flag=;}
     while(c>=''&&c<='')x=x*+(c-),c=getchar();
     if(flag)n=-x;
     else n=x;
 }
 void deal_dis()
 {
     //for(int i=1;i<=n;i++)
     //    for(int j=1;j<i;j++)
     //        dis[i][j]=INF;

     for(int i=;i<=n;i++)
         for(int j=;j<i;j++)
             dis[i][j]=sqrt((where[i].x-where[j].x)*(where[i].x-where[j].x)+(where[i].y-where[j].y)*(where[i].y-where[j].y)); 

     /*for(int i=1;i<=n;i++)
     {
         for(int j=1;j<i;j++)
             printf("%.2lf ",dis[i][j]);
         printf("\n");
     }*/
     //cout<<dis[219][65];
 }
 void add_edge(int p,int q,double we)
 {
     edge[num].u=p;
     edge[num].v=q;
     edge[num].w=we;
 //    cout<<edge[num].w<<endl;
     num++;
 }
 int find(int x)
 {
     if(f[x]!=x)
     f[x]=find(f[x]);
     return f[x];
 }
 int unionn(int x,int y)
 {
     int fx=find(x);
     int fy=find(y);
     f[fx]=fy;
 }
 int comp(const node & a,const node & b)
 {
     return a.w<b.w;
 }
 void kruskal()
 {
     sort(edge+,edge+num+,comp);
     int k=;
     for(int i=;i<=num;i++)
     {
         if(find(edge[i].u)!=find(edge[i].v))
         {
             unionn(edge[i].u,edge[i].v);
             k++;
             ans=ans+edge[i].w;
         }
         if(k==n)
         break;
     }
     printf("%.2lf",ans);
 }
 int main()
 {
     freopen("roads.in","r",stdin);
     freopen("roads.out","w",stdout);
     read(n);read(m);
     for(int i=;i<=n;i++)
         cin>>where[i].x>>where[i].y,f[i]=i;

     for(int i=;i<=m;i++)
     {
         int x,y;
         read(x);read(y);
         add_edge(x,y,0.00);
         vis[x][y]=;
     }

     deal_dis();

     for(int i=;i<=n;i++)
         for(int j=;j<i;j++)
             if(i!=j&&vis[i][j]==)
             add_edge(i,j,dis[i][j]);

     kruskal();
     return ;
 }