【POJ 1273】Drainage Ditches&Dinic详解

想来上一次认真写网络流还是。。。。一年多之前了。。。。
dinic算法写起来简单理解起来简单时间复杂度还棒棒的真是完美啊！
很少能见到时间、空间、编程复杂度都这么低的题目了，业界良心啊！（然而建图业界毒瘤2333）

直接写步骤吧：
1、bfs找出每个点的层数（从源点出发的最短距离）
2、dfs搜索，注意每次必须走到下一层！下一层！同层不能走！
3、找到终点，减去流量，加上反向流量
4、回溯。回溯点u必须满足：
(i)u->v这条边已经流满了；(ii)满足(i)的情况下取最上层
5、dfs结束，回到1
若bfs无法到达汇点，则算法结束。

其中这个4啊听起来很玄乎，但是写的时候似乎是随便乱搞就行？dfs果然方便。

至于时间复杂度，Dinic的上限和EK算法一样是O(V^2*E)，但是Dinic在很多情况下的速度是很快的，一些没有特地卡的题目甚至能跑出线性来。（讲道理啊，Dinic跑不出来还有什么能跑出来的？）

还有就是这个反向边，如果下标从0开始存，那么偶数一定是正向边，奇数一定是反向边，直接异或即可。但如果用指针写惯了之后，从1开始存，就会反一下。。。然后就不能直接异或了。还有我很想知道网上某些反向边直接+1而不是取异或的程序是怎么对的。。

还有的还有就是我还没加上Dinic弧优化。。

#include<cmath>
#include<cstdio>
#include<vector>
#include<queue>
#include<cstring>
#include<iomanip>
#include<stdlib.h>
#include<iostream>
#include<algorithm>
#define ll long long
#define inf 1000000000
#define mod 1000000007
#define N 1000
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
struct NetWorkFlow{int tar,nxt,flow;} edge[N*2];
int head[N],list[N],dep[N];
int tot,n,m,i,x,y,z,res;
void add(int x,int y,int z)
{
    edge[++tot].tar = y; edge[tot].nxt = head[x];
    edge[tot].flow = z; head[x] = tot; 
}
int bfs()
{
    int h = 0; int t = 1; list[1] = 1; dep[1] = 1;
    while (h < t)
        {
            h++; int nw = list[h];
            for (i = head[nw]; i ; i = edge[i].nxt)
                {
                    int nt = edge[i].tar;
                    if (dep[nt] == 0 && edge[i].flow)
                        list[++t] = nt,dep[nt] = dep[nw] + 1;
                }
        }
    if (dep[n] == 0) return 0; else return 1;
}
int dfs(int nw,int MinFlow)
{
    int ResFlow = 0; int i;
    if (nw == n || MinFlow <= 0) return MinFlow; //流量小于等于0退出 
    for (i = head[nw]; i ; i = edge[i].nxt)
        {
            int nt = edge[i].tar;
            if (dep[nt] == dep[nw] + 1) //严格dfs 
                {
                    int temp = dfs(nt,min(MinFlow,edge[i].flow));
                    int f;
                    edge[(i-1)^1+1].flow += temp; //反向边和原边相邻存储 
                    edge[i].flow -= temp;
                    MinFlow -= temp; //目前还能增广的流量 
                    ResFlow += temp;
                    if (MinFlow == 0) break; //回溯的条件 
                }
        }
    return ResFlow;
}
int main()
{
    while (scanf("%d%d",&m,&n) != EOF)
        {
            res = 0; tot = 0; 
            memset(edge,0,sizeof(edge));
            memset(head,0,sizeof(head));
            fo(i,1,m) {scanf("%d%d%d",&x,&y,&z); add(x,y,z); add(y,x,0);}
            memset(dep,0,sizeof(dep));
            while (bfs())
                {
                    res += dfs(1,inf);
                    memset(dep,0,sizeof(dep));
                }
            cout<<res<<endl;
        }
    return 0;
}