poj 3164 Command Network

http://poj.org/problem?id=3164

比着别人的代码写的 有些地方还是不理解 但是不能老在这上面耽误 先放一下 以后遇到再慢慢研究

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<map>
#include<queue>
#include<cmath>
#define LL long long

using namespace std;

const int N=105;
const double M=pow(2.0,32);
struct point
{
    double x,y;
}mem[N];
struct ss
{
    int i,j;
    double d;
}side[N*N];
bool lin[N][N];
double in[N];
int vi[N];
int pre[N];
int f[N];
int n,m;
double Fdist(int i,int j)
{
    return sqrt((mem[i].x-mem[j].x)*(mem[i].x-mem[j].x)+
                (mem[i].y-mem[j].y)*(mem[i].y-mem[j].y));
}
double min_utree()
{
    double ans=0.0;
    int v=n;
    int root=1;
    while(1)
    {
        for(int i=1;i<=v;++i)
        in[i]=M;
        for(int l=0;l<m;++l)
        {
            int x=side[l].i;
            int y=side[l].j;
            if(x==y)
            continue;
            if(side[l].d<in[y])
            {
                pre[y]=x;
                in[y]=side[l].d;
            }
        }
        memset(f,-1,sizeof(f));
        memset(vi,-1,sizeof(vi));
        int num=0;
        in[root]=0;
        for(int i=1;i<=v;++i)
        {
            ans=ans+in[i];
            int l=i;
            while(l!=root&&vi[l]!=i&&f[l]==-1)
            {
                vi[l]=i;
                l=pre[l];
            }
            if(l!=root&&f[l]==-1)
            {
                ++num;
                int k=pre[l];
                while(k!=l)
                {
                    f[k]=num;
                    k=pre[k];
                }
                f[l]=num;
            }
        }
        if(num==0)
        break;
        for(int i=1;i<=v;++i)
        if(f[i]==-1)
        f[i]=(++num);
        for(int l=0;l<m;++l)
        {
            int x=side[l].i;
            int y=side[l].j;
            side[l].i=f[x];
            side[l].j=f[y];
            if(f[x]!=f[y])
            {
                side[l].d=(side[l].d-in[y]);
            }
        }
        root=f[root];
        v=num;
    }
    return ans;
}
void dfs(int x)
{
    vi[x]=1;
    for(int i=1;i<=n;++i)
    {
        if(vi[i]==0&&lin[x][i])
        dfs(i);
    }
}
int main()
{
    //freopen("data.txt","r",stdin);
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        for(int i=1;i<=n;++i)
        scanf("%lf %lf",&mem[i].x,&mem[i].y);
        memset(lin,false,sizeof(lin));
        for(int l=0;l<m;++l)
        {
            scanf("%d %d",&side[l].i,&side[l].j);
            side[l].d=Fdist(side[l].i,side[l].j);
            if(side[l].i==side[l].j)
            side[l].d=M;
            lin[side[l].i][side[l].j]=true;
        }
        memset(vi,0,sizeof(vi));
        dfs(1);
        int l;
        for(l=1;l<=n;++l)
        if(vi[l]==0)
        break;
        if(l<=n)
        {
            printf("poor snoopy\n");
            continue;
        }
        double ans=min_utree();
        printf("%.2f\n",ans);
    }
    return 0;
}