最小费用最大流——ZKW
对于最小费用最大流,我们的通常做法是EK+SPFA。
然而,卡常界大佬ZKW发明了一个求解最小费用最大流的方法,很强啊。
在学ZKW费用流前,先说说KM算法。
KM算法
为啥要先提这个呢?因为ZKW费用流用了一个和它非常类似的做法。
KM算法求的是二分图最大权完美匹配。
在此,我来口胡一下(这个算法其实我并未打过,只懂思想)。
和匈牙利算法差不多,区别在于标号。
对于左右两边的点各自有个标号D" role="presentation">DD,一开始,左边的标号设为最大的连出去的边的权值,右边的设为0。
每次只走,满足Di+Dj=v" role="presentation">Di+Dj=vDi+Dj=v的边。
当发生冲突时,就将右边的点的标号全部加一,左边的点的标号全部减一,然后继续搞。
这就是KM算法的大概思路。
ZKW费用流
到正题了。
类似于KM算法,ZKW算法也有个标号。
当Dx=Dy+w" role="presentation">Dx=Dy+wDx=Dy+w时,才会走这条边。
ZKW算法中,是类似于Dinic的多路增广,在每次增广后,就修改距离标号。
求出距离标号的最小增加值,然后修改距离标号。
当找不出增广路时,算法结束。
代码
using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
int n,m,S,T;
struct EDGE
{
int to,c,w;
EDGE *las;
} e[100001];
int ne;
EDGE *last[5001];
#define rev(ei) (e+(int((ei)-e)^1))
int maxflow,mincost;
bool vis[5001];
int dis[5001];
int dfs(int,int);
bool change();
void flow();
int main()
{
scanf("%d%d%d%d",&n,&m,&S,&T);
for (int i=1;i<=m;++i)
{
int u,v,c,w;
scanf("%d%d%d%d",&u,&v,&c,&w);
e[ne]={v,c,w,last[u]};
last[u]=e+ne;
++ne;
e[ne]={u,0,-w,last[v]};
last[v]=e+ne;
++ne;
}
flow();
printf("%d %d\n",maxflow,mincost);
return 0;
}
int dfs(int x,int s)
{
if (x==T)
{
maxflow+=s;
mincost+=dis[S]*s;
return s;
}
vis[x]=1;
int have=0;
for (EDGE *ei=last[x];ei;ei=ei->las)
if (!vis[ei->to] && ei->c && dis[ei->to]+ei->w==dis[x])
{
int t=dfs(ei->to,min(s-have,ei->c));
ei->c-=t;
rev(ei)->c+=t;
have+=t;
}
return have;
}
bool change()
{
int d=INT_MAX;
for (int i=1;i<=n;++i)
if (vis[i])
for (EDGE *ei=last[i];ei;ei=ei->las)
if (!vis[ei->to] && ei->c)
d=min(d,dis[ei->to]-dis[i]+ei->w);//找出最小增加值
if (d==INT_MAX)
return 0;
for (int i=1;i<=n;++i)
if (vis[i])
dis[i]+=d;//更改距离标号
return 1;
}
void flow()
{
maxflow=0;
mincost=0;
do
do
memset(vis,0,sizeof vis);
while (dfs(S,INT_MAX));
while (change());
}
最小费用最大流——ZKW的更多相关文章
- 洛谷.3381.[模板]最小费用最大流(zkw)
题目链接 Update:我好像刚知道多路增广就是zkw费用流.. //1314ms 2.66MB 本题优化明显 #include <queue> #include <cstdio&g ...
- 【BZOJ-3876】支线剧情 有上下界的网络流(有下界有源有汇最小费用最大流)
3876: [Ahoi2014]支线剧情 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 821 Solved: 502[Submit][Status ...
- BZOJ-1061 志愿者招募 线性规划转最小费用最大流+数学模型 建模
本来一眼建模,以为傻逼题,然后发现自己傻逼...根本没想到神奇的数学模型..... 1061: [Noi2008]志愿者招募 Time Limit: 20 Sec Memory Limit: 162 ...
- BZOJ-1070 修车 最小费用最大流+拆点+略坑建图
1070: [SCOI2007]修车 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 3624 Solved: 1452 [Submit][Status] ...
- 【BZOJ】1221: [HNOI2001] 软件开发(最小费用最大流)
http://www.lydsy.com/JudgeOnline/problem.php?id=1221 先吐槽一下,数组依旧开小了RE:在spfa中用了memset和<queue>的版本 ...
- bzoj 1927 [Sdoi2010]星际竞速(最小费用最大流)
1927: [Sdoi2010]星际竞速 Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 1576 Solved: 954[Submit][Statu ...
- Luogu P3381 (模板题) 最小费用最大流
<题目链接> 题目大意: 给定一张图,给定条边的容量和单位流量费用,并且给定源点和汇点.问你从源点到汇点的最带流和在流量最大的情况下的最小费用. 解题分析: 最小费用最大流果题. 下面的是 ...
- POJ - 2516 Minimum Cost(最小费用最大流)
1.K种物品,M个供应商,N个收购商.每种物品从一个供应商运送到一个收购商有一个单位运费.每个收购商都需要K种物品中的若干.求满足所有收购商需求的前提下的最小运费. 2.K种物品拆开来,分别对每种物品 ...
- LOJ#3097 [SNOI2019]通信 最小费用最大流+cdq分治/主席树/分块优化建图
瞎扯 我们网络流模拟赛(其实是数据结构模拟赛)的T2. 考场上写主席树写自闭了,直接交了\(80pts\)的暴力,考完出来突然发现: woc这个题一个cdq几行就搞定了! 题意简述 有\(n\)个哨站 ...
随机推荐
- error LNK2001: unresolved external symbol _main解决办法(zz)
error LNK2001: unresolved external symbol _main解决办法 解决外部符号错误:_main,_WinMain@16,__beginthreadex -!t ...
- 运算符的基本概念以及常用Scanner、随机数Random、选择结构的初步了解
运算符 分类 算术运算符 位运算符 关系运算符|比较运算符 逻辑运算符 条件运算符 赋值运算符 其中优先级顺序从上到下,可以记忆口诀:单目乘除位关系,逻辑三目后赋值 操作数: 运算符左右两边的数 表达 ...
- java_网络编程之BS(web案例)
package BsServersocket; import java.io.*; import java.net.ServerSocket; import java.net.Socket; publ ...
- Activiti学习笔记2 — HelloWorld
一. Maven的POM配置文件 <project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="htt ...
- springcloud Finchley 版本hystrix 和 hystrix Dashboard
hystrix的断路功能 引用上个项目,创建新的model ,cloud-hystrix pom.xml <?xml version="1.0" encoding=" ...
- 关于Button控件的CommandName属性用法的一个实例
注:本文分享于悠闲的博客,地址:http://www.cnblogs.com/9999/archive/2009/11/24/1609234.html 1.前台的代码 <%@ Page Lang ...
- uoj37 主旋律
题意:一个班级n个人,如果a爱b,那么a->b一条有向边.问有多少种删边集合使得图仍然强联通? n<=15. 标程: #include<cstdio> #include&l ...
- jmeter做bbs作业时提示404错误
在用jemter做bbs作业时候,注册成功后会跳转到主页,在写主页的脚本的时候,将fiddler抓到的url复制到路径下方“/bbs/forum.php”,但是第一次复制的时候在/bbs/forum. ...
- Python中死锁的形成示例及死锁情况的防止
死锁示例搞多线程的经常会遇到死锁的问题,学习操作系统的时候会讲到死锁相关的东西,我们用Python直观的演示一下.死锁的一个原因是互斥锁.假设银行系统中,用户a试图转账100块给用户b,与此同时用户b ...
- resful风格
package com.atguigu.springboot.controller; import com.atguigu.springboot.dao.DepartmentDao; import c ...