刘汝佳的蓝书上已经给出了大部分,先给上完整代码(以草地排水为例)。

 #include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cstdio> using namespace std;
#define e edges[i] const int Maxn=,Maxm=,INF=0x7f7f7f7f;
int n,m,s,t;
int d[Maxn],cur[Maxn],num[Maxn],p[Maxn];
int q[Maxn],ql,qr; struct Edge{
int to,cap,flow,next;
int adv(){return cap-flow;}
}edges[Maxm*];int tot=,fir[Maxn];
void AddEdge(int from,int to,int cap){
edges[++tot]=(Edge){to,cap,,fir[from]};fir[from]=tot;
edges[++tot]=(Edge){from,,,fir[to]};fir[to]=tot;
}
void init(){
scanf("%d%d",&m,&n);s=;t=n;
for(int u,v,w,i=;i<=m;i++){
scanf("%d%d%d",&u,&v,&w);
AddEdge(u,v,w);
}
}
void BFS(){
for(int i=;i<=n;i++)d[i]=n;//这里最好不要用-1或者0x7fffffff等值,否则在统计num的时候可能会出问题
q[qr=(ql=)+]=t;d[t]=;
for(int x;ql<qr;){
x=q[++ql];
for(int i=fir[x];i;i=e.next){
if(!e.adv()&&d[e.to]==n){//注意这里的!e.adv(),也可以写为!e.cap表示只走反向边,如果不写应该不会影响正确性,只会影响效率
d[e.to]=d[x]+;
q[++qr]=e.to;
}
}
}
}
int Augment(){
int a=INF;
for(int x=t;x!=s;x=edges[p[x]^].to){
a=min(a,edges[p[x]].adv());
}
for(int x=t;x!=s;x=edges[p[x]^].to){
edges[p[x]].flow+=a;
edges[p[x]^].flow-=a;
}
return a;
}
int Maxflow(){
memcpy(cur,fir,sizeof cur);
BFS();
int flow=;
for(int i=;i<=n;i++)num[d[i]]++;//如果d初始为-1等,访问d[-1],d[0x7fffffff]可能会导致未知错误
for(int x=s;d[s]<n;){
if(x==t){
flow+=Augment();
x=s;
}
int ok=;
for(int&i=cur[x];i;i=e.next){
if(e.adv()&&d[e.to]+==d[x]){
p[x=e.to]=i;
ok=;
break;
}
}
if(!ok){
int M=n;
cur[x]=fir[x];
for(int i=cur[x];i;i=e.next){
if(e.adv())M=min(M,d[e.to]+);
}
if(!--num[d[x]])break;
num[d[x]=M]++;
if(x!=s)x=edges[p[x]^].to;
}
}
return flow;
}
int main(){
freopen("input.txt","r",stdin);
freopen("","w",stdout); init();
cout<<Maxflow(); return ;
}

原始EK的效率浪费在每次都要一遍BFS,而Dinic选择多路增广来优化效率,ISAP则还是单路增广,但是只用做一次BFS。

这里的BFS是反向的,表示到汇点的距离(当然正向也没什么问题),当我们流完一条增广路,有些点的d值可能会改变,这时我们要修改这些点的d值,把他设为min{d[v]+1,(u,v)没有满流}(u是这个点),这样的正确性是显然的。然后我们只需要每次只走距离相差为1的点就行了QuQ.

还有就是退出条件,当d[s]>=n时,显然不存在增广路了,因为每走一个点,到汇点的距离会-1,这样假设能走到汇点,那么d[t]>=1,是矛盾的。

还有,如果某个d值没有任何一个点是,那么就会出现断层,此时显然也走不到汇点,这就是所谓的gap优化?

还有个当前弧优化,就是用cur来记录当前走到哪条边了,下次不走以前的了,接着往下走。

