【無定根朱劉演算法】hdu 3072

 解法:

由於沒有根,所以我們可以虛擬一個根,到每個點的權值很大很大(所有權值+1即可),權值很大可以保證最後只有一個點連這個虛擬根,然後最大答案減去這個很大的權值,而要輸出最小的根的話則有點噁心- -由於我的模板會改變點的ID號,所以要儲存原來的ID,如果和最初的跟節點連得話,那麼更新一下最小的ID(也可不用儲存,根據點的儲存順序直接算出來)

#include <list>#include <map>#include <set>#include <queue>#include <string>#include <deque>#include <stack>#include <algorithm>#include <iostream>#include <iomanip>#include <cstdio>#include <math.h>#include <cstdlib>#include <limits.h>#include <time.h>#include <string.h>using namespace std;#define LL long long#define PI acos(-1.0)#define FRE freopen("a.txt","r",stdin)#define N 1005#define type int#define MAX INT_MAXstruct Edge{    int u , v;    type cost;}E[11005];int pre[N],ID[N],vis[N];type In[N];int pos;type zhuliu(int root,int NV,int NE) {    type ret = 0;    while(true) {        //1.找最小入邊        for(int i=0;i<NV;i++) In[i] = MAX;        for(int i=0;i<NE;i++){            int u = E[i].u;            int v = E[i].v;            if(E[i].cost < In[v] && u != v) {  //這一步可以把自環切掉，找出每個點的最小入邊                pre[v] = u;                if(u==root)                pos=i;                In[v] = E[i].cost;            }        }        for(int i=0;i<NV;i++) {            if(i == root) continue;            if(In[i] == MAX)    return -1;//除了跟以外有點沒有入邊,則根無法到達它        }        //2.找環        int cntnode = 0;    memset(ID,-1,sizeof(ID));    memset(vis,-1,sizeof(vis));        In[root] = 0;        for(int i=0;i<NV;i++) {//標記每個環            ret += In[i];            int v = i;            while(vis[v] != i && ID[v] == -1 && v != root) {                vis[v] = i; //vis的作用就是把環中各節點都標記為有入邊的那個點，如，把1,2,3,4,5,6都標記為1                v = pre[v];            }            if(v != root && ID[v] == -1) {                for(int u = pre[v] ; u != v ; u = pre[u]) {                    ID[u] = cntnode;                }                ID[v] = cntnode ++;            }        }        if(cntnode == 0)    break;//無環        for(int i=0;i<NV;i++) if(ID[i] == -1) {            ID[i] = cntnode ++;        }        //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].cost -= In[v];            }        }        NV = cntnode;        root = ID[root];    }    return ret;}int main(){FRE;    int n,m;    while(scanf("%d%d",&n,&m)!=EOF){        int i,j;        int sum=0;        int cnt=0;        for(i=0;i<m;i++){            scanf("%d%d%d",&E[cnt].u,&E[cnt].v,&E[cnt].cost);            E[cnt].u++;            E[cnt].v++;            sum+=E[cnt].cost;            cnt++;        }        sum++;        for(i=1;i<=n;i++){            E[cnt].u=0;            E[cnt].v=i;            E[cnt].cost=sum;            cnt++;        }        int ans=zhuliu(0,n+1,cnt);        if(ans==-1 || ans>=2*sum)puts("impossible");//噹噹前最小樹兩次經過虛擬邊時就不能構成最小樹形圖        else        printf("%d %d\n",ans-sum,pos-m);        printf("\n");    }    return 0;}