最大流/ISAP

话说ISAP是真快...(大多数情况)吊打dinic,而且还好写...

大概思路就是: 在dinic的基础上, 动态修改层数, 如果终点层数 \(>\) 点数, break.

暂时并不知道isap有没有可能被卡. 应该不会

const int ninf=(int)1e9+500,psz=(int)450,esz=(int)1e5+50;

//network flow
struct te{int t,pr,fl;}edge[esz*2];
int hd[psz],pe=1,np,ss,tt; //np: number of points
void adde(int f,int t,int fl){edge[++pe]=(te){t,hd[f],fl};hd[f]=pe;}
void addsg(int f,int t,int fl){adde(f,t,fl);adde(t,f,0);} //add an unidirectional edge
void adddb(int f,int t,int fl){adde(f,t,fl);adde(t,f,fl);} //add a bi-directional edge int d[psz],gap[psz],cur[psz];
void init(){ //bfs
rep(i,1,np)cur[i]=hd[i],gap[i]=0,d[i]=0;
static int que[psz],qh,qt;
qh=1,qt=0;
que[++qt]=tt,d[tt]=1,gap[1]=1;
while(qh<=qt){
int u=que[qh++];
for(int i=hd[u],v=edge[i].t;i;i=edge[i].pr,v=edge[i].t){
if(d[v])continue;
d[v]=d[u]+1,++gap[d[v]];
que[++qt]=v;
}
}
}
int dfs(int p,int mi){
if(mi==0||p==tt)return mi;
int fl=0,tmp;
for(int &i=cur[p],v=edge[i].t;i;i=edge[i].pr,v=edge[i].t){
if(d[v]+1!=d[p])continue;
tmp=dfs(v,min(edge[i].fl,mi));
fl+=tmp,mi-=tmp,edge[i].fl-=tmp,edge[i^1].fl+=tmp;
if(mi==0)return fl; //这个时候不需要修改d[p], 已经没有流量可以增广; 如果修改会wa,并不知道为什么,欢迎神仙指出
}
if(gap[d[p]]==1)d[ss]=np+1;
--gap[d[p]],++d[p],++gap[d[p]];
cur[p]=hd[p];
return fl;
}
int maxfl(){
int res=0;
init();
while(d[ss]<=np)res+=dfs(ss,ninf);
return res;
}

费用流/ZKW

跑的慢

const int psz=1e4+50,esz=1e5+50,ninf=1e9;
int n,m,s,t;
struct te{int t,c,v,pr;}edge[esz*2];
int hd[psz],pe=1,ss,tt,np;
void adde(int f,int t,int c,int v){edge[++pe]=(te){t,c,v,hd[f]};hd[f]=pe;}
void addsg(int f,int t,int c,int v){adde(f,t,c,v);adde(t,f,0,-v);} int que[psz],qh=1,qt=0;
int mind[psz],cur[psz],vi[psz];
bool init(){
memcpy(cur,hd,(np+2)*4);
qh=1,qt=0;
rep(i,1,np)mind[i]=ninf,vi[i]=0;
que[++qt]=ss,mind[ss]=0,vi[ss]=1;
while(qh<=qt){
int u=que[qh%np];++qh,vi[u]=0;
for(int i=hd[u],v;i;i=edge[i].pr){
v=edge[i].t;
if(edge[i].c==0||mind[v]<=edge[i].v+mind[u])continue;
mind[v]=edge[i].v+mind[u];
if(vi[v]==0)++qt,que[qt%np]=v,vi[v]=1;
}
}
return mind[tt]<ninf;
}
int dfs(int p,int mi){
if(p==tt||mi==0)return mi;
vi[p]=1;
int fl=0,tmp;
for(int &i=cur[p],v;i;i=edge[i].pr){
v=edge[i].t;
if(mind[v]!=mind[p]+edge[i].v||edge[i].c==0||vi[v])continue;
tmp=dfs(v,min(mi,edge[i].c));
fl+=tmp,mi-=tmp,edge[i].c-=tmp,edge[i^1].c+=tmp;
if(mi==0)break;//important
}
vi[p]=0;return fl;
}
int ansc=0,ansf=0;
void mcmf(){
while(init()){
int tmp=dfs(ss,ninf);
ansf+=tmp,ansc+=tmp*mind[tt];
}
}

