当前位置:首页 > 芯闻号 > 充电吧
[导读]膜板题 ->排水沟poj1273 网络流弄了将近一个月,寒假玩v家去了导致信息竞赛进度不大,不过还是先总结一下吧 一些概念 一 网络流满足三个性质: 1 容量限制 每条边不能够提供大

膜板题 ->排水沟poj1273
网络流弄了将近一个月,寒假玩v家去了导致信息竞赛进度不大,不过还是先总结一下吧

一些概念

一 网络流满足三个性质:
1 容量限制
每条边不能够提供大于其流量的边,(反向边要加上对应的增广流量是为了满足流量守恒)
2 对称性
f(u,v)=-f(v,u)
3 流量守恒
所有点的流入量等于流出量
通俗点说就是你不能从1号城市运3个冰球到2号,而运到2号的再传到3号城市的冰球变成了1个(好吧很意义不明)

二 增广路
能走的路一定是还有流量的路,
而增广路是从s->t还能多增加流的路。

三 交替路
二分图从一端到另一端,走一边到另一边再到原来的一边的非同一个顶点的路是交替路(比较通俗)

存图方式

用vector模拟邻接表

struct edge
{
int to,cap,rev;
};

其中cap是剩余容量,rev是反向边(和刘汝佳方法不一样)推荐白书上那种 方便

关于两个模板,fulkerson(我老是打错成**)和dinic算法
我们先来说一下fulkerson算法吧

fulkerson算法 复杂度O(F|E|)

用dfs找增广路 为了不走回头路,我们使用数组标记
一条路能走当且仅当它的剩余容量大于0,那么我们只需要不断地沿着路dfs就行了,唯一注意的是,我们要在增广一条路后去掉对应的容量 并把反向边加上容量(根据容量守恒原则)
没什么很难的地方
代码 可以练好后交poj1273

#include 
#include 
#include 
#include 
#include 
#include 

using namespace std;

const int maxn=500;
const int inf=0x3f3f3f3f;
struct edge
{
    int to;
    int cap,rev;
    void ini()
    {
        to=0;
        cap=0;
        rev=0;
    }
};
vector g[maxn];
bool used[maxn];
int s,t,n,m;
int ans=0;


void addedge(int from,int to,int cap)
{
    g[from].push_back((edge){to,cap,g[to].size()});
    g[to].push_back((edge){from,0,g[from].size()-1});
}

int dfs(int v,int t,int f)
{
    if(v==t) return f;
    used[v]=true;
    for(unsigned i=0;i0)
        {
            int d;
            d=dfs(buffer.to,t,min(f,buffer.cap));
            if(d>0)
            {
                g[v][i].cap-=d;
                g[buffer.to][buffer.rev].cap+=d;
                return d;
            }
        }
    }
    return 0;
}

int max_flow(int st,int ed)
{
    ans=0;
    while(1)
    {
        memset(used,0,sizeof(used));
        int buff=dfs(st,ed,inf);
        if(buff==0) return ans;
        else
        {
            ans+=buff;
        }
    }
}

