ISAP 简介

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

 1 #include<iostream>
 2 #include<algorithm>
 3 #include<cstring>
 4 #include<cstdlib>
 5 #include<cstdio>
 6
 7 using namespace std;
 8 #define e edges[i]
 9
10 const int Maxn=210,Maxm=210,INF=0x7f7f7f7f;
11 int n,m,s,t;
12 int d[Maxn],cur[Maxn],num[Maxn],p[Maxn];
13 int q[Maxn],ql,qr;
14
15 struct Edge{
16     int to,cap,flow,next;
17     int adv(){return cap-flow;}
18 }edges[Maxm*2];int tot=1,fir[Maxn];
19 void AddEdge(int from,int to,int cap){
20     edges[++tot]=(Edge){to,cap,0,fir[from]};fir[from]=tot;
21     edges[++tot]=(Edge){from,0,0,fir[to]};fir[to]=tot;
22 }
23 void init(){
24     scanf("%d%d",&m,&n);s=1;t=n;
25     for(int u,v,w,i=1;i<=m;i++){
26         scanf("%d%d%d",&u,&v,&w);
27         AddEdge(u,v,w);
28     }
29 }
30 void BFS(){
31     for(int i=1;i<=n;i++)d[i]=n;//这里最好不要用-1或者0x7fffffff等值,否则在统计num的时候可能会出问题
32     q[qr=(ql=0)+1]=t;d[t]=0;
33     for(int x;ql<qr;){
34         x=q[++ql];
35         for(int i=fir[x];i;i=e.next){
36             if(!e.adv()&&d[e.to]==n){//注意这里的!e.adv(),也可以写为!e.cap表示只走反向边,如果不写应该不会影响正确性,只会影响效率
37                 d[e.to]=d[x]+1;
38                 q[++qr]=e.to;
39             }
40         }
41     }
42 }
43 int Augment(){
44     int a=INF;
45     for(int x=t;x!=s;x=edges[p[x]^1].to){
46         a=min(a,edges[p[x]].adv());
47     }
48     for(int x=t;x!=s;x=edges[p[x]^1].to){
49         edges[p[x]].flow+=a;
50         edges[p[x]^1].flow-=a;
51     }
52     return a;
53 }
54 int Maxflow(){
55     memcpy(cur,fir,sizeof cur);
56     BFS();
57     int flow=0;
58     for(int i=1;i<=n;i++)num[d[i]]++;//如果d初始为-1等,访问d[-1],d[0x7fffffff]可能会导致未知错误
59     for(int x=s;d[s]<n;){
60         if(x==t){
61             flow+=Augment();
62             x=s;
63         }
64         int ok=0;
65         for(int&i=cur[x];i;i=e.next){
66             if(e.adv()&&d[e.to]+1==d[x]){
67                 p[x=e.to]=i;
68                 ok=1;
69                 break;
70             }
71         }
72         if(!ok){
73             int M=n;
74             cur[x]=fir[x];
75             for(int i=cur[x];i;i=e.next){
76                 if(e.adv())M=min(M,d[e.to]+1);
77             }
78             if(!--num[d[x]])break;
79             num[d[x]=M]++;
80             if(x!=s)x=edges[p[x]^1].to;
81         }
82     }
83     return flow;
84 }
85 int main(){
86     freopen("input.txt","r",stdin);
87     freopen("","w",stdout);
88
89     init();
90     cout<<Maxflow();
91
92     return 0;
93 }

原始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来记录当前走到哪条边了,下次不走以前的了,接着往下走。

时间: 2024-10-13 11:59:06

ISAP 简介的相关文章

Android网络通讯简介

网络通信应该包含三部分的内容:发送方.接收方.协议栈.发送方和接收方是参与通信的主体,协议栈是发送方和接收方进行通信的契约.按照服务类型,网络通信可分为面向连接和无连接的方式.面向连接是在通信前建立通信链路,而通信结束后释放该链路.无连接的方式则不需要在通信前建立通信连接,这种方式不保证传输的质量. Android提供了多种网络通信的方式,如Java中提供的网络编程,在Android中都提供了支持.Android中常用的网络编程方式如下: 针对TCP/IP协议的Socket和ServerSock

微信红包的架构设计简介