[模板] 网络流相关/最大流ISAP/费用流zkw的更多相关文章

  1. 【BZOJ1834】网络扩容(最大流,费用流)

    [BZOJ1834]网络扩容(最大流,费用流) 题面 Description 给定一张有向图,每条边都有一个容量C和一个扩容费用W.这里扩容费用是指将容量扩大1所需的费用.求: 1. 在不扩容的情况下 ...

  2. 洛谷P3381 【模板】最小费用最大流(dijstra费用流)

    题目描述 如题,给出一个网络图,以及其源点和汇点,每条边已知其最大流量和单位流量费用,求出其网络最大流和在最大流情况下的最小费用. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表 ...

  3. 【BZOJ1834】network 网络扩容(最大流,费用流)

    题意:给定一张有向图,每条边都有一个容量C和一个扩容费用W.这里扩容费用是指将容量扩大1所需的费用. 求: 1. 在不扩容的情况下,1到N的最大流: 2. 将1到N的最大流增加K所需的最小扩容费用. ...

  4. BZOJ1834 [ZJOI2010]network 网络扩容 【最大流,费用流】

    1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec  Memory Limit: 64 MB Submit: 3394  Solved: 1774 [Subm ...

  5. luogu P5470 [NOI2019]序列 dp 贪心 费用流 模拟费用流

    LINK:序列 考虑前20分 容易想到爆搜. 考虑dp 容易设\(f_{i,j,k,l}\)表示前i个位置 选了j对 且此时A选择了k个 B选择了l个的最大值.期望得分28. code //#incl ...

  6. 【洛谷 P2604】 [ZJOI2010]网络扩容(最大流,费用流)

    题目链接 第一问就是简单的最大流. 第二问,保留第一问求完最大流的残量网络. 然后新建一个源点,向原源点连一条流量为k,费用为0的边. 然后所有边重新连一起(原来的边保留),费用为题目所给,最小费用即 ...

  7. HDU 5045(Contest-费用流)[template:费用流]

    Contest Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submi ...

  8. hdu1533 费用流模板

    Going Home Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total ...

  9. 【Codevs1237&网络流24题餐巾计划】(费用流)

    题意:一个餐厅在相继的 N 天里,每天需用的餐巾数不尽相同. 假设第 i 天需要 ri块餐巾(i=1,2,…,N).餐厅可以购买新的餐巾,每块餐巾的费用为 p 分: 或者把旧餐巾送到快洗部,洗一块需 ...

随机推荐

  1. vue中的路由的跳转的参数

    vue中的路由跳转传参 params 与 query this.$router.push({ name:"detail", params:{ name:'nameValue', c ...

  2. location.origin兼容IE

    if (window["context"] == undefined) { if (!window.location.origin) { window.location.origi ...

  3. Android为TV端助力 播放视频卡顿问题

    问题分析: 1.连接服务器,ping IP 看看有没有丢包 2.写一个只有mediaplayer的Demo,放到同款的盒子或者电视,保证网咯也是一样的情况下,看看Demo卡不卡 如果一直不卡,那有可能 ...

  4. MUI 实用教程

    MUI 实用JS教程: https://www.kancloud.cn/benhailong/mui/319751  MUI 实用教程: https://www.kancloud.cn/benhail ...

  5. Powershell-获取命令和帮助

    Get-Help 获取命令的帮助文档 Update -Help 更新帮助文档 Save-Help 保存文档 Get-Help Get-VM 加上-Full参数获取详细说明 [-name] <st ...

  6. PowerDesigner导出SQL,注释为空时以name代替

    版本 操作步骤 打开Edit Current DBMS 选中Script->Objects->Column->Add 将Value中的内容全部替换为如下 %20:COLUMN% [% ...

  7. LeetCode算法题-Minimum Absolute Difference in BST(Java实现)

    这是悦乐书的第253次更新,第266篇原创 01 看题和准备 今天介绍的是LeetCode算法题中Easy级别的第120题(顺位题号是530).给定具有非负值的二叉搜索树,找到任意两个节点的值之间的最 ...

  8. IDEA包名显示设置

    项目结构视图右上角那个齿轮 选择[Compact Empty Middle Packages],包会合并显示 [Hide Empty Middle Packages]去掉前面的√,不分层级显示

  9. 一篇文章搞定mongodb

    一 安装 1 安装目录下新建文件夹data,etc,logs #在bin文件下启动cmd,指定数据存储的路径mongod --dbpath D:\MongoDB\data\db 2 etc文件夹中新建 ...

  10. 爬虫系列二(数据清洗--->正则表达式)

    一 正则常识 模式 描述 \w 匹配字母数字及下划线 \W 匹配非字母数字及下划线 \s 匹配任意空白字符,等价于 [\t\n\r\f]. \S 匹配任意非空字符 \d 匹配任意数字,等价于 [0-9 ...