Hdu 2121 Ice_cream’s world II 无根最小树形图

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2121

题目大意：给你N个城市，选择一个城市建首都，城市编号为0~N-1，给你M条路，每条路包括u，v，w，即从u到v的花费为w，现在问你那个城市建首都的话，从这个城市到其余城市的花费最小，最小花费是多少，如果有多个适合的城市，输出编号小的城市。无法到达输出impossible。

解题思路：这是无根最小树形图的题，我们可以先假定一个虚根，即源点，并设源点到每个城市的权值为   全图所有权值之和大一点，然后就是树形图模板，细节在代码中。

代码如下：

//本题边长都为__int64,其实int型也可以过，但长整型应该更广，所以本题都为__int64
#include <cstdio>
#include<iostream>
#include <cstring>
#include <cstdlib>
#include<cmath>
using namespace std;
typedef double type;
#define INF 2000000000
#define N 1005

struct edge //图的结构体
{
    int u,v;__int64 w;
}e[10005];

int m,n,pre[N],id[N],visit[N],xroot;
__int64 in[N],sum;
//eCnt为图中的边数
//n为图中的顶点数
//pre[i]为顶点i的前驱节点
//id[i]为缩环，形成新图的中间量
//in[i]为点i的最小入边
//visit[i]遍历图时记录顶点是否被访问过
__int64 directedMST(int root,int nv,int ne)
{
    __int64 ans=0;
    while(1)
    {
        //1.找最小入边
        for(int i=0;i<nv;i++) in[i]=INF;
        for(int i=0;i<ne;i++)
        {
            int u=e[i].u;
            int v=e[i].v;
            if(u!=v&&e[i].w<in[v])
            {
                if(u==root)  //此处标记与源点相连的最小边
                    xroot=i;
                in[v]=e[i].w;
                pre[v]=u;
            }
        }
        for(int i=0;i<nv;i++)                          //判断图是否连通
            if(i!=root&&in[i]==INF) return -1;  //除了跟以外有点没有入边,则根无法到达它
         //2.找环
        int nodeCnt=0;          //图中环的数目
        memset(id, -1, sizeof(id));
        memset(visit, -1, sizeof(visit));
        in[root]=0;
        for(int i=0;i<nv;i++)
        {
            ans+=in[i];
            int v=i;
            while(visit[v]!=i&&id[v]==-1&&v!=root)//每个点寻找其前序点，要么最终寻找至根部，要么找到一个环
            {
                visit[v]=i;
                v=pre[v];
            }
            if(v!=root&&id[v]==-1)//缩点
            {
                for(int u=pre[v];u!=v;u=pre[u])
                    id[u]=nodeCnt;
                id[v]=nodeCnt++;
            }
        }
        if(nodeCnt==0) break;//如果无环，跳出循环
        for(int i=0; i<nv; i++)
            if(id[i]==-1)
                id[i]=nodeCnt++;
        //3.缩点，重新标记
        for(int i=0;i<ne;i++)
        {
            int v=e[i].v;
            e[i].u=id[e[i].u];
            e[i].v=id[e[i].v];
            if(e[i].u!=e[i].v)
                e[i].w-=in[v];
        }
        nv=nodeCnt;
        root=id[root];
    }
    return ans;
}

int main()
{
    int m;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        sum=0;
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%I64d",&e[i].u,&e[i].v,&e[i].w);
            e[i].u++;e[i].v++;     //都++之后，把0设为超级源点，联通各点
            sum+=e[i].w;
            if(e[i].u==e[i].v)
                e[i].w=INF;//消除自环
        }
        sum++;      //此处必须++，因为需要权值比总权值大，因为这个w几次，，，
        for(int i=m;i<n+m;i++)
        {
            e[i].u=0;
            e[i].v=i-m+1;
            e[i].w=sum;
        }
        __int64 ans=directedMST(0,n+1,m+n);
        if(ans==-1 || ans-sum>=sum) printf("impossible\n");//ans-sum是除去虚根的最小树形图的最短路径，如果这个距离比所有的边权值和sum还大，说明还有另外的边由虚点发出，故说明此图不连通
        else printf("%I64d %d\n", ans-sum, xroot-m);
        printf("\n");
    }
    return 0;
}