Uva-1440-Inspection (只有下界的网络最小流)

来源：互联网发布：重庆宝钻资产管理知乎编辑：程序博客网时间：2024/06/07 02:27

Inspection

题目：

题目还是比较裸的没有上界只有下界的网络最小流问题。输出路劲太烦人了。！不想写了。给出别人的代码。

算法分析：

网络最小流一般用二分查找判断就好了。要注意的就是在最后输出路径的时候要在调用一次最大流。因为，之前结束的流量不一定是最后答案的流量。

#include<cstdio>#include<cstring>#include<iostream>#include<queue>using namespace std;const int inf=1<<29;const int maxn=130;const int maxm=maxn*maxn*3;int n,e,st,s_st,s_des,des,head[maxn],nxt[maxm],pnt[maxm],flow[maxm],level[maxn];int c_in[maxn],c_out[maxn],out[maxn];bool go[maxn][maxn];queue<int> q;void AddEdge(int u,int v,int f){    //printf("Add %d %d %d\n",u,v,f);    pnt[e]=v;nxt[e]=head[u];flow[e]=f;head[u]=e++;    pnt[e]=u;nxt[e]=head[v];flow[e]=0;head[v]=e++;}bool BFS(){    memset(level,0,sizeof(level));    level[st]=1;    q.push(st);    while(!q.empty())    {        int u=q.front();        q.pop();        for(int i=head[u];i!=-1;i=nxt[i])            if(flow[i]&&!level[pnt[i]])            {                level[pnt[i]]=level[u]+1;                q.push(pnt[i]);            }    }    return level[des];}int DFS(int u,int maxf){    if(u==des||!maxf)        return maxf;    for(int i=head[u],t;i!=-1;i=nxt[i])        if(level[pnt[i]]==level[u]+1&&(t=DFS(pnt[i],min(maxf,flow[i]))))        {            flow[i]-=t;            flow[i^1]+=t;            return t;        }    return level[u]=0;}int maxflow(int s,int t){    st=s;des=t;    int ans=0;    while(BFS())        while(1)        {            int val=DFS(st,inf);            if(!val)                break;            ans+=val;        }    return ans;}bool check(int val){    e=st=0,des=n+1,s_st=n+2,s_des=n+3;    memset(head,-1,sizeof(head));    for(int i=1;i<=n;i++)        for(int j=1;j<=n;j++)            if(go[i][j])                AddEdge(i,j,inf);    for(int i=1;i<=n;i++)    {        AddEdge(st,i,inf);        AddEdge(i,des,inf);    }    AddEdge(des,st,val);    int sum=0;    for(int i=1;i<=n;i++)    {        if(c_in[i]>c_out[i])        {            AddEdge(s_st,i,c_in[i]-c_out[i]);            sum+=c_in[i]-c_out[i];        }        if(c_in[i]<c_out[i])            AddEdge(i,s_des,c_out[i]-c_in[i]);    }    int f=maxflow(s_st,s_des);    return f==sum;}int cap[maxn][maxn];int dfs(int u,int cnt){    out[cnt]=u;    if(u==des)    {        for(int i=1;i<cnt;i++)            printf("%d%c",out[i],i==cnt-1?'\n':' ');        return 1;    }    for(int i=1;i<=des;i++)        if(cap[u][i])        {            cap[u][i]--;            if(dfs(i,cnt+1))                return 1;        }    return 0;}void solve(){    int sl=0,sr=inf,ans=0;    while(sl<=sr)    {        int mid=(sl+sr)>>1;        if(check(mid))        {            sr=mid-1;                    ans=mid;        }        else            sl=mid+1;    }    printf("%d\n",ans);    check(ans);    st=0;des=n+1;    for(int i=1;i<=n;i++)        cap[st][i]=cap[i][des]=1;    for(int u=0;u<=n;u++)        for(int i=head[u];i!=-1;i=nxt[i])            if(cap[u][pnt[i]])                cap[u][pnt[i]]=flow[i^1]+1-(u==st||pnt[i]==des);    for(int i=1;i<=des;i++)        while(cap[st][i])        {            cap[st][i]--;            dfs(i,1);        }}int main(){    while(scanf("%d",&n)!=EOF)    {        memset(c_in,0,sizeof(c_in));        memset(c_out,0,sizeof(c_out));        memset(go,0,sizeof(go));        for(int i=1;i<=n;i++)        {            int m;            scanf("%d",&m);            for(int j=0;j<m;j++)            {                int v;                scanf("%d",&v);                go[i][v]=1;                cap[i][v]=1;                c_out[i]++;                c_in[v]++;            }        }        solve();    }    return 0;}

1 0