hdu 2121 最小树形图 +虚根

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

 

题意：

        n个城市，m条单向边，每条边都有一个权值，现在要你求出这样一个节点，由这个节点出发能走遍所有的城市，同时使总权值最小。如果找不到就输出impossible，能的话输出最小权值和这个节点的编号。

 

做法：

       最小树形图虚根的第一道题，因为不能确定哪个点开始是最好的，所以要建立一个虚根，这个根和所有的节点都要连一条无穷大（其实只要大于所有路的总权值就好啦）的边，算出答案之后把这个无穷大减掉就好了，这样的话每个点都有机会可以作为根，如果这个答案>=2*sum的话，很明显是由虚根多连了一条边，那么就是不合法的。

        这里找真正的根这里有个小技巧，我开始的时候是想找节点的pre，找到第一个pre为n+1这个点的就是真正的根，但是wa了，为什么我也。。不知道。然后去网上找了别人的做法，都是在算法中找点的时候就改掉了用pos=i然后再最后的时候减掉m。

        我的理解大概就是，因为我们的边是在加完图里面之后再加的虚根的边，所以如果图内部已经有连向点x的边，那么for到虚根这一块的边的时候就不会进入这个if语句，也就不会被修改了（因为虚根连出来的点都是要大于总权值的），如果图内有环，那么在缩点之后，for的点从小到大开始，那么也是一定会先到较小的点的，也是符合条件的。至于最后要减掉m就比较好理解了，因为在加虚根之前就有m条边，从点1开始for，pos是哪个就是哪个喽。

 

代码如下：

     

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
int vis[1005],idx[1005],pre[3005],cnt,pos;
ll in[1005];
struct node{
    int u,v;
    ll w;
}edge[1000005];
ll solve(int t,int n,int m){
    int x,y,i,j;
    ll ans=0;
    while(1){
        for(i=1;i<=n;i++)in[i]=inf;
        for(i=1;i<=m;i++){
            node e=edge[i];
            if(e.u!=e.v&&e.w<in[e.v]){
                in[e.v]=e.w;
                pre[e.v]=e.u;
                if(e.u==t) pos=i;
            }
        }
        for(i=1;i<=n;i++){
            if(in[i]==inf&&i!=t)return -1;
            if(i!=t)ans+=in[i];
        }
        int tot=0;
        memset(vis,-1,sizeof(vis));
        memset(idx,-1,sizeof(idx));
        for(i=1;i<=n;i++){
            y=i;
            while(vis[y]!=i&&idx[y]==-1&&y!=t){
                vis[y]=i;
                y=pre[y];
            }
            if(y!=t&&idx[y]==-1){
                tot++;
                idx[y]=tot;
                for(x=pre[y];x!=y;x=pre[x])idx[x]=tot;
            }
        }
        if(!tot)break;
        for(i=1;i<=n;i++)if(idx[i]==-1)idx[i]=++tot;
        for(i=1;i<=m;i++){
            int y=edge[i].v;
            edge[i].u=idx[edge[i].u];
            edge[i].v=idx[edge[i].v];
            if(edge[i].u!=edge[i].v)edge[i].w-=in[y];
        }
        n=tot;
        t=idx[t];
    }
    return ans;
}
int main(){
    int n,m,cas=0;
    while(scanf("%d%d",&n,&m)!=EOF){
        int i,j;
        cnt=0;
        ll sum=0;
        for(i=1;i<=m;i++){
            int u,v;
            ll c;
            scanf("%d%d%lld",&u,&v,&c);
            u++,v++;
            edge[i].u=u;
            edge[i].v=v;
            edge[i].w=c;
            sum+=c;
        }
        sum+=1;
        int now=0;
        for(int i=1;i<=n;i++){
            now++;
            edge[m+now].u=n+1;
            edge[m+now].v=i;
            edge[m+now].w=sum;
        }
        ll ans=solve(n+1,n+1,m+now);
        if(ans==-1||ans-sum>=sum)printf("impossible\n\n");
        else {
            printf("%lld %d\n\n",ans-sum,pos-m-1);
        }
    }
    return 0;
}