POJ 3164最小树形图

/*POJ 3164
一个通信网络被破坏了
给定每个点的坐标，和两两之间可修复的边和长度
现在选择部分修通，使从1点出发可到达所有点，并且代价最短

可看出是一个生成树的问题，做了UVA的那道，这道很简单了

PS:unidirection 是有向
   undirecttion 才是无向
C++交才可以
*/
#include<iostream>
#include<string.h>
#include<stdio.h>
#include<stdlib.h>
#include<cmath>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#define maxn 100+5
#define maxe 20000+5
#define typec double
using namespace std;

const typec inf=99999999;
int nextint(){int x;scanf("%d",&x);return x;}
int nextdou(){double x;scanf("%lf",&x);return x;}
double X[maxn],Y[maxn];
double Dis(int i,int j){
    double x1=X[i],x2=X[j],y1=Y[i],y2=Y[j];
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
struct Edge{
    int u,v;
    typec c;//权值
}edge[maxe];

int vis[maxn],pre[maxn],id[maxn];
typec in[maxn];
typec zhuliu(int root,int nv,int ne){//最小树形图
      int i,j,u,v;
      typec ans=0;
      while(1){
              for(i=0;i<nv;i++)in[i]=inf;
              for(i=0;i<ne;i++){
                      u=edge[i].u;
                      v=edge[i].v;
                      if(u!=v&&edge[i].c<in[v]){
                              pre[v]=u;
                              in[v]=edge[i].c;
                      }
              }
              in[root]=0;
              for(i=0;i<nv;i++)if(fabs(in[i]-inf)<1e-6)return -1;
              int cnt=0;
              memset(id,-1,sizeof(id));
              memset(vis,-1,sizeof(vis));
              for(i=0;i<nv;i++){
                      v=i;
                      ans+=in[i];
                      while(vis[v]!=i&&id[v]==-1&&v!=root){
                              vis[v]=i;
                              v=pre[v];
                      }
                      if(v!=root&&id[v]==-1){
                              for(u=pre[v];u!=v;u=pre[u])id[u]=cnt;
                              id[v]=cnt++;
                      }
              }
              if(cnt==0)break;
              for(i=0;i<nv;i++)if(id[i]==-1)id[i]=cnt++;
              for(i=0;i<ne;i++){
                      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].c-=in[v];
              }
              nv=cnt;root=id[root];
      }
        return ans;
}
int n,m;
int cnt;
void built_tree(){
    cnt=0;
    for(int i=0;i<m;i++){
        int u=nextint(),v=nextint();
        u--,v--;
        double dis=Dis(u,v);
//        cout<<u<<","<<v<<" dis="<<dis<<endl;
        edge[cnt].u=u;
        edge[cnt].v=v;
        edge[cnt++].c=dis;
//        edge[cnt++]=(Edge){v,u,dis};
    }
}
int main(){
//    freopen("out.txt","w",stdout);
    while(cin>>n>>m){
        for (int i=0;i<n;i++){
            X[i]=nextdou();
            Y[i]=nextdou();
        }
        built_tree();
        double ans=zhuliu(0,n,cnt);
        if (ans<0) printf("poor snoopy
");else printf("%.2lf
",ans);
    }
    return 0;
}