HDU 6141 I am your Father!（最小树形图）

Description

给出一个$n$个点$m$条边的有向图，求以$1$为根的最大费用有向树，在权值和相同的情况下选择$n$节点的父亲节点编号最小的方案，要求输出最小费用和$n$节点的父亲节点编号

Input

第一行一整数$T$表示用例组数，每组用例首先输入两整数$n$和$m$分别表示点数和边数，之后$m$行每行输入三个整数$x,y,w$表示$x$到$y$有一条权值为$w$的边$(1\le n\le 10^3,m\le 10^4,1\le w\le 100)$

Output

对于每组用例，输出最大费用有向树的权值和和$n$节点的父亲节点编号

Sample Input

2
3 3
1 2 10
1 3 10
2 3 10
3 3
1 2 10
1 3 10
2 3 11

Sample Output

20 1
21 2

Solution

模版题，通过改变边权使得$n$节点父亲的字典序最小，把点编号$0$~$n-1$，把每条边的边权乘上$-n$，把$u$到$n$的边边权加上$u$，设新图的最小费用为$-ans$，则最大费用为$\lceil \frac{ans}{n}\rceil$，$n$节点的父亲节点编号为$\lceil \frac{ans}{n}\rceil\cdot n-ans+1$

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 1111
struct Edge
{
    int u,v,c;
}edge[10*maxn];
int pre[maxn],id[maxn],vis[maxn],in[maxn];
int zhuliu(int root,int n,int m)
{
    int ans=0,u,v;
    while(1)
    {
        for(int i=0;i<n;i++)in[i]=INF;
        for(int i=0;i<m;i++)
            if(edge[i].u!=edge[i].v&&edge[i].c<in[edge[i].v])
            {
                pre[edge[i].v]=edge[i].u;
                in[edge[i].v]=edge[i].c;
            }
        for(int i=0;i<n;i++)
            if(i!=root&&in[i]==INF) 
                return -1;
        int res=0;
        memset(vis,-1,sizeof(vis));
        memset(id,-1,sizeof(id));
        in[root]=0;
        for(int i=0;i<n;i++)
        {
            ans+=in[i];
            v=i;
            while(vis[v]!=i&&id[v]==-1&&v!=root)
            {
                vis[v]=i;
                v=pre[v];
            }
            if(id[v]==-1&&v!=root)
            {
                for(u=pre[v];u!=v;u=pre[u])
                    id[u]=res;
                id[v]=res++;
            }
        }   
        if(res==0)break;
        for(int i=0;i<n;i++)
            if(id[i]==-1)
                id[i]=res++;
        for(int i=0;i<m;i++)
        {
            v=edge[i].v;
            edge[i].u=id[edge[i].u];
            edge[i].v=id[edge[i].v];
            if(edge[i].u!=edge[i].v)    
                edge[i].c-=in[v];
        }
        n=res;
        root=id[root];
    }
    return ans;
}
int T,n,m;
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%d",&edge[i].u,&edge[i].v,&edge[i].c);
            edge[i].u--,edge[i].v--,edge[i].c*=-n;
            if(edge[i].v==n-1)edge[i].c+=edge[i].u;

        }
        int ans=-zhuliu(0,n,m);
        printf("%d %d\n",(ans+n-1)/n,(ans+n-1)/n*n-ans+1);
    }
    return 0;
}