zkw费用流+当前弧优化

 var
o,v:array[..] of boolean;
f,s,d,dis:array[..] of longint;
next,p,c,w:array[-..] of longint;
i,j,k,l,y,t,ss,tt,n,ans,imp,flow:longint;
procedure link(i,j,k,l:longint);
begin
inc(t);
next[t]:=d[i]; d[i]:=t; p[t]:=j; c[t]:=k; w[t]:=l;
next[-t]:=d[j]; d[j]:=-t; p[-t]:=i; w[-t]:=-l;
end;
function dfs(i,flow:longint):longint;
var j,k,l,min:longint;
begin
if i=tt then
begin
inc(ans,dis[i]*flow); exit(flow);
end;
k:=s[i]; j:=p[k]; dfs:=;
o[i]:=true; v[i]:=true;
while k<> do
begin
l:=dis[i]+w[k]-dis[j]; min:=flow;
if c[k]<min then min:=c[k];
if(min>)and(l<f[j])then f[j]:=l;
if(min>)and(l=)and(not o[j])then
begin
l:=dfs(j,min);
inc(dfs,l); dec(flow,l);
dec(c[k],l); inc(c[-k],l);
end;
if flow= then break;
s[i]:=next[s[i]];
k:=s[i]; j:=p[k];
end;
o[i]:=false;
end;
begin
build;
repeat
for i:=ss to tt do s[i]:=d[i];
fillchar(v,sizeof(v),false);
fillchar(f,sizeof(f),);
inc(flow,dfs(ss, shl ));
imp:= shl ;
for i:=ss to tt do
if(not v[i])and(f[i]<imp)then imp:=f[i];
for i:=ss to tt do if not v[i] then inc(dis[i],imp);
until imp= shl ;
writeln(ans);
end.

PASCAL

 #include<bits/stdc++.h>
using namespace std;
const int INF=;
int o[],v[],f[],c[],s[],dis[];
int b[],p[],d[],w[];
int j,k,ss,tt,n,ans,imp,flow,cnt;
void link(int i,int j,int k,int l)
{
cnt++; b[cnt]=c[i]; c[i]=cnt; p[cnt]=j; d[cnt]=k; w[cnt]=l;
cnt++; b[cnt]=c[j]; c[j]=cnt; p[cnt]=i; d[cnt]=; w[cnt]=-l;
}
int dfs(int x,int flow)
{
if(x==tt){ ans+=dis[x]*flow; return flow; }
int ans2=; o[x]=; v[x]=;
for(int i=s[x];i;i=b[i])
{
int j=p[i]; int l=dis[x]+w[i]-dis[j]; int mi=min(flow,d[i]);
if((mi>)and(l<f[j]))f[j]=l;
if((mi>)and(l==)and(!o[j]))
{
l=dfs(j,mi);
ans2+=l; flow-=l; d[i]-=l; d[i^]+=l;
}
if(flow==)break; s[x]=b[s[x]];
}
o[x]=; return ans2;
}
int main()
{
cnt=; build;
while()
{
for(int i=ss;i<=tt;i++)s[i]=c[i];
memset(v,,sizeof(v)); for(int i=ss;i<=tt;i++)f[i]=INF;
flow+=dfs(ss,INF); imp=INF;
for(int i=ss;i<=tt;i++)if((!v[i])and(f[i]<imp))imp=f[i];
if(imp==INF)break;
for(int i=ss;i<=tt;i++)if(!v[i])dis[i]+=imp;
}
printf("%d %d\n",flow,ans);
}

C++

