[COGS 0011] 运输问题1

11. 运输问题1

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

【问题描述】

    一个工厂每天生产若干商品，需运输到销售部门进行销售。从产地到销地要经过某些城镇，有不同的路线可以行走，每条两城镇间的公路都有一定的流量限制。请你计算，在不考虑其它车辆使用公路的前提下，如何充分利用所有的公路，使产地运输到销地的商品最多，最多能运输多少商品。

【输入格式】

输入文件有若干行
第一行，一个整数n，表示共有n个城市$(2 leq n leq 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

好了对于这么水的题我表示只是来存代码的(逃

注意读入建图的时候的反向边标记然后大力Dinic就好了w

GitHub

#include <queue>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>

const int MAXV=110;
const int MAXE=10010;
const int INF=0x7FFFFFFF;

struct Edge{
    int from;
    int to;
    int flow;
    Edge* rev;
    Edge* next;
};
Edge E[MAXE];
Edge* head[MAXV];
Edge* top=E;

int v;
int depth[MAXV];
int flow[MAXV][MAXV];

int Dinic(int,int);
void Initialize();
void Insert(int,int);
int DFS(int,int,int);
bool BFS(int,int);

int main(){
    freopen("maxflowa.in","r",stdin);
    freopen("maxflowa.out","w",stdout);
    Initialize();
    printf("%d
",Dinic(1,v));
    return 0;
}

int Dinic(int s,int t){
    int ans=0;
    while(BFS(s,t)){
        ans+=DFS(s,INF,t);
    }
    return ans;
}

int DFS(int s,int flow,int t){
    if(s==t||flow==0)
        return flow;
    int tmp=flow;
    int k;
    for(Edge* i=head[s];i!=NULL;i=i->next){
        if(i->flow!=0&&tmp!=0&&depth[i->to]==depth[s]+1){
            k=DFS(i->to,std::min(tmp,i->flow),t);
            if(k==0){
                depth[i->to]=0;
                continue;
            }
            tmp-=k;
            i->flow-=k;
            i->rev->flow+=k;
            if(tmp==0)
                break;
        }
    }
    return flow-tmp;
}

bool BFS(int s,int t){
    memset(depth,0,sizeof(depth));
    std::queue<int> q;
    depth[s]=1;
    q.push(s);
    while(!q.empty()){
        s=q.front();
        q.pop();
        for(Edge* i=head[s];i!=NULL;i=i->next){
            if(depth[i->to]==0&&i->flow!=0){
                q.push(i->to);
                depth[i->to]=depth[s]+1;
                if(i->to==t)
                    return true;
            }
        }
    }
    return false;
}

void Initialize(){
    scanf("%d",&v);
    for(int i=1;i<=v;i++){
        for(int j=1;j<=v;j++){
            scanf("%d",&flow[i][j]);
        }
    }
    for(int i=1;i<=v;i++){
        for(int j=1;j<i;j++){
            Insert(i,j);
        }
    }
}

inline void Insert(int a,int b){
    top->from=a;
    top->to=b;
    top->flow=flow[a][b];
    top->next=head[a];
    head[a]=top;
    top->rev=top+1;
    top++;
    top->from=b;
    top->to=a;
    top->flow=flow[b][a];
    top->next=head[b];
    head[b]=top;
    top->rev=top-1;
    top++;
}