UVA 11865 Stream My Contest（最小树形图、朱刘算法）

题目：http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2965

题目大意：有n个点，m条有向边，和一个定值cost。n个点，编号0~n-1，0为服务器，只能发数据，其他为主机。每条边有一个带宽和价格，问你在不超过cost 的情况下，整个网络的最小带宽的最大值。

解题思路：这种题目，先一个，就是二分，确定这个最小带宽（筛选边），然后再判断最小的总价值是不是在 cost 以内即可。怎么判断呢？那就是确定 root 的最小树形图了，用朱刘算法（看到了中国人发明的算法，好开心 = =）。关于这个算法，看下面这个博客的三个链接吧：我是链接。时间复杂度为 O（VE）。

        不得不说这个算法的换边和迭代过程很精妙啊！关于这个实现过程，也很有讲究的，找环给编号那里需要好好地琢磨琢磨，其实就和 dfs 一样，当然用 dfs 实现也没有问题，但是这样写比较简单，也不用再写一个 dfs 函数了。起先我一直有一个地方搞不清楚，就是代码中换边那里，如果 id 不一样就要减，那么如果本身连接这两个点的边就不是在环里呢，不是不应该减吗，只有环里的才需要换边呀？其实这里还有地方，那就是 ret ，它只在函数刚进来的时候 = 0。为什么 while 里面每次算的时候不需要初始化为 0 呢？这两个问题连起来就能想明白了。如果这两个点的不是在环里，那么把减掉，但是前面并没有缩点，这样这条边的值就为 0 了，以后加 in[ v ] 也会是 + 0。所以这样写是没有问题的。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int INF = 0x0fffffff;
const int MAXN = 66;
const int MAXM = 11111;

struct Edge
{
    int s,t,val,b;
    void read()
    {
        scanf("%d%d%d%d",&s,&t,&b,&val);
    }
} edge[MAXM<<1];

int in[MAXN],pre[MAXN],id[MAXN],vis[MAXN];

int directed_mst(int root,int n,int m,int mid)
{
    int ret = 0;
    while(1)
    {   
        //1.找最小入边
        for(int i = 0;i < n;i++)
            in[i] = INF;
        for(int i = 0;i < m;i++)
        {
            if(edge[i].b < mid) continue;
            int u = edge[i].s;
            int v = edge[i].t;
            if(edge[i].val < in[v] && u != v)
            {
                in[v] = edge[i].val;
                pre[v] = u;
            }
        }
        for(int i = 0;i < n;i++)
        {
            if(i == root) continue;
            else if(in[i] == INF) return -1;//除了根以外，有点没有入边，则根无法到达它
        }
        //2.找环
        int cnt_node = 0;
        memset(id,-1,sizeof(id));
        memset(vis,-1,sizeof(vis));
        in[root] = 0;
        for(int i = 0;i < n;i++)
        {
            ret += in[i];
            int v = i;
            while(vis[v] != i && id[v] == -1 && v != root)
            {
                vis[v] = i;
                v = pre[v];
            }
            if(id[v] == -1 && v != root)
            {
                for(int u = pre[v]; u != v;u = pre[u])
                    id[u] = cnt_node;
                id[v] = cnt_node++;
            }
        }
        if(cnt_node == 0) break;//无环，即为答案
        for(int i = 0;i < n;i++)
            if(id[i] == -1) id[i] = cnt_node++;
        //3.缩点，重新建图
        for(int i = 0;i < m;i++)
        {
            if(edge[i].b < mid) continue;
            int v = edge[i].t;
            edge[i].s = id[edge[i].s];
            edge[i].t = id[edge[i].t];
            if(edge[i].s != edge[i].t)
                edge[i].val -= in[v];//换边
        }
        n = cnt_node;
        root = id[root];
    }
    return ret;
}

int main()
{
    int _;
    scanf("%d",&_);
    while(_--)
    {
        int n,m,cost;
        scanf("%d%d%d",&n,&m,&cost);
        for(int i = 0; i < m;i++)
        {
            edge[i].read();
            edge[i+m] = edge[i];
        }
        int ans = 0;
        int l = 0,r = 1e6;
        while(l <= r)
        {
            for(int i = 0;i < m;i++)
                edge[i] = edge[i+m];

            int mid = (l+r)>>1;
            int tmp = directed_mst(0,n,m,mid);
            //printf("mid = %d,tmp = %d\n",mid,tmp);
            if(tmp != -1 && tmp <= cost)
            {
                ans = mid;
                l = mid+1;
            }
            else r = mid-1;
        }
        if(ans == 0)
            puts("streaming not possible.");
        else printf("%d kbps\n",ans);
    }
    return 0;
}