网络流最大流 Dinic算法模板

Dinic算法跟EK算法的优化地方就在于，EK是每次bfs记录前驱结点，然后从汇点找到一条回到源点的路。

而dinic是bfs记录深度，不断从源点dfs到汇点并且更新路径流量，直到dfs不到汇点。

邻接矩阵代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
#define pr pair<int,int>
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxn=405;
int dis[maxn],c[maxn][maxn];
int n,m;
int bfs(int st,int ed)
{
    memset(dis,-1,sizeof(dis));//深度
    dis[st]=0;
    queue<int>q;
    q.push(1);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int v=1;v<=n;v++)
        {
            if(dis[v]<0&&c[u][v]>0)//没有遍历过的点
            {
                dis[v]=dis[u]+1;
                q.push(v);
            }
        }
    }
    if(dis[ed]>0) return true;//如果都没有跑的汇点，意思没有增广路
    return false;
}
int find(int u,int low)//其实是一个dfs的过程，上边bfs记录深度也是因为这里通过dfs跑出多次增广路
{
    int a=0;
    if(u==n) return low;
    for(int i=1;i<=n;i++)
    {
        if(c[u][i]>0&&dis[i]==dis[u]+1&&(a=find(i,min(low,c[u][i]))))//3个条件
            //1.这条路有残留量2.这两点相连，通俗的讲就是这次找增广路的过程中，水就是沿着这条路的，但是下次就不一定.3就是递归下一个点看能不能到汇点
            //同时在过程中找到这条路最小流量，这种方法在稠密的图中有很大的优化作用
        {
            c[u][i]-=a;//更新残留量
            c[i][u]+=a;
            return a;
        }
    }
    return 0;
}
int main()
{
    int x,y,z;
    while(cin>>m>>n)
    {
        memset(c,0,sizeof(c));
        for(int i=0;i<m;i++)
            scanf("%d %d %d",&x,&y,&z),c[x][y]+=z;
        int ans=0,temp;
        while(bfs(1,n))
        {
            while(temp=find(1,inff))
                ans+=temp;
        }
        printf("%d\n",ans);
    }
    return 0;
}

邻接表板子：

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
#define pr pair<int,int>
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxn=205;
int d[maxn];
int head[maxn],sign;
struct node
{
    int to,p,val;
}edge[maxn<<1];
int n,m;
void add(int u,int v,int val)
{
    edge[sign]=node{v,head[u],val};
    head[u]=sign++;
}
void init()
{
    sign=0;
    for(int i=0;i<=n;i++)
        head[i]=-1;
}
int find(int u,int low)
{
    int a=0;
    if(u==n) return low;
    for(int i=head[u];~i;i=edge[i].p)
    {
        int v=edge[i].to;
        if(d[v]==d[u]+1&&edge[i].val>0&&(a=find(v,min(low,edge[i].val))))
        {
            edge[i].val-=a;
            edge[i^1].val+=a;
            return a;
        }
    }
    return 0;
}
int bfs(int st,int ed)
{
    memset(d,0,sizeof(d));
    d[st]=1;
    queue<int> q;
    q.push(st);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];~i;i=edge[i].p)
        {
            int v=edge[i].to;
            if(!d[v]&&edge[i].val>0)
            {
                d[v]=d[u]+1;
                q.push(v);
            }
        }
    }
    if(d[ed]>0) return true;
    return false;
}
int main()
{
    int x,y,z;
    while(cin>>m>>n)
    {
        init();
        for(int i=0;i<m;i++)
        {
            scanf("%d %d %d",&x,&y,&z);
            add(x,y,z),add(y,x,0);
        }
        int ans=0,temp;
        while(bfs(1,n))
        {
            while(temp=find(1,inff))
                ans+=temp;
        }
        printf("%d\n",ans);
    }
}