ZOJ 3496 Assignment | 二分+有上下界网络流

Posted MSPqwq

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了ZOJ 3496 Assignment | 二分+有上下界网络流相关的知识,希望对你有一定的参考价值。

题目:

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3496

大概意思:给你一个网络,有源汇,在保证最大流的情况下求下面两个问题答案

1.所有边中流量最大的边流量最小

2.所有边中流量最小的边流量最大


题解:

De了一下午啊啊,之前学的上下界网络流有问题!

对于问题一,我们二分最大流量,每次建图跑最大流,看是不是和之前一样即可

对于问题二,我们同样二分答案lim,这样每条边满足流量限制w[i]∈[lim,c[i]]

这样建图跑有源汇最大流即可.

下面是坑点

之前写的有源汇最大流都是先跑可行流然后把超级源和超级汇删掉再跑,答案还得加加减减,但是就是Wa

后来发现其实并没有那么难,在可行流的残余网络直接再跑一边最大流就是答案

感性证明如下:

因为如果可行的话超级源的所有出边都满流,这样算最大流的时候只会算到t->s的反边上,又因为原来的可行流就在这条边上,这样直接就能算出答案

别忘了输出答案*p

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
typedef long long ll;
#define N 505
#define M 400005
#define INF 0x3f3f3f3f
using namespace std;
ll Case,n,m,t,s,P,ans1,ans2,S,T,ecnt,Maxflow,l,r,u[M],v[M],c[M],mid,Maxc,lim;
ll head[N],lev[N],cur[N],du[N];
queue <ll> q;
struct adj
{
    ll nxt,v,w;
} e[M];
ll read()
{
    ll ret=0,neg=1;
    char j=getchar();
    for (; j>9 || j<0; j=getchar())
        if (j==-) neg=-1;
    for (; j>=0 && j<=9; j=getchar())
        ret=ret*10+j-0;
    return ret*neg;
}
void add(ll u,ll v,ll w)
{
    e[++ecnt].v=v;e[ecnt].w=w;e[ecnt].nxt=head[u];head[u]=ecnt;
    e[++ecnt].v=u;e[ecnt].w=0;e[ecnt].nxt=head[v];head[v]=ecnt;
}
void init()
{
    ecnt=1;
    memset(head,0,sizeof(head));
}
bool Bfs()
{
    while (!q.empty()) q.pop();
    for (ll i=1; i<=lim; i++)
        cur[i]=head[i],lev[i]=-1;
    q.push(S),lev[S]=1;
    while (!q.empty())
    {
        ll u=q.front();
        q.pop();
        for (ll i=head[u],v; i; i=e[i].nxt)
            if (lev[v=e[i].v]==-1 && e[i].w>0)
            {
                q.push(v),lev[v]=lev[u]+1;
                if (v==T) return 1;
            }
    }
    return 0;
}
ll Dfs(ll u,ll flow)
{
    if (u==T) return flow;
    ll ret=0,delta;
    for (ll &i=cur[u],v; i; i=e[i].nxt)
        if (e[i].w>0 && lev[v=e[i].v]==lev[u]+1)
        {
            delta=Dfs(v,min(e[i].w,flow-ret));
            if (delta)
            {
                e[i].w-=delta;
                e[i^1].w+=delta;
                ret+=delta;
                if (ret==flow) break;
            }
        }
    return ret;
}
ll getG(ll lim)
{
    ll sum=0,tmp=0;
    init();
    add(t,s,INF);S=n+1,T=n+2;
    memset(du,0,sizeof(du));
    for (ll i=1; i<=m; i++)
        if (c[i]-lim<0) return -1;
        else add(u[i],v[i],c[i]-lim),du[u[i]]-=lim,du[v[i]]+=lim;
    for (ll i=1; i<=n; i++)
    {
        if (du[i]>0) add(S,i,du[i]),sum+=du[i];
        if (du[i]<0) add(i,T,-du[i]);
    }
    while (Bfs()) tmp+=Dfs(S,INF);
    if (tmp!=sum) return -1;
    tmp=0;S=s;T=t;
    while (Bfs()) tmp+=Dfs(S,INF);
    return tmp;
}
int main()
{
    Case=read();
    while (Case--)
    {
        n=read(),m=read(),s=read(),t=read(),P=read();
        init();
        S=++s;T=++t;
        Maxflow=Maxc=l=0;
        lim=n+2;
        for (ll i=1; i<=m; i++)
            u[i]=read(),v[i]=read(),c[i]=read(),add(++u[i],++v[i],c[i]),r=Maxc=max(Maxc,c[i]);
        while (Bfs()) Maxflow+=Dfs(S,INF);
        while (l<r)
        {
            init();
            ll mid=l+r>>1,tmp=0;
            for (ll i=1; i<=m; i++) add(u[i],v[i],min(c[i],mid));
            while (Bfs()) tmp+=Dfs(S,INF);
            if (tmp==Maxflow) r=mid;
            else l=mid+1;
        }
        ans1=l;l=0;r=Maxc;
        while (l<r)
        {
            ll mid=l+r+1>>1,tmp=0;
            if (getG(mid)==Maxflow) l=mid;
            else r=mid-1;
        }
        ans2=l;
        printf("%lld %lld\n",1ll*ans1*P,1ll*ans2*P);
    }
    return 0;
}

 

以上是关于ZOJ 3496 Assignment | 二分+有上下界网络流的主要内容,如果未能解决你的问题,请参考以下文章

[zoj] 3496 Assignment || 有源汇上下界最大流

ZOJ3496 Assignment

POJ3189_Steady Cow Assignment(二分图多重匹配/网络流+二分构图)

POJ3189 Steady Cow Assignment(二分图多重匹配)

HDU - 5289 Assignment (RMQ+二分)

POJ3189 Steady Cow Assignment —— 二分图多重匹配/最大流 + 二分