COGS 11. 运输问题1

★★☆   输入文件:maxflowa.in   输出文件:maxflowa.out   简单对比
时间限制:1 s   内存限制:128 MB

【问题描述】
    一个工厂每天生产若干商品,需运输到销售部门进行销售。从产地到销地要经过某些城镇,有不同的路线可以行走,每条两城镇间的公路都有一定的流量限制。请你计算,在不考虑其它车辆使用公路的前提下,如何充分利用所有的公路,使产地运输到销地的商品最多,最多能运输多少商品。
【输入格式】
输入文件有若干行
第一行,一个整数n,表示共有n个城市(2<=n<=100),产地是1号城市,销地是n号城市。
下面有n行,每行有n个数字。第p行第q列的数字表示城镇p与城镇q之间有无公路连接。数字为0表示无,大于0表示有公路,且该数字表示该公路流量。
【输出格式】
输出文件有一行
第一行,1个整数max,表示最大流量为max。
【输入输出样例】
输入文件名: maxflowa.in
6
0 4 8 0 0 0
0 0 4 4 1 0
0 0 0 2 2 0
0 0 0 0 0 7
0 0 0 6 0 9
0 0 0 0 0 0
输出文件名:maxflowa.out
8
 
 
#include <cstdio>
#include <queue>
#define N 60005
#define inf 0x3f3f3f3f
using namespace std;
int n,cnt=1,to[N<<1],dep[N],cur[N],head[N],nextt[N<<1],flow[N<<1];
inline void ins(int u,int v,int w)
{
    nextt[++cnt]=head[u];to[cnt]=v;flow[cnt]=w;head[u]=cnt;
    nextt[++cnt]=head[v];to[cnt]=u;flow[cnt]=0;head[v]=cnt;
}
bool bfs()
{
    for(int i=1;i<=n;++i) dep[i]=-1,cur[i]=head[i];
    queue<int>q;
    dep[1]=0;
    q.push(1);
    for(int u;!q.empty();)
    {
        u=q.front();
        q.pop();
        for(int i=head[u];i;i=nextt[i])
        {
            int v=to[i];
            if(dep[v]==-1&&flow[i])
            {
                dep[v]=dep[u]+1;
                if(v==n) return true;
                q.push(v); 
            }
        }
    }
    return false;
}
inline int min(int a,int b) {return a>b?b:a;}
int dfs(int u,int limit)
{
    if(u==n||!limit) return limit;
    int f,res=0;
    for(int &i=cur[u];i;i=nextt[i])
    {
        int v=to[i];
        if(dep[v]==dep[u]+1&&flow[i]&&(f=dfs(v,min(flow[i],limit))))
        {
            flow[i]-=f;
            flow[i^1]+=f;
            limit-=f;
            res+=f;
            if(!limit) break;
        }
    }
    if(res!=limit) dep[u]=-1;
    return res;
}
int main(int argc,char *argv[])
{
    freopen("maxflowa.in","r",stdin);
    freopen("maxflowa.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
     for(int a,j=1;j<=n;++j)
     {
         scanf("%d",&a);
         if(a) ins(i,j,a);
     }
    int ans=0;
    while(bfs()) ans+=dfs(1,inf);
    printf("%d
",ans);
    return 0;
}
原文地址:https://www.cnblogs.com/ruojisun/p/7684177.html