@来源于QCon某高可用架构群整理,整理朱玉华. 背景:有某个朋友在朋友圈咨询微信红包的架构,于是乎有了下面的文字(有误请提出,谢谢) 概况:2014年微信红包使用数据库硬抗整个流量,2015年使用cache抗流量. 微信的金额什么时候算? 答:微信金额是拆的时候实时算出来,不是预先分配的,采用的是纯内存计算,不需要预算空间存储.. 采取实时计算金额的考虑:预算需要占存储,实时效率很高,预算才效率低. 实时性:为什么明明抢到红包,点开后发现没有? 答:2014年的红包一点开就知道金额,分两次操作

JSON 简介

ylbtech-JSON: JSON 简介 JSON:JavaScript Object Notation(JavaScript 对象表示法) JSON是存储和交换文本信息的语法,类似 XML. JSON 比 XML 更小.更快.更易解析. JSON 实例 { "employee":[ {"firstName":"John","lastName":"Doe"}, {"firstName"

Docker简介

Docker简介 什么是Docker: 正所谓Docker的英文本意为"搬运工",所以在我们的世界里,可以理解为Docker搬运的是装满任意类型的APP的集装箱,开发者可以通过Docker将APP变成一种标准化的.可移动植的.自动管理的组件.它用一种新的方式实现了轻量级的虚拟机,专业术语成为应用容器(Application Container) Docker的优势: 1.利用率高 ·Docker对系统资源的利用率很高,一台主机可以同时运行数千个Docker容器 2.可以快速的交付应用程

kafka入门:简介、使用场景、设计原理、主要配置及集群搭建(转)

问题导读: 1.zookeeper在kafka的作用是什么? 2.kafka中几乎不允许对消息进行"随机读写"的原因是什么? 3.kafka集群consumer和producer状态信息是如何保存的? 4.partitions设计的目的的根本原因是什么? 一.入门 1.简介 Kafka is a distributed,partitioned,replicated commit logservice.它提供了类似于JMS的特性,但是在设计实现上完全不同,此外它并不是JMS规范的实现.k

Quartz.NET简介及入门指南

Quartz.NET简介 Quartz.NET是一个功能完备的开源调度系统,从最小的应用到大规模的企业系统皆可适用. Quartz.NET是一个纯净的用C#语言编写的.NET类库,是对非常流行的JAVA开源调度框架 Quartz 的移植. 入门指南 本入门指南包括以下内容: 下载 Quartz.NET 安装 Quartz.NET 根据你的特定项目配置 Quartz 启动一个样例程序 下载和安装 你可以下载 zip 文件或使用 Nuget 程序包.Nuget 程序包只包含 Quartz.NET 运

ASP.Net简介、IIS服务器和Repeater重复器

简介:ASP.NET - 制作网站应用程序的技术 WebForm -出来时间比较早,敏捷.便捷开发,封装一些控件,慢慢发现一些控件做的挺好,真正使用没有那么敏捷 MVC -出来时间比较晚 什么东西? winform 界面 - 后台 - 数据库 共同组合出来的程序:ASP.NET 界面(HTML+CSS+JS) - 后台 - 数据库 运行机制:winform - 程序是安装在用户的电脑上,程序是运行在用户电脑上的.net Framework框架上的 ASP.NET - 通过浏览器向服务器发送请求,

CloudFoundry in 1 Box简介:PCF-Dev篇

在<CloudFoundry in 1 Box简介:Bosh-lite篇>我们介绍了Bosh-lite的架构和部署.在本篇中,我们将详细描述另一个CloudFoundry in 1 Box解决方案PCF-Dev. 1PCF-dev简介 PCF是Pivotal发行的Cloud Foundry商业版,PCF-Dev原名MicroPCF,是Pivotal为PCF的应用开发人员准备的一款App单虚拟机版的CloudFoundry.但是,麻雀虽小,五脏俱全.PCF-Dev虽然可以在仅仅一台虚拟式上即可运

1、elasticsearch简介

1.elasticsearch简介 中文帮助文档地址:http://es.xiaoleilu.com/ • Elasticsearch是一个基于Lucene的实时的分布式搜索和分析引擎.设计用于云计算中,能够达到实时搜索,稳定,可靠,快速,安装使用方便.基于RESTful接口. – 普通请求是...get?a=1 – rest请求....get/a/1 • Elasticsearch的用户 – GitHub,Wikipedia,ebay等... • ES VS Solr – 接口 • 类似web