Poj 3164 Command Network

描述
一种文字持久的战争后，对武器的战争终于爆发了littleken之间和KnuthOcean的王国。突然和猛烈攻击KnuthOcean的部队已呈现littleken指挥网络彻底失败。必须马上建立临时网络。littleken订单史努比项目负责。
这种情况对每一个细节，史努比认为，当务之急是使littenken的命令到达每一个断开的节点在破坏网络和决定计划建立一个单向通信网络。节点分布在一个平面上。如果littleken的命令可以直接传送从一个节点到另一个节点B，

电线要沿直线段连接的两个节点建立。因为这是在战时，没有之间的所有节点对CAN线建成。史努比想要计划要求导线的最短的总长度，建设可以很快完成。

输入
输入包含多个测试案例。每个测试案例的第一行包含两个整数n（n≤100），在摧毁网络节点数，和M（M≤104），数量的节点对之间的导线，可建。下一行n包含一对有序的XI和易，

给出节点的笛卡尔坐标。

然后按照mlines每个都包含两个整数i和j之间的1和n（含）意义的线可以结我与节点j为前者对后者的单向的命令传递之间建立。littleken总部都位于节点1。进程到文件结束。
输出
对于每个测试用例，输出一行，其中包含最短的总长度到小数点后两位数的总长度。

题解:

最小树形图模板题，直接跑朱刘算法

朱刘算法：http://blog.csdn.net/wsniyufang/article/details/6747392

这一题没有什么别的

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
struct Node
{
    int u,v;
    double dis;
}edge[20001];
int n,m;
double X[10001],Y[10001],ans,in[10001];
int pre[10001],vis[10001],id[10001];
void add(int u,int v,double d,int num)
{
    edge[num].u=u;
    edge[num].v=v;
    edge[num].dis=d;
}
int main()
{int i,j,u,v,root;
//freopen("canal0.in","r",stdin);
//freopen("canal.out","w",stdout);
    cin>>n>>m;
     for (i=1;i<=n;i++)
     {
         scanf("%lf%lf",&X[i],&Y[i]);
     }
     for (i=1;i<=m;i++)
     {
         scanf("%d%d",&u,&v);
         add(u,v,sqrt((X[u]-X[v])*(X[u]-X[v])+(Y[u]-Y[v])*(Y[u]-Y[v])),i);
     }
     ans=0;root=1;
      while (1)
      {//cout<<root<<endl;
          memset(in,127,sizeof(in));
           for (i=1;i<=m;i++)
           if (edge[i].u!=edge[i].v&&edge[i].dis<in[edge[i].v])
           {
               in[edge[i].v]=edge[i].dis;
               pre[edge[i].v]=edge[i].u;
         }
         in[root]=0;
         int cnt=0;
         memset(id,-1,sizeof(id));
        memset(vis,-1,sizeof(vis));
          for (i=1;i<=n;i++)
          {
              ans+=in[i];
               v=i;
               while (vis[v]!=i&&id[v]==-1&&v!=root)
               {
                   vis[v]=i;
                   v=pre[v];
             }
              if (id[v]==-1&&v!=root)
              {
                  id[v]=++cnt;
                   for (j=pre[v];j!=v;j=pre[j]) id[j]=cnt;
              }
          }
           if (cnt==0) break;
           for (i=1;i<=n;i++) 
           if (id[i]==-1) id[i]=++cnt;
            for (i=1;i<=m;i++)
            {
                int v=edge[i].v;
                edge[i].u=id[edge[i].u];
                edge[i].v=id[edge[i].v];
                if (edge[i].u!=edge[i].v) edge[i].dis=edge[i].dis-in[v];
            }
        n=cnt;
        root=id[root];
      }
    printf("%.2lf",ans);
}