zkw费用流+当前弧优化的更多相关文章

  1. 图论-zkw费用流

    图论-zkw费用流 模板 这是一个求最小费用最大流的算法,因为发明者是神仙zkw,所以叫zkw费用流(就是zkw线段树那个zkw).有些时候比EK快,有些时候慢一些,没有比普通费用流算法更难,所以学z ...

  2. 学习了ZKW费用流

    所谓ZKW费用流,其实就是Dinic. 若干年前有一个人发明了最小增广路算法,每次用BFS找一条增广路,时间O(nm^2) 然后被DinicD飞了:我们为什么不可以在长度不变时多路增广呢?时间O(n^ ...

  3. zkw费用流

    期末结束,竞赛生活继续开始,先怒刷完寒假作业再说 至于期末考试,数学跪惨,各种哦智障错,还有我初中常用的建系大法居然被自己抛至脑后,看来学的还是不扎实,以后数学要老老实实学.物理被永哥黑了两分,然后很 ...

  4. 【zkw费用流】[网络流24题]餐巾计划问题

    题目描述 一个餐厅在相继的N天里,第i天需要Ri块餐巾(i=l,2,-,N).餐厅可以从三种途径获得餐巾. (1)购买新的餐巾,每块需p分: (2)把用过的餐巾送到快洗部,洗一块需m天,费用需f分(f ...

  5. CSU 1948: 超级管理员(普通费用流&&zkw费用流)

    Description 长者对小明施加了膜法,使得小明每天起床就像马丁的早晨一样. 今天小明早上醒来发现自己成了一位仓管员.仓库可以被描述为一个n × m的网格,在每个网格上有几个箱子(可能没有).为 ...

  6. BZOJ2673 [Wf2011]Chips Challenge 费用流 zkw费用流 网络流

    https://darkbzoj.cf/problem/2673 有一个芯片,芯片上有N*N(1≤N≤40)个插槽,可以在里面装零件. 有些插槽不能装零件,有些插槽必须装零件,剩下的插槽随意. 要求装 ...

  7. zkw费用流 学习笔记

    分析 记\(D_i\)为从\(S\)出发到\(i\)的最短路 最短路算法保证, 算法结束时 对于任意存在弧\((i,j)\)满足\(D_i + c_{ij}\ge D_j\) ① 且对于每个 \(j\ ...

  8. P4015 运输问题【zkw费用流】

    输入输出样例 输入 #1复制 2 3 220 280 170 120 210 77 39 105 150 186 122 输出 #1复制 48500 69140zuixiaofeiyo 说明/提示 1 ...

  9. 最大流当前弧优化Dinic模板

    最大流模板: 普通最大流 无向图限制:将无向图的边拆成2条方向相反的边 无源汇点有最小流限制的最大流:理解为水管流量形成循环,每根水管有流量限制,并且流入量等于流出量 有源汇点的最小流限制的最大流 顶 ...

随机推荐

  1. Ubuntu下安装Naginx, PHP5(及PHP-FPM),MySQL

    一:安装前做个简单的说明 二:安装MySQL 三:安装Nginx 四:安装PHP5 五:配置 nginx,以下是我本机的配置文件. 六:让MySQL支持PHP5 七:配置PHP-FPM 八:在/etc ...

  2. PHP学习之输出语句、注释、算数运算符

    今天学习了PHP的输出语句:

  3. iOS开发和localStorage/sessionStorage

    一.前言 在近期的工作中,有前端同学告诉我要清除localStorage,我当时对localStorage完全没有概念,所以就在w3c看了一下相关的内容,下面简单的介绍一下.算是对iOS开发者普及H5 ...

  4. zendstudio文件编码修改问题

    转载:http://blog.csdn.net/kunlong0909/article/details/7818620 朋友,在zendstudio ide中,你是否碰到导入一个项目后,发现项目中文件 ...

  5. 使用sencha cmd 一键生成你的应用程序代码

    一键生成你的应用程序代码: ------------------------------------------------------------ 我们的出发点就是使用命令来产生一个应用程序,执行以 ...

  6. win7+iis7.5+asp.net下 CS0016: 未能写入输出文件“c:\Windows\Microsoft.NET\Framework\v2.0.50727\Temporary ASP.NET Files 解决方案

    http://www.cnblogs.com/finesite/archive/2011/01/28/1946940.html 给C:\Windows\temp 文件夹赋予完全控制权限,如图:

  7. git技巧记录--blame

    git blame [-L<m,n>] FilePath 可以查看代码每一行是谁写的(根据该行最后一次改动情况), -L表示要查看的行数范围, m: 起始行数, n:结束行数. 方便快速定 ...

  8. Sublime更换默认字体的方法

    Sublime是一款很不错的编辑器,不过默认安装后的字体却不尽人意,并且Sublime竟然连个完整的设置页面都没有(直接让你编辑配置文件).于是很多人对这字体就忍气吞声了.其实只要添加一行代码就可以完 ...

  9. Linux下防火墙开启相关端口及查看已开启端口

    最近利用Apache Mina实现了一个http服务,发布到linux下发现无法访问,通过HttpClient来发送http请求时,报如下错误: Exception in thread "m ...

  10. SQL Server 2008 R2——VC++ ADO 操作 重复利用_ParameterPtr

    ==================================声明================================== 本文原创,转载在正文中显要的注明作者和出处,并保证文章的完 ...