网络流

最新推荐文章于 2025-09-12 17:31:31 发布

weixin_34247299

最新推荐文章于 2025-09-12 17:31:31 发布

阅读量57

点赞数

CC 4.0 BY-SA版权

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/Lance1ot/p/8962864.html

本文介绍了一种高效的网络流算法——Dinic算法，并提供了详细的代码实现。该算法首先使用BFS进行分层，然后在此基础上运行DFS来寻找增广路径。

dinic
网络流

先利用bfs进行分层，然后再在上面跑dfs.

用的是dinic

具体还是看代码把

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
struct node
{
    int point;
    int value;
    int nxt;
};
node line[200100];
int head[50100],tail=-1;//利用的yi
int cur[50100],maxflow;
int dep[50100];
int inf=0x7fffffff;
int n,m,s,t;
queue<int>q;
void add(int x,int y,int dis,int base)
{
    line[++tail].point=y;
    if(base)    
        line[tail].value=dis;
    line[tail].nxt=head[x];
    head[x]=tail;
}
bool bfs(int start,int end)
{
    memset(dep,0,sizeof(dep));
    while(!q.empty())
        q.pop();
    for(int i=1;i<=n;i++)
        cur[i]=head[i];
    dep[start]=1;
    q.push(start);
    while(!q.empty())
    {
        int now=q.front();
        q.pop();
        for(int i=head[now];i!=-1;i=line[i].nxt)
        {
            if(!dep[line[i].point]&&line[i].value)
            {
                dep[line[i].point]=dep[now]+1;
                q.push(line[i].point);
            }
        }
    }
    if(dep[end])
        return true;
    else
        return false;
}
int dfs(int now,int aim,int limte)
{
    if(now==aim||!limte)
        return limte;
    int flow=0,f;
    for(int i=cur[now];i!=-1;i=line[i].nxt)
    {
        cur[now]=i;
        if(dep[line[i].point]==dep[now]+1&&(f=dfs(line[i].point,aim,min(limte,line[i].value))))
        {
            flow+=f;
            limte-=f;
            line[i].value-=f;
            line[i^1].value+=f;
            if(!limte)
                break;
        }
    }
    return flow;
    
}
int dinic(int s,int e)
{
    int ans=0;
    while(bfs(s,e))//利用bfs分层
        ans+=dfs(s,e,inf);//dfs寻找增广路
    return ans;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t);//读入数据
        int a,b,c;
    for(int i=1;i<=n;i++)
        head[i]=-1;//初始化链表
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c,1);//最后一个参数是判断是否是正向边还是反向边
        add(b,a,c,0);
    }
    printf("%d",dinic(s,t));
}

转载于:https://www.cnblogs.com/Lance1ot/p/8962864.html