• hdu1532 网络流最大流问题


      这是我的网络流第一题,用的是增广路算法即EK算法,算法思想:每次用BFS找一条最短的增广路径,然后沿着这条路径修改流量值(实际修改的是残量网络的边权),顺带修改反图的残余网络的边权。当没有增广路时,算法停止,此时的流就是最大流。

    。还有这里之所以要建反图是因为,bfs找最短的增广路的时候是随机找的,找的可能不是最’优‘的,那么就会影响最后我们能得到的最大流,这个时候反图就派上用场了,它可以弥补我们之前找最短增广路时犯下的错误,也就是说即使我们之前找的增广路有些可能不是最优的,但是没关系,有反图在,我们最后依然可以得到最大流。

      推荐几篇入门博客:

    1.大佬的博客

    2.大佬的博客

    3.大佬的博客

    #include<stdio.h>//ek算法
    #include<queue>
    #include<string.h>
    #include<iostream>
    using namespace std;
    #define maxn=300;
    #define inf=0x3f3f3f3f;
    int m,n;
    int mp[maxn][maxn];//用于保存正图和反图
    int maxflow;
    int flow[maxn];//flow表示找增广路的过程中的节点的流入值
    bool vis[maxn];
    int father[maxn];//用于回溯建反图用
    void init()
    {
        memset(mp,0,sizeof(mp));
        maxflow=0;
    }
    void solve(int s,int e)
    {
        queue<int>q;
        while(1)
        {
            memset(vis,false,sizeof(vis));
            while(!q.empty()) q.pop();
            memset(flow,0,sizeof(flow));
            flow[s]=inf;
            q.push(s);
            vis[s]=true;
            while(!q.empty())
            {
                int u=q.front();
                q.pop();
                for(int i=s;i<=e;i++)
                {
                    if(!vis[i]&&mp[u][i]>0)
                    {
                        vis[i]=true;
                        q.push(i);
                        flow[i]=min(flow[u],mp[u][i]);
                        father[i]=u;
                    }
                }
                if(flow[e]>0)
                {
                    for(int i=e;i!=1;i=father[i])
                    {
                        mp[father[i]][i]-=flow[e];
                        mp[i][father[i]]+=flow[e];
                    }
                    break;
                }
            }
            if(flow[e]==0) break;
            maxflow+=flow[e];
        }
    }
    int main()
    {
        while(scanf("%d%d",&m,&n)!=EOF)
        {
            init();
            int u,v,w;
            for(int i=1;i<=m;i++)
            {
                scanf("%d%d%d",&u,&v,&w);
                mp[u][v]+=w;
            }
            solve(1,n);
            printf("%d
    ",maxflow);
        }
        return 0;
    }
    
  • 相关阅读:
    ES6模块开发+单文件组件
    Vue路由学习
    Vuex学习
    Vue组件
    Vue事件处理
    Git下载前后端代码步骤
    小黑记事本
    简单计算器
    ubuntu的基础命令
    拓扑排序以及求解关键路径
  • 原文地址:https://www.cnblogs.com/eason9906/p/11755040.html
Copyright © 2020-2023  润新知