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

HDU - 2121 ：http://acm.hdu.edu.cn/showproblem.php?pid=2121

比较好的朱刘算法blog：https://blog.csdn.net/txl199106/article/details/62045479

题意：

　　在一个有向图中，找一个点，使得这个点到其他点的距离和最小，输出距离和，和这个点的坐标。

思路：

　　无根最小树形图，设所有的有向图的距离和为sum。自己建立一个虚拟的原点（n+1），向每一个节点连一条距离为sum+1的边。以n+1为根结点跑一遍最小树形图（复杂度O（VE）），如果求出的ans == -1 或者 ans >=2*（sum + 1），无解，因为这么大的ans，只可能用了两条我们自己建立的边。由于每条边的端点在跑最小树形图的时候会改变，所以记录这是第几条边rtt，结果就是rtt - m，代码中由减了1是因为原图是Base0的。

/*
* @Author: chenkexing
* @Date:   2018-09-05 11:05:14
* @Last Modified by:   chenkexing
* @Last Modified time: 2018-09-10 20:21:22
*/
#include <algorithm>
#include  <iterator>
#include  <iostream>
#include   <cstring>
#include   <cstdlib>
#include   <iomanip>
#include    <bitset>
#include    <cctype>
#include    <cstdio>
#include    <string>
#include    <vector>
#include     <stack>
#include     <cmath>
#include     <queue>
#include      <list>
#include       <map>
#include       <set>
#include   <cassert>
using namespace std;
//#pragma GCC optimize(3)
//#pragma comment(linker, "/STACK:102400000,102400000")  //c++
#define lson (l , mid , rt << 1)
#define rson (mid + 1 , r , rt << 1 | 1)
#define debug(x) cerr << #x << " = " << x << "\n";
#define pb push_back
#define pq priority_queue

typedef long long ll;
typedef unsigned long long ull;

typedef pair<ll ,ll > pll;
typedef pair<int ,int > pii;
typedef pair<int,pii> p3;

//priority_queue<int> q;//这是一个大根堆q
//priority_queue<int,vector<int>,greater<int> >q;//这是一个小根堆q
#define fi first
#define se second
//#define endl '\n'

#define OKC ios::sync_with_stdio(false);cin.tie(0)
#define FT(A,B,C) for(int A=B;A <= C;++A)  //用来压行
#define REP(i , j , k)  for(int i = j ; i <  k ; ++i)
//priority_queue<int ,vector<int>, greater<int> >que;

const ll mos = 0x7FFFFFFF;  //
const ll nmos = 0x80000000;  //-2147483648
const int inf = 0x3f3f3f3f;
const ll inff = 0x3f3f3f3f3f3f3f3f; //
const int mod = 1e9+;
const double esp = 1e-;
const double PI=acos(-1.0);

template<typename T>
inline T read(T&x){
    x=;int f=;char ch=getchar();
    while (ch<''||ch>'') f|=(ch=='-'),ch=getchar();
    while (ch>=''&&ch<='') x=x*+ch-'',ch=getchar();
    return x=f?-x:x;
}

/*-----------------------showtime----------------------*/
            const int maxm = ;
            const int maxn = ;
            struct Edge
            {
                int from,to,c;
            }e[maxm];
            int in[maxn],vis[maxn],pre[maxn],id[maxn];
            int rtt;
            int zhuliu(int root,int n,int m){
                int res = ;
                while(true){
                    memset(in, inf, sizeof(in));
                    for(int i=; i<=m; i++){
                        if(e[i].from != e[i].to && e[i].c < in[e[i].to]){
                            pre[e[i].to] = e[i].from;
                            in[e[i].to] = e[i].c;
                            if(e[i].from == root)rtt = i;
                        }
                    }
                    for(int i=; i<=n; i++){
                        if(i!=root&&in[i] == inf)
                            return -;
                    }
                    int tn = ,v;
                    memset(id,-,sizeof(id));
                    memset(vis,-,sizeof(vis));

                    in[root] = ;
                    for(int i=; i<=n; i++){
                        res += in[i];
                        v = i;
                        while(v !=root && id[v] == - && vis[v] != i){
                            vis[v] = i;
                            v = pre[v];
                        }
                        if(v!=root && id[v] == -){
                            id[v] = ++tn;
                            for(int u = pre[v]; u!=v; u = pre[u]){
                                id[u] = tn;
                            }
                        }
                    }
                    if(tn == )break;
                    for(int i=; i<=n; i++){
                        if(id[i] == -)id[i] = ++tn;
                    }

                    for(int i=; i<=m; i++){
                        int v = e[i].to;
                        e[i].to = id[e[i].to];
                        e[i].from = id[e[i].from];
                        if(e[i].to != e[i].from){
                            e[i].c -= in[v];
                        }
                    }

                    n = tn;root = id[root];
                }
                return res;

            }
int main(){
            int n,m,r;
            while(~scanf("%d%d", &n, &m))
            {
                    int sum = ;
                    for(int i=; i<=m; i++){
                        int u,v,c;
                        scanf("%d%d%d", &u, &v, &c);
                        u++,v++;
                        e[i].from = u;e[i].to = v;
                        e[i].c = c;
                        sum += c;
                    }
                    sum++;
                    for(int i=m+; i<=n+m; i++){
                        e[i].from = n+;
                        e[i].to = i-m;
                        e[i].c = sum;
                    }

                    int ans = zhuliu(n+,n+,n+m);
                    // debug(ans);
                    if(ans == - || ans - sum >= sum){
                        puts("impossible");
                    }
                    else {
                        printf("%d %d\n", ans - sum, rtt - m - );
                    }
                    printf("\n");
            }

            return ;
}

HDU - 2121