网络流(自行理解的Dinic)

自行理解的Dinic 注释即讲解

#include<bits/stdc++.h>
const int maxn=;
using namespace std;
int read()
{
    char ch=getchar();int f=,w=;
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch<=''&&ch>=''){w=w*+ch-'';ch=getchar();}
    return f*w;
}                            //读入优化

struct sj{
    int next;
    int to;
    int w;                    //w是剩余的流量
}a[maxn*];
int head[maxn],size=;        //注意size 要赋值为1 否则会死循环
int m,n;
int s,t;
int cent[maxn];                //层数
int ans;

void add(int x,int y,int z)
{
    a[++size].to=y;
    a[size].next=head[x];
    head[x]=size;
    a[size].w=z;
}                            //加边

bool bfs()
{
    memset(cent,-,sizeof(cent));        //cent 层数先清零
    queue <int> q;
    q.push(s);
    cent[s]=;                            //起点进栈
    while(q.empty()!=)
    {
        int ks=q.front();
        q.pop();
        for(int i=head[ks];i;i=a[i].next)
        {
            int ts=a[i].to;
            if(cent[ts]==-&&a[i].w>)
            {
              q.push(ts);
              cent[ts]=cent[ks]+;        //层数加1
            }
        }
    }
    if(cent[t]<)                        //如果已经找不到增广路
    return ;
    return ;
}

int dfs(int x,int low)                    //dfs搜索 参数:当前的点,当前最小的流量
{
    if(x==t)return low;                    //如果已经搜到终点了 就直接返回
    int flow=,r;                        //避免重复流 每次都新定义一个统计的流量 每层都各自求和
    for(int i=head[x];i;i=a[i].next)
    {
        int st=a[i].to;
        if((a[i].w>)&&(cent[st]==cent[x]+))
        {
            r=dfs(st,min(a[i].w,low-flow));            //递归在这里
            if(r>)
            {
                a[i].w-=r;                            //当前增广路走过的边 残量减掉这么多
                a[i^].w+=r;                        //反向边则加上这么多
                flow+=r;                            //流量加上r 这次找已经流了这么多
            }
        }
    }
    return flow;                                    //递归返回时 则返回已经流完的流量
}

int main()                                //主函数
{
    n=read();m=read();
    s=read();t=read();
    for(int i=;i<=m;i++)
    {
        int x,y,z;
        x=read();
        y=read();
        z=read();
        add(x,y,z);
        add(y,x,);
    }                                    //建立反向边
    while(bfs())
    {
        int row=dfs(s,0x3f3f3f);        //row代表这次搜索找到的流量总和
        ans+=row;
        cout<<row<<endl;
    }                                    //还能找到增广路
    cout<<ans;
    return ;
}