POJ3164 Command Network 最小树形图

　　题目链接：http://poj.org/problem?id=3164

　　裸的最小树形图，用朱—刘算法解决，具体实现过程如下：算法一开始先判断从固定根开始是否可达所有原图中的点，若不可，则一定不存在最小树形图。这一步是一个很随便的搜索，写多搓都行，不加废话。第二步，遍历所有的边，从中找出除根结点外各点的最小入边，累加权值，构成新图。接着判断该图是否存在环。若不存在，则该图便是所求最小树型图，当前权为最小权。否则对环缩点，然后回到第二步继续判断。简化就是三个过程：找边—>找环—>缩点；

　　算法流程图：

                        

　　参考样例：

            

　　代码实现有很多种方法，为了兼顾代码的长度和效率，我选择了O(VE)的写法。

//STATUS:C++_AC_110MS_244KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
#define LL __int64
#define pdi pair<double,int>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int N=110,M=1000000,INF=0x3f3f3f3f,MOD=1999997;
const LL LLNF=0x3f3f3f3f3f3f3f3fLL;
const double DNF=100000000000;

struct Node{
    double x,y;
}nod[N];
struct Edge{
    int u,v;
    double w;
}e[N*N];
int pre[N],id[N],vis[N];
double minw[N];
int n,m;

inline double dist(Node &a,Node &b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double zhu_liu(int root)
{
    int i,cou,u,v,k;
    double ans=0;
    while(1)
    {
        //init
        mem(pre,-1);
        for(i=1;i<=n;i++)minw[i]=DNF;
        for(i=0;i<m;i++){
            u=e[i].u;
            v=e[i].v;
            if(e[i].w<minw[v] && u!=v){
                pre[v]=u;
                minw[v]=e[i].w;
            }
        }
        pre[root]=-1;minw[root]=0;
        for(cou=0,i=1;i<=n;i++)
            if(pre[i]==-1 && i!=root)cou++;
            else ans+=minw[i];
        if(cou)return -1;
        //cheack the circle
        mem(vis,0);
        mem(id,0);
        for(i=1,k=0;i<=n;i++){
            if(id[i])continue;
            u=i;
            while(u!=-1 && !id[u] && vis[u]!=i){
                vis[u]=i;
                u=pre[u];
            }
            if(u!=-1 && !id[u] && vis[u]==i){
                k++;
                while(id[u]!=k){
                    id[u]=k;
                    u=pre[u];
                }
            }
        }
        if(!k)break;
        for(i=1;i<=n;i++)if(!id[i])id[i]=++k;
        //eliminate circle
        for(i=0;i<m;i++){
            e[i].w-=minw[e[i].v];
            e[i].u=id[e[i].u];
            e[i].v=id[e[i].v];
        }
        n=k;
        root=id[root];
    }
    return ans;
}

int main()
{
 //   freopen("in.txt","r",stdin);
    int i,j,a,b;
    double c,ans;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(i=1;i<=n;i++)
            scanf("%lf%lf",&nod[i].x,&nod[i].y);
        for(i=0;i<m;i++){
            scanf("%d%d",&a,&b);
            c=dist(nod[a],nod[b]);
            e[i].u=a,e[i].v=b,e[i].w=c;
        }

        ans=zhu_liu(1);
        if(ans>=0)printf("%.2lf\n",ans);
        else printf("poor snoopy\n");
    }
    return 0;
}