int main()
{
    while(scanf("%d%d",&m,&n)==2)
    {
        s=1;
        t=n;
        for(int i=0;i
dinic算法 复杂度O(|V||V||E|) 上界比较松

将图分层,
对于一个图,每一次增广一定有一条路会断(即没有剩余容量)
这时我们记一个dis 数组来记起点到各个点的距离(我们暂时只考虑所有cost为1的情况)
这样可以保证不对一条没用的边进行多次检查,节省了时间,其余与深搜类似

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
//手写队列
//AC
using namespace std;

const int maxn=5000000;
const int inf=0x3f3f3f3f;
struct edge
{
    int to;
    int cap;
    int rev;
};
vector g[maxn];
int dis[maxn];
int cur[maxn];
int q[maxn];
int s,t,n,m;
int res=0;

int readint()
{
    char c;
    int ans=0;
    while(c=getchar())
    {
        if(isdigit(c))
        {
            ans=ans*10+c-'0';
        }
        else c=getchar();
    }
    return ans;
}

void addedge(int from,int to,int cap)
{
    g[from].push_back((edge){to,cap,g[to].size()});
    g[to].push_back((edge){from,0,g[from].size()-1});
}

void bfs()
{
    memset(dis,-1,sizeof(dis));
    memset(q,0,sizeof(q));
    int hd=0;
    int tl=0;
    q[hd]=s;
    dis[s]=0;
    while(hd<=tl)
    {
        int head=q[hd];
        hd++;
        for(unsigned i=0;i0 && dis[e.to]<0)
            {
                dis[e.to]=dis[head]+1;
                tl++;
                q[tl]=e.to;
            }
        }
    }
}

int dfs(int v,int ed,int flow)
{
    int f;
    if(v==ed)
    {
        return flow;
    }
    else
    {
        for(int &i=cur[v];i0 && dis[e.to]==dis[v]+1)
            {
                f=dfs(e.to,ed,min(flow,e.cap));
                if(f>0)
                {
                    e.cap-=f;
                    g[e.to][e.rev].cap+=f;//注意结构体后要带成员
                    return f;
                }
            }
        }
    }
    return 0;
}

int max_flow(int st,int ed)
{
    int ans=0;
    while(1)
    {
        bfs();
        if(dis[t]<0) return ans;
        else
        {
            memset(cur,0,sizeof(cur));
            int f;
            while((f=dfs(st,ed,inf))>0)
            {
                ans+=f;
            }
        }
    }
}

int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t);
        res=0;
        for(int i=0;i
二分图匹配的网络流变形

对于一个二分图,我们想到如果将最大匹配转化为最大流就好了。
那么如何转化呢?对于一个匹配,一定是数条没有公共顶点的边,那么就想到了数条这样相同流的增广路
我们添加一个源点和汇点,它们分别向两边的两点各连一条容量为1,cost=0的边,这样就转化为了最大流来求这些边的问题了

#include 
#include 
#include 
#include 
#include 
#include 
//AC
using namespace std;

const int maxn=10000;
const int inf=0x3f3f3f3f;
int dis[maxn];
int cur[maxn];
int s,t,n,m,e;

struct edge
{
    int to,cap,rev;
};
vector g[maxn];

void addedge(int from,int to,int cap)
{
    g[from].push_back((edge){to,cap,g[to].size()});
    g[to].push_back((edge){from,0,g[from].size()-1});
}
//dinic
bool bfs(int st,int ed)
{
    memset(dis,inf,sizeof(dis));
    int q[maxn];
    int head=0;
    int tail=0;
    q[head]=st;
    dis[st]=0;
    while(head<=tail)
    {
        int h=q[head];
        head++;
        for(unsigned i=0;i0 && dis[e.to]==inf)
            {
                dis[e.to]=dis[h]+1;
                tail++;
                q[tail]=e.to;
            }
        }
    }
    return (dis[ed]!=inf);
}

int dfs(int v,int ed,int leftflow)
{
    if(leftflow==0 || v==ed)
    {
        return leftflow;
    }
    else
    {
        for(int &i=cur[v];i0 && dis[e.to]==dis[v]+1)/////
            {
                int f=dfs(e.to,ed,min(leftflow,e.cap));
                if(f>0)
                {
                    e.cap-=f;///////////////////
                    g[e.to][e.rev].cap+=f;///////////////
                    return f;
                }
            }
        }
        return 0;
    }
}

int main()
{
    scanf("%d%d%d",&n,&m,&e);
    s=0;
    t=n+m+1;
    memset(dis,inf,sizeof(dis));
    memset(cur,0,sizeof(cur));
    for(int i=0;i=1 && from<=n && to>=1 && to<=m)
        {
            to+=n;
            addedge(from,to,1);
        }
    }
    for(int i=1;i<=n;i++)
    {
        addedge(s,i,1);
    }
    for(int j=n+1;j<=n+m;j++)
    {
        addedge(j,t,1);
    }
    int ans=0;
    while(bfs(s,t))
    {
        int f;
        memset(cur,0,sizeof(cur));/////////////
        while((f=dfs(s,t,inf))>0)
        {
            ans+=f;
        }
    }
    printf("%dn",ans);
    return 0;
}
固定流最小费用流问题

给你10个冰球,你要运完,但是你要经过长度最小的道路
如何做呢?这里推荐一个用连续最短路的贝尔曼福德实现的算法:
流固定,每次以最短路增广,从流中减去对应流
这样写看上去比较好懂的


