深度优先搜索求最大流

深度优先搜索求最大流

---

题目

Luogu P2740 [USACO4.2]草地排水Drainage Ditches
题目背景

在农夫约翰的农场上，每逢下雨，贝茜最喜欢的三叶草地就积聚了一潭水。这意味着草地被水淹没了，并且小草要继续生长还要花相当长一段时间。因此，农夫约翰修建了一套排水系统来使贝茜的草地免除被大水淹没的烦恼（不用担心，雨水会流向附近的一条小溪）。作为一名一流的技师，农夫约翰已经在每条排水沟的一端安上了控制器，这样他可以控制流入排水沟的水流量。

题目描述

农夫约翰知道每一条排水沟每分钟可以流过的水量，和排水系统的准确布局（起点为水潭而终点为小溪的一张网）。需要注意的是，有些时候从一处到另一处不只有一条排水沟。

根据这些信息，计算从水潭排水到小溪的最大流量。对于给出的每条排水沟，雨水只能沿着一个方向流动，注意可能会出现雨水环形流动的情形。

输入输出格式

输入格式：
第1行: 两个用空格分开的整数N (0 <= N <= 200) 和 M (2 <= M <= 200)。N是农夫John已经挖好的排水沟的数量，M是排水沟交叉点的数量。交点1是水潭，交点M是小溪。

第二行到第N+1行: 每行有三个整数，Si, Ei, 和 Ci。Si 和 Ei (1 <= Si, Ei <= M) 指明排水沟两端的交点，雨水从Si 流向Ei。Ci (0 <= Ci <= 10,000,000)是这条排水沟的最大容量。

输出格式：
输出一个整数，即排水的最大流量。

输入输出样例

输入样例#1：

5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10

输出样例#1：

50

大致思路

题意简化：给定一个有n个点m条边的图，图中的边有以下属性：

两个端点u(起点)和v(终点)，和边的容量w。有方向。

现在要从原点1把一些物资运输到终点n，每条边在通过的物资数量大于w时就会粉碎（被删除失效），求最多能把多少物资运送到终点？

思路
对于每个边，我们把他的容量记为$C(u,v)$，实际运送的物品量叫流量，记为$f(u,v)$，剩下的可用容量叫残量$c(u,v)=C(u,v)-f(u,v)$
那么我们可以推知以下结论：

1. $f(u,v)$和$f(v,u)$中，最多有1个正数。且$f(u,v)+f(v,u)=0$（斜对称性）
2. $f(u,v)<=C(u,v)$（容量限制）
3. 对于除了起点和终点的结点u，$\sum_{(u,v)\in Edge} f(u,v)=0$(流量平衡)

算法思想，就是从$f(all)=0$开始，在满足上述三个性质的情况下增加流量。残量网络中的任何一条从起点到终点的路径都对应一条原图中的增广路。只要求出改路径中的残量最小值$d$，把对应边的流量增加$d$即可。当残量网络里没有增广路，则当前流就是最大流。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=1000+10;
inline void read(int &misaka) {
    int ans=0,flag=1;
    char ch=getchar();
    while(ch<'0' || ch>'9') {
        if(ch=='-')
            flag=-1;
        ch=getchar();
    }
    while(ch>='0' && ch<='9') {
        ans=ans*10+(int)(ch-'0');
        ch=getchar();
    }
    misaka=ans*flag;
}
struct E {
    int    next,to;
    int c;
} edge[maxn];
int head[maxn],cnt=0;
int n,m;
int ans=0,pf;
bool vis[maxn];
inline void add(int u,int v,int w) {
    edge[cnt].to=v;
    edge[cnt].c=w;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}
inline bool dfs(int x,int f) {
    if(x==n) {
        ans+=f;
        pf=f;
        return true;
    }
    vis[x]=1;
    for(int i=head[x]; ~i; i=edge[i].next) {
        if(!vis[edge[i].to]&&edge[i].c>0) {
            if(dfs(edge[i].to,min(f,edge[i].c))) {
                edge[i].c-=pf;
                edge[i^1].c+=pf;
                return true;
            }
        }
    }
    return false;
}
int main() {
    memset(head,-1,sizeof(head));
    ios::sync_with_stdio(false);
    read(m);read(n);
    int x,y,z;
    for(int i=1; i<=m; i++) {
        read(x);read(y);read(z);
        add(x,y,z);
        add(y,x,0);
    }
    while(dfs(1,100000000)) {
        pf=0;
        memset(vis,0,sizeof(vis));
    }
    cout<<ans<<endl;
    return 0;
}