POJ 3164 Command Network(最小树形图)

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

本文链接:http://www.cnblogs.com/Ash-ly/p/5532260.html

题意:

　　有一个指挥网络被战争破坏了,你的任务是建立一个临时的通信网络.这个网络最重要的的一点是能够把命令传到被摧毁的网络的每个点上,所以要建立一个单向的传输网络.假设所有的传输节点都分别分布在一个平面上.如果命令要想从节点A传送到结点B上,必须建立一条单向电缆从节点Ａ连接到结点Ｂ.因为是战争时期,不能任意两个节点之间都要连接,所以让你给出满足要求的通信网络所需要的最少的电缆的长度.

思路:

　　电缆的传输信息是单向的,并且命令只要能够传下去就可以,不需要反馈回来.所以可以建立有向图,很显然求得是这个有向图上的最小树形图.这里直接用朱刘算法求最小树形图了.这道题WA了好几次,最后看到一句话"用double时,输入要用%lf输出要用%f".......这题有毒~

代码:

#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <string>
#include <queue>
#include <stack>
#include <vector>
#include <map>
typedef long long LL;
using namespace std;
const int MAXV = 100;
const int MAXE = 10000;
const int INF = 0x3f3f3f3f;

struct Point{
    int x;
    int y;
}pt[MAXV + 7];

double zhuliu(int root, int V, double map[MAXV + 7][MAXV + 7]){
    bool visited[MAXV + 7];
    bool flag[MAXV + 7];
    int pre[MAXV + 7];
    double sum = 0;
    int i, j, k;
    for(i = 0; i <= V; i++) flag[i] = false, map[i][i] = INF;
    pre[root] = root;
    while(true){
        for(i = 1; i <= V; i++){
            if(flag[i] || i == root) continue;
            pre[i] = i;
            for(j = 1; j <= V; j++)
                if(!flag[j] && map[j][i] < map[pre[i]][i])
                    pre[i] = j;
            if(pre[i] == i) return -1;
        }
        for(i = 1; i <= V; i++){
            if(flag[i] || i == root) continue;
            for(j = 1; j <= V; j++) visited[j] = false;
            visited[root] = true;
            j = i;
            do{
                visited[j] = true;
                j = pre[j];
            }while(!visited[j]);
            if(j != i)continue;
            i = j;
            do{
                sum += map[pre[j]][j];
                j = pre[j];
            }while(j != i);
            j = i;
            do{
                for(k = 1; k <= V; k++)
                    if(!flag[k] && map[k][j] < INF && k != pre[j])
                        map[k][j] -= map[pre[j]][j];
                j = pre[j];
            }while(j != i);
            for(j = 1; j <= V; j++){
                if(j == i) continue;
                for(k = pre[i]; k != i; k = pre[k]){
                    if(map[k][j] < map[i][j]) map[i][j] = map[k][j];
                    if(map[j][k] < map[j][i]) map[j][i] = map[j][k];
                }
            }
            for(j = pre[i]; j != i; j = pre[j]) flag[j] = true;
            break;
        }
        if(i > V){
            for(i = 1; i <= V; i++)
                if(!flag[i] && i != root) sum += map[pre[i]][i];
            break;
        }
    }
    return sum;
}

double dist(struct Point a, struct Point b)
{
    return hypot(a.x - b.x, a.y - b.y);
}

int main()
{
    //freopen("input.txt", "r", stdin);
    int n, m;
    while(~scanf("%d%d", &n, &m)){
        double map[MAXV + 7][MAXV + 7] = {0};
        for(int i = 0; i <= n; i++) for(int j = 0; j <= n; j++) map[i][j] = INF;
        memset(&pt, 0, sizeof(Point));
        for(int i = 1; i <= n; i++) scanf("%d%d", &pt[i].x, &pt[i].y);
        for(int i = 0; i < m; i++){
            int u, v;
            scanf("%d%d", &u, &v);
            map[u][v] = dist(pt[u], pt[v]);
        }
        double ans = zhuliu(1, n, map);
        if(ans != -1)    printf("%.2f
", ans);
        else printf("poor snoopy
");
    }
    return 0;
}