uva11865 二分+最小树形图（朱刘算法）

/*uva11865 二分+最小树形图（朱刘算法）
注意二分的写法，树形图可作为模板，ps：杭电的模板有问题
*/
#include<iostream>
#include<string.h>
#include<stdio.h>
#include<stdlib.h>
#include<cmath>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#define typec int
#define maxn 100+5
#define maxe 10000+5
using namespace std;

const typec inf=0x3f3f3f3f;
int nextint(){int x;scanf("%d",&x);return x;}
struct Edge{
    int u,v;
    int b;
    typec c;//权值
}edges[maxe],edge[maxe];

int vis[maxn],pre[maxn],id[maxn],in[maxn],cost;
int zhuliu(int root,int nv,int ne){//最小树形图
      int i,j,u,v;
      int ans=0;
      while(1){
              for(i=0;i<nv;i++)in[i]=inf;
              for(i=0;i<ne;i++){
                      u=edge[i].u;
                      v=edge[i].v;
                      if(u!=v&&edge[i].c<in[v]){
                              pre[v]=u;
                              in[v]=edge[i].c;
                      }
              }
              in[root]=0;
              for(i=0;i<nv;i++)if(in[i]==inf)return -1;
              int cnt=0;
              memset(id,-1,sizeof(id));
              memset(vis,-1,sizeof(vis));
              for(i=0;i<nv;i++){
                      v=i;
                      ans+=in[i];
                      while(vis[v]!=i&&id[v]==-1&&v!=root){
                              vis[v]=i;
                              v=pre[v];
                      }
                      if(v!=root&&id[v]==-1){
                              for(u=pre[v];u!=v;u=pre[u])id[u]=cnt;
                              id[v]=cnt++;
                      }
              }
              if(cnt==0)break;
              for(i=0;i<nv;i++)if(id[i]==-1)id[i]=cnt++;
              for(i=0;i<ne;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];
              }
              nv=cnt;root=id[root];
      }
        return ans;
}
int T,Cost,n,m;
int cnt;
void built_tree(int M){
    cnt=0;
    for(int i=0;i<m;i++)
    if (edges[i].b>=M) edge[cnt++]=edges[i];
}
int main(){
//    freopen("out.txt","w",stdout);
    cin>>T;
    while(cin>>n>>m>>Cost){
        int L=0,R=0;
        for(int i=0;i<m;i++) {
            int u,v,b,c;
            u=nextint();v=nextint();b=nextint();c=nextint();
            edges[i]=(Edge){u,v,b,c};
            R=max(R,edges[i].b);
        }
        bool ok=false;int ans=0;
        while(L<R){
            int M=(L+R+1)/2;
            built_tree(M);
            int mc=zhuliu(0,n,cnt);
            if (mc>0 && mc<=Cost) {
                ok=true;
                ans=L=M;
            }else R=M-1;
        }
        if (ok) printf("%d kbps
",ans);else printf("streaming not possible.
");

    }
    return 0;
}