朱刘算法,二分法(比赛网络,UVA 11865)

最小值最大     =====>     二分法


二分最小宽带值m,建图,图中的有向边宽带值均不小于m。然后跑朱刘算法,如果无法建图或者花费过高,那么就放宽条件(减小最小宽带值),如果条件无限宽还不行(m==1),那就不可能了,否则输出结果。


朱刘算法模板:http://blog.csdn.net/leolin_/article/details/6751867


    #define N  1005  
    #define type int  
    #define MAX INT_MAX  
    struct Edge{  
        int u , v;  
        type cost;  
    }E[40005];  
    int pre[N],ID[N],vis[N];  
    type In[N];  
    int n,m;  
    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;  
                    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;  
    }   

代码

    #include<bits/stdc++.h>
    using namespace std;
    typedef long long ll;
    #define N  1005
    #define type ll
    #define MAX LONG_LONG_MAX
    const ll maxm = 2e4;
    struct Edge{
        ll u , v;
        type cost;
    }E[40005];
    ll pre[N],ID[N],vis[N];
    type In[N];

    type zhuliu(ll root,ll NV,ll NE) {//根、结点数、边数
        type ret = 0;
        while(true) {
            //1.找最小入边
            for(ll i=0;i<NV;i++) In[i] = MAX;
            for(ll i=0;i<NE;i++){
                ll u = E[i].u;
                ll v = E[i].v;
                if(E[i].cost < In[v] && u != v) {  //这一步可以把自环切掉,找出每个点的最小入边
                    pre[v] = u;
                    In[v] = E[i].cost;
                }
            }
            for(ll i=0;i<NV;i++) {
                if(i == root) continue;
                if(In[i] == MAX)    return -1;//除了跟以外有点没有入边,则根无法到达它
            }
            //2.找环
            ll cntnode = 0;
        memset(ID,-1,sizeof(ID));
        memset(vis,-1,sizeof(vis));
            In[root] = 0;
            for(ll i=0;i<NV;i++) {//标记每个环
                ret += In[i];
                ll 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(ll u = pre[v] ; u != v ; u = pre[u]) {
                        ID[u] = cntnode;
                    }
                    ID[v] = cntnode ++;
                }
            }
            if(cntnode == 0)    break;//无环
            for(ll i=0;i<NV;i++) if(ID[i] == -1) {
                ID[i] = cntnode ++;
            }
            //3.缩点,重新标记,调整进入环的边的权值,以便下一轮循环找出最小入边
            for(ll i=0;i<NE;i++) {
                ll 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;
    }

    ll n,M,C;
    ll u[maxm],v[maxm],d[maxm],c[maxm];

int main()
{
    ll T;
    scanf("%lld",&T);
    while(T--)
    {
        scanf("%lld %lld %lld",&n,&M,&C);
        ll l=1,r=0;
        for(ll i=0;i<M;i++)
        {
            scanf("%lld %lld %lld %lld",u+i,v+i,d+i,c+i);
            r=max(r,d[i]);
        }
        r++;
        while(l<r)
        {
            ll m=l+(r-l)/2;
            ll xxx=0;
            for(ll i=0;i<M;i++) if(d[i]>=m) E[xxx++]=(Edge){u[i],v[i],c[i]};
            ll yyy=zhuliu(0,n,xxx);
            if(yyy==-1||yyy>C) r=m;
            else l=m+1;
        }
        if(l==1) puts("streaming not possible.");
        else printf("%lld kbps\n",l-1);

    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值