int min_cost_flow(int st,int ed,int f)
{
    int ans=0;
    while(f>0)
    {
        bool update=true;
        memset(dis,inf,sizeof(dis));
        dis[st]=0;
        while(update)
        {
            update=false;
            for(int i=0;i<=n+m+1;i++)
            {
                if(dis[i]!=inf)
                {
                    for(unsigned j=0;j0 && dis[e.to]>dis[i]+e.cost)
                        {
                            dis[e.to]=dis[i]+e.cost;
                            prevv[e.to]=i;
                            preve[e.to]=j;
                            update=true;
                        }
                    }
                }
            }
        }

        if(dis[ed]==inf)
        {
            return -1;
        }

        int d=inf;
        for(int i=ed;i!=st;i=prevv[i])
        {
            d=min(d,g[prevv[i]][preve[i]].cap);
        }

        f-=d;
        ans+=dis[ed]*d;
        for(int i=ed;i!=st;i=prevv[i])
        {
            edge & e=g[prevv[i]][preve[i]];
            e.cap-=d;
            g[i][e.rev].cap+=d;
        }
    }
    return ans;
}

应用当然不止这些 最主要是建模,这里有些题(书上的例题)
poj1273排水沟
poj3281奶牛的餐饮
poj3686玩具
luogu3386二分图匹配模板

本站声明: 本文章由作者或相关机构授权发布,目的在于传递更多信息,并不代表本站赞同其观点,本站亦不保证或承诺内容真实性等。需要转载请联系该专栏作者,如若文章内容侵犯您的权益,请及时联系本站删除。
换一批
延伸阅读

9月2日消息,不造车的华为或将催生出更大的独角兽公司,随着阿维塔和赛力斯的入局,华为引望愈发显得引人瞩目。

关键字: 阿维塔 塞力斯 华为

加利福尼亚州圣克拉拉县2024年8月30日 /美通社/ -- 数字化转型技术解决方案公司Trianz今天宣布,该公司与Amazon Web Services (AWS)签订了...

关键字: AWS AN BSP 数字化

伦敦2024年8月29日 /美通社/ -- 英国汽车技术公司SODA.Auto推出其旗舰产品SODA V,这是全球首款涵盖汽车工程师从创意到认证的所有需求的工具,可用于创建软件定义汽车。 SODA V工具的开发耗时1.5...

关键字: 汽车 人工智能 智能驱动 BSP

北京2024年8月28日 /美通社/ -- 越来越多用户希望企业业务能7×24不间断运行,同时企业却面临越来越多业务中断的风险,如企业系统复杂性的增加,频繁的功能更新和发布等。如何确保业务连续性,提升韧性,成...

关键字: 亚马逊 解密 控制平面 BSP

8月30日消息,据媒体报道,腾讯和网易近期正在缩减他们对日本游戏市场的投资。

关键字: 腾讯 编码器 CPU

8月28日消息,今天上午,2024中国国际大数据产业博览会开幕式在贵阳举行,华为董事、质量流程IT总裁陶景文发表了演讲。

关键字: 华为 12nm EDA 半导体

8月28日消息,在2024中国国际大数据产业博览会上,华为常务董事、华为云CEO张平安发表演讲称,数字世界的话语权最终是由生态的繁荣决定的。

关键字: 华为 12nm 手机 卫星通信

要点: 有效应对环境变化,经营业绩稳中有升 落实提质增效举措,毛利润率延续升势 战略布局成效显著,战新业务引领增长 以科技创新为引领,提升企业核心竞争力 坚持高质量发展策略,塑强核心竞争优势...

关键字: 通信 BSP 电信运营商 数字经济

北京2024年8月27日 /美通社/ -- 8月21日,由中央广播电视总台与中国电影电视技术学会联合牵头组建的NVI技术创新联盟在BIRTV2024超高清全产业链发展研讨会上宣布正式成立。 活动现场 NVI技术创新联...

关键字: VI 传输协议 音频 BSP

北京2024年8月27日 /美通社/ -- 在8月23日举办的2024年长三角生态绿色一体化发展示范区联合招商会上,软通动力信息技术(集团)股份有限公司(以下简称"软通动力")与长三角投资(上海)有限...

关键字: BSP 信息技术
关闭
关闭