ISAP 简介的更多相关文章

  1. ASP.NET Core 1.1 简介

    ASP.NET Core 1.1 于2016年11月16日发布.这个版本包括许多伟大的新功能以及许多错误修复和一般的增强.这个版本包含了多个新的中间件组件.针对Windows的WebListener服 ...

  2. MVVM模式和在WPF中的实现(一)MVVM模式简介

    MVVM模式解析和在WPF中的实现(一) MVVM模式简介 系列目录: MVVM模式解析和在WPF中的实现(一)MVVM模式简介 MVVM模式解析和在WPF中的实现(二)数据绑定 MVVM模式解析和在 ...

  3. Cassandra简介

    在前面的一篇文章<图形数据库Neo4J简介>中,我们介绍了一种非常流行的图形数据库Neo4J的使用方法.而在本文中,我们将对另外一种类型的NoSQL数据库——Cassandra进行简单地介 ...

  4. REST简介

    一说到REST,我想大家的第一反应就是“啊,就是那种前后台通信方式.”但是在要求详细讲述它所提出的各个约束,以及如何开始搭建REST服务时,却很少有人能够清晰地说出它到底是什么,需要遵守什么样的准则. ...

  5. Microservice架构模式简介

    在2014年,Sam Newman,Martin Fowler在ThoughtWorks的一位同事,出版了一本新书<Building Microservices>.该书描述了如何按照Mic ...

  6. const,static,extern 简介

    const,static,extern 简介 一.const与宏的区别: const简介:之前常用的字符串常量,一般是抽成宏,但是苹果不推荐我们抽成宏,推荐我们使用const常量. 执行时刻:宏是预编 ...

  7. HTTPS简介

    一.简单总结 1.HTTPS概念总结 HTTPS 就是对HTTP进行了TLS或SSL加密. 应用层的HTTP协议通过传输层的TCP协议来传输,HTTPS 在 HTTP和 TCP中间加了一层TLS/SS ...

  8. 【Machine Learning】机器学习及其基础概念简介

    机器学习及其基础概念简介 作者:白宁超 2016年12月23日21:24:51 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本系列文章是作者结 ...

  9. Cesium简介以及离线部署运行

    Cesium简介 cesium是国外一个基于JavaScript编写的使用WebGL的地图引擎,一款开源3DGIS的js库.cesium支持3D,2D,2.5D形式的地图展示,可以自行绘制图形,高亮区 ...

随机推荐

  1. hdoj 1251 字典树

    代码: #include <stdio.h>#define  MAX    26 typedef struct TrieNode{     int nCount;      struct ...

  2. 关于call 与 apply 那些事

    1.定义: call : 调用一个对象的一个方法,以另一个对象替换当前的对象. apply : 应用某一对象的一个方法,用另一个对象替换当前的对象. 2.用法: call:call(thisObj, ...

  3. C、C++、Java异或运算交换变量变量值的区别

    今天看到一位大神的博客,深受感触.决定也发一篇博客,证明一下我还活着. 于是我翻看以前学习时做的一些笔记,整理了一下,得到了一个关于异或运算交换变量变量值的笔记. 首先来看下面三组表达式,看起来他们都 ...

  4. Android Fragment基础及使用

    同一个app内的界面切换 用Fragment比较合适,因为Activity比较重量级 Fragment 轻量级,切换灵活 --------------------------------------- ...

  5. JQuery树形目录制作

    <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1//EN" "http://www.w3.org/TR/xhtml11/DT ...

  6. MFC之MessageBox用法

    一    函数原型及参数 function MessageBox(hWnd: HWND; Text, Caption: PChar; Type: Word): Integer; hWnd:对话框父窗口 ...

  7. Web NFC API

    W3C Editor's Draft 29 December 2014 This version: http://www.w3.org/2012/nfc/web-api/ Latest publish ...

  8. Asynchronous

    http://sporto.github.io/blog/2012/12/09/callbacks-listeners-promises/ http://www.cnblogs.com/liuning ...

  9. Delphi调用WINAPI时到底应该是指针还是结构体(注意是Delphi变量本身就是指针)

    看MSDN,GetWindowRect的说明如下: BOOL WINAPI GetWindowRect( _In_  HWND   hWnd, _Out_ LPRECT lpRect // 注意,没* ...

  10. php pack、unpack、ord 函数使用方法

    string pack ( string $format [, mixed $args [, mixed $... ]] ) Pack given arguments into a binary st ...