hdu 1532 Drainage Ditches

最新推荐文章于 2020-01-27 13:13:08 发布

--子非鱼--

最新推荐文章于 2020-01-27 13:13:08 发布

阅读量744

点赞数

CC 4.0 BY-SA版权

分类专栏：网络流编程

本文链接：https://blog.youkuaiyun.com/hp_satan/article/details/8868801

编程同时被 2 个专栏收录

71 篇文章

订阅专栏

网络流

3 篇文章

订阅专栏

本文深入探讨了最大流问题的概念及其在计算机科学中的应用，特别关注使用SAP算法解决裸最大流问题的过程。通过实例代码，详细解释了如何实现并优化算法以求得最优解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

点击打开题目链接

裸的最大流问题，用的是SAP算法。

#include<iostream>
#include<cstdio>
#include<cstring>
#define clr(a,b) memset(a,b,sizeof(a))
#define M 0xfffff
using namespace std;
const int I=300;

int map[I][I],pre[I],dis[I],gap[I];

int SAP(int s,int t)
{
    clr(pre,-1);
    clr(dis,0);
    clr(gap,0);
    gap[0]=t;
    int u=s,ans(0);
    while(dis[s]<t)
    {
        int v;
        for(v=1;v<=t;v++)
        {
            if(map[u][v]>0 && dis[u]==dis[v]+1)
                break;
        }
        if(v<=t)
        {
            pre[v]=u;
            u=v;
            if(v==t)
            {
                int minflow=M;
                for(int i=v;i!=s;i=pre[i])
                {
                    if(minflow>map[pre[i]][i])
                        minflow=map[pre[i]][i];
                }
                ans+=minflow;
                for(int i=v;i!=s;i=pre[i])
                {
                    map[pre[i]][i]-=minflow;
                    map[i][pre[i]]+=minflow;
                }
                u=s;
            }
        }
        else
        {
            int mindis=t;
            for(int i=1;i<=t;i++)
            {
                if(map[u][i]>0 && mindis>dis[i]+1)
                    mindis=dis[i]+1;
            }
            --gap[dis[u]];
            if(gap[dis[u]]==0) return ans;
            dis[u]=mindis;
            gap[mindis]++;
            if(u!=s) u=pre[u];
        }
    }
    return ans;
}

int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        clr(map,0);
        for(int i=0;i<n;i++)
        {
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            map[a][b]+=c;
        }
        int ans=SAP(1,m);
        printf("%d\n",ans);
    }
    return 0;
}