index.html 37 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697
  1. <!DOCTYPE html><html lang="en"><head><meta charset="UTF-8"><meta name="viewport" content="width=device-width,initial-scale=1,maximum-scale=2"><meta name="theme-color" content="#222"><meta name="generator" content="Hexo 4.2.0"><link rel="apple-touch-icon" sizes="180x180" href="/blog/blog/images/apple-touch-icon-next.png"><link rel="icon" type="image/png" sizes="32x32" href="/blog/blog/images/favicon-frog.png"><link rel="icon" type="image/png" sizes="16x16" href="/blog/blog/images/favicon-frog.png"><link rel="mask-icon" href="/blog/blog/images/logo.svg" color="#222"><link rel="stylesheet" href="/blog/css/main.css"><link rel="stylesheet" href="//fonts.googleapis.com/css?family=Comic Sans MS:300,300italic,400,400italic,700,700italic|Consolas:300,300italic,400,400italic,700,700italic&display=swap&subset=latin,latin-ext"><link rel="stylesheet" href="/blog/lib/font-awesome/css/font-awesome.min.css"><link rel="stylesheet" href="/blog/lib/pace/pace-theme-minimal.min.css"><script src="/blog/lib/pace/pace.min.js"></script><script id="hexo-configurations">var NexT=window.NexT||{},CONFIG={hostname:"schtonn.github.io",root:"/blog/",scheme:"Muse",version:"7.8.0",exturl:!1,sidebar:{position:"left",display:"post",padding:18,offset:12,onmobile:!1},copycode:{enable:!0,show_result:!0,style:"flat"},back2top:{enable:!0,sidebar:!1,scrollpercent:!0},bookmark:{enable:!1,color:"#222",save:"auto"},fancybox:!1,mediumzoom:!1,lazyload:!0,pangu:!1,comments:{style:"tabs",active:"valine",storage:!0,lazyload:!1,nav:null,activeClass:"valine"},algolia:{hits:{per_page:10},labels:{input_placeholder:"Search for Posts",hits_empty:"We didn't find any results for the search: ${query}",hits_stats:"${hits} results found in ${time} ms"}},localsearch:{enable:!1,trigger:"auto",top_n_per_article:1,unescape:!1,preload:!1},motion:{enable:!0,async:!1,transition:{post_block:"fadeIn",post_header:"slideDownIn",post_body:"slideDownIn",coll_header:"slideLeftIn",sidebar:"slideUpIn"}}}</script><meta name="description" content="前置知识 存图方式(邻接表,邻接矩阵),图的遍历(dfs,bfs) 引入 我们举个例子吧: 有一个毒瘤水管工,他会造水管,有一天他造了一个水管网络,就像一个图。其中有一个点只有出边,就是源点(起点),还有一个点只有入边,是汇点(终点)。 点之间有一些管子,管子 \(i\) 在一个单位时间内能流 \(c_i\) 的水,现在毒瘤水管工想知道,他的水管在一个单位时间内,从源点最多能运多少水到汇点"><meta property="og:type" content="article"><meta property="og:title" content="网络最大流-Dinic"><meta property="og:url" content="https://schtonn.github.io/blog/posts/dinic/index.html"><meta property="og:site_name" content="schtonn"><meta property="og:description" content="前置知识 存图方式(邻接表,邻接矩阵),图的遍历(dfs,bfs) 引入 我们举个例子吧: 有一个毒瘤水管工,他会造水管,有一天他造了一个水管网络,就像一个图。其中有一个点只有出边,就是源点(起点),还有一个点只有入边,是汇点(终点)。 点之间有一些管子,管子 \(i\) 在一个单位时间内能流 \(c_i\) 的水,现在毒瘤水管工想知道,他的水管在一个单位时间内,从源点最多能运多少水到汇点"><meta property="og:locale" content="en_US"><meta property="og:image" content="https://img-blog.csdnimg.cn/20200111112801923.png"><meta property="og:image" content="https://img-blog.csdnimg.cn/20200111114148208.png"><meta property="og:image" content="https://img-blog.csdnimg.cn/20200111112801923.png"><meta property="og:image" content="https://img-blog.csdnimg.cn/20200111114953988.png"><meta property="og:image" content="https://schtonn.github.io/blog/images/tumor_fixed.png"><meta property="article:published_time" content="2020-03-02T03:31:22.000Z"><meta property="article:modified_time" content="2022-10-19T15:02:06.678Z"><meta property="article:author" content="Alex"><meta property="article:tag" content="graph"><meta name="twitter:card" content="summary"><meta name="twitter:image" content="https://img-blog.csdnimg.cn/20200111112801923.png"><link rel="canonical" href="https://schtonn.github.io/blog/posts/dinic/"><script id="page-configurations">CONFIG.page={sidebar:"",isHome:!1,isPost:!0,lang:"en"}</script><title>网络最大流-Dinic | schtonn</title><noscript><style>.sidebar-inner,.use-motion .brand,.use-motion .collection-header,.use-motion .comments,.use-motion .menu-item,.use-motion .pagination,.use-motion .post-block,.use-motion .post-body,.use-motion .post-header{opacity:initial}.use-motion .site-subtitle,.use-motion .site-title{opacity:initial;top:initial}.use-motion .logo-line-before i{left:initial}.use-motion .logo-line-after i{right:initial}</style></noscript></head><body itemscope itemtype="http://schema.org/WebPage"><div class="container use-motion"><div class="headband"></div><header class="header" itemscope itemtype="http://schema.org/WPHeader"><div class="header-inner"><div class="site-brand-container"><div class="site-nav-toggle"><div class="toggle" aria-label="Toggle navigation bar"><span class="toggle-line toggle-line-first"></span><span class="toggle-line toggle-line-middle"></span><span class="toggle-line toggle-line-last"></span></div></div><div class="site-meta"><a href="/blog/" class="brand" rel="start"><span class="logo-line-before"><i></i></span><h1 class="site-title">schtonn</h1><span class="logo-line-after"><i></i></span></a><p class="site-subtitle" itemprop="description">schtonn</p></div><div class="site-nav-right"><div class="toggle popup-trigger"></div></div></div><nav class="site-nav"><ul id="menu" class="menu"><li class="menu-item menu-item-home"><a href="/blog/" rel="section"><i class="fa fa-fw fa-home"></i> Home</a></li><li class="menu-item menu-item-tags"><a href="/blog/tags/" rel="section"><i class="fa fa-fw fa-tags"></i> Tags</a></li><li class="menu-item menu-item-archives"><a href="/blog/archives/" rel="section"><i class="fa fa-fw fa-archive"></i> Archives</a></li><li class="menu-item menu-item-games"><a href="/blog/games/" rel="section"><i class="fa fa-fw fa-gamepad"></i> Games</a></li></ul></nav></div></header><div class="back-to-top"><i class="fa fa-arrow-up"></i> <span>0%</span></div><main class="main"><div class="main-inner"><div class="content-wrap"><div class="content post posts-expand"><article itemscope itemtype="http://schema.org/Article" class="post-block" lang="en"><link itemprop="mainEntityOfPage" href="https://schtonn.github.io/blog/posts/dinic/"><span hidden itemprop="author" itemscope itemtype="http://schema.org/Person"><meta itemprop="image" content="/blog/images/avatar.gif"><meta itemprop="name" content="Alex"><meta itemprop="description" content="blog"></span><script type="text/javascript" src="/blog/js/md5.js"></script><script></script><script>document.oncopy=function(e){window.event&&(e=window.event);try{var t=e.srcElement;return"INPUT"==t.tagName&&"text"==t.type.toLowerCase()||"TEXTAREA"==t.tagName}catch(e){return!1}}</script><span hidden itemprop="publisher" itemscope itemtype="http://schema.org/Organization"><meta itemprop="name" content="schtonn"></span><header class="post-header"><h1 class="post-title" itemprop="name headline"> 网络最大流-Dinic</h1><div class="post-meta"><span class="post-meta-item"><span class="post-meta-item-icon"><i class="fa fa-calendar-o"></i></span> <span class="post-meta-item-text">Posted on</span> <time title="Created: 2020-Mar-02 11:31:22" itemprop="dateCreated datePublished" datetime="2020-03-02T11:31:22+08:00">2020-Mar-02</time></span><span class="post-meta-item"><span class="post-meta-item-icon"><i class="fa fa-calendar-check-o"></i></span> <span class="post-meta-item-text">Edited on</span> <time title="Modified: 2022-Oct-19 23:02:06" itemprop="dateModified" datetime="2022-10-19T23:02:06+08:00">2022-Oct-19</time></span><span class="post-meta-item"><span class="post-meta-item-icon"><i class="fa fa-comment-o"></i></span> <span class="post-meta-item-text">Valine:</span><a title="valine" href="/blog/posts/dinic/#valine-comments" itemprop="discussionUrl"><span class="post-comments-count valine-comment-count" data-xid="/blog/posts/dinic/" itemprop="commentCount"></span></a></span></div></header><div class="post-body" itemprop="articleBody"><h2 id="前置知识">前置知识</h2><p>存图方式(邻接表,邻接矩阵),图的遍历(dfs,bfs)</p><h2 id="引入">引入</h2><p>我们举个例子吧: 有一个毒瘤水管工,他会造水管,有一天他造了一个水管网络,就像一个图。其中有一个点只有出边,就是源点(起点),还有一个点只有入边,是汇点(终点)。</p><p>点之间有一些管子,管子 <span class="math inline">\(i\)</span> 在一个单位时间内能流 <span class="math inline">\(c_i\)</span> 的水,现在毒瘤水管工想知道,他的水管在一个单位时间内,从源点最多能运多少水到汇点。<a id="more"></a></p><h2 id="概念">概念</h2><ul><li>源点:起点,一般用s表示;</li><li>汇点:终点,一般用t表示;</li><li>残余网络:进行增广后剩余的网络;</li><li>増广: &gt; 増广就是在残余网络中寻找从源点到汇点的可行路径(増广路),并将该路径上的所有边的容量减去路径中的最小容量,形成新的残余网络,<strong>人话就是找一条能走的路,看看这条路上最窄的水管是那条,然后让最窄的水管满载运行。</strong> 例如: <img data-src="https://img-blog.csdnimg.cn/20200111112801923.png"> 如果当前有这样一个残余网络,那么 <span class="math inline">\(s\rightarrow4\rightarrow1\rightarrow t\)</span> 就是一条増广路,最小容量是4,进行増广过后就形成了这样一张图: <img data-src="https://img-blog.csdnimg.cn/20200111114148208.png"> 如何寻找増广路?直接dfs即可。</li><li>反向边: &gt; 有时候,増广的时会出现假的答案,例如还是那个图: <img data-src="https://img-blog.csdnimg.cn/20200111112801923.png"> 有 <span class="math inline">\(s\rightarrow4\rightarrow1\rightarrow t\)</span> 和 <span class="math inline">\(s\rightarrow1\rightarrow t\)</span> 两条増广路,万一程序选错了怎么办? 这时就要请出反向边了。 每次増广的时候,<strong>在残余网络上逆着増广路径建容量与増广路径容量相等的反向边</strong>,比如刚才那张图,就顺着<span class="math inline">\(t\rightarrow1\rightarrow4\rightarrow s\)</span>建容量为4的边。如果接下来増广时走过了反向边,就相当于把原来的増广撤回去了。 这就相当于给了程序一个反悔的机会。</li></ul><h2 id="朴素算法">朴素算法</h2><ol type="1"><li>増广</li><li>沿着増广路径建反向边</li><li>如果源点与汇点依然连通,回到1</li></ol><h2 id="dinic优化">Dinic优化</h2><p>Dinic的优化就是用bfs建立了由s开始的一个分层图,每次寻找増广路时必须让边上的层数严格递增,就可以确保每一步都离汇点近了一些这样就不会陷入毒瘤数据卡成的死循环,比如这样的著名毒瘤数据:</p><p><img data-src="https://img-blog.csdnimg.cn/20200111114953988.png"></p><p>在这个数据中,如果用朴素算法,就会让中间容量为1的边上下抖动抽搐,等到他抽了999次的时候才把上面和下面的999减没。如果用Dinic,则不会走中间的边,因为层数都是1。所以两次直接求出999+999。</p><p><img data-src="/blog/images/tumor_fixed.png"></p><h2 id="代码">代码</h2><h3 id="朴素">朴素</h3><figure class="highlight cpp"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br><span class="line">35</span><br><span class="line">36</span><br><span class="line">37</span><br><span class="line">38</span><br><span class="line">39</span><br><span class="line">40</span><br><span class="line">41</span><br><span class="line">42</span><br><span class="line">43</span><br><span class="line">44</span><br><span class="line">45</span><br><span class="line">46</span><br><span class="line">47</span><br><span class="line">48</span><br><span class="line">49</span><br><span class="line">50</span><br><span class="line">51</span><br></pre></td><td class="code"><pre><span class="line"><span class="meta">#<span class="meta-keyword">include</span> <span class="meta-string">"bits/stdc++.h"</span></span></span><br><span class="line"><span class="keyword">using</span> <span class="keyword">namespace</span> <span class="built_in">std</span>;</span><br><span class="line"><span class="keyword">const</span> <span class="keyword">long</span> <span class="keyword">long</span> N=<span class="number">10010</span>,inf=<span class="number">0x3f3f3f3f3f3f3f3f</span>;</span><br><span class="line"><span class="class"><span class="keyword">struct</span> <span class="title">node</span>&#123;</span></span><br><span class="line"> <span class="keyword">long</span> <span class="keyword">long</span> v,c,nxt;</span><br><span class="line">&#125;e[N&lt;&lt;<span class="number">1</span>];</span><br><span class="line"><span class="keyword">long</span> <span class="keyword">long</span> tot=<span class="number">1</span>,h[N];</span><br><span class="line"><span class="function"><span class="keyword">void</span> <span class="title">adde</span><span class="params">(<span class="keyword">long</span> <span class="keyword">long</span> u,<span class="keyword">long</span> <span class="keyword">long</span> v,<span class="keyword">long</span> <span class="keyword">long</span> c)</span></span>&#123;</span><br><span class="line"> e[++tot].v=v;</span><br><span class="line"> e[tot].c=c;</span><br><span class="line"> e[tot].nxt=h[u];</span><br><span class="line"> h[u]=tot;</span><br><span class="line">&#125;</span><br><span class="line"><span class="keyword">long</span> <span class="keyword">long</span> n,m,ans,vis[N];</span><br><span class="line"><span class="function"><span class="keyword">long</span> <span class="keyword">long</span> <span class="title">path</span><span class="params">(<span class="keyword">long</span> <span class="keyword">long</span> u,<span class="keyword">long</span> <span class="keyword">long</span> f)</span></span>&#123;</span><br><span class="line"> <span class="keyword">if</span>(vis[u])<span class="keyword">return</span> <span class="number">0</span>;</span><br><span class="line"> vis[u]=<span class="literal">true</span>;</span><br><span class="line"><span class="comment">// cout&lt;&lt;u&lt;&lt;' '&lt;&lt;f&lt;&lt;endl;</span></span><br><span class="line"> <span class="keyword">if</span>(u==n)<span class="keyword">return</span> f;</span><br><span class="line"> <span class="keyword">for</span>(<span class="keyword">long</span> <span class="keyword">long</span> i=h[u];i;i=e[i].nxt)&#123;</span><br><span class="line"> <span class="keyword">long</span> <span class="keyword">long</span> v=e[i].v;</span><br><span class="line"> <span class="keyword">if</span>(e[i].c)&#123;</span><br><span class="line"> <span class="keyword">long</span> <span class="keyword">long</span> fv=path(v,min(f,e[i].c));</span><br><span class="line"> <span class="keyword">if</span>(fv)&#123;</span><br><span class="line"> e[i].c-=fv;</span><br><span class="line"> e[i^<span class="number">1</span>].c+=fv;</span><br><span class="line"> <span class="keyword">return</span> fv;</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> <span class="keyword">return</span> <span class="number">0</span>;</span><br><span class="line">&#125;</span><br><span class="line"><span class="function"><span class="keyword">void</span> <span class="title">solve</span><span class="params">()</span></span>&#123;</span><br><span class="line"> <span class="keyword">long</span> <span class="keyword">long</span> d=<span class="number">0</span>;</span><br><span class="line"> <span class="keyword">while</span>(d=path(<span class="number">1</span>,inf))&#123;</span><br><span class="line"> ans+=d;</span><br><span class="line"> <span class="built_in">memset</span>(vis,<span class="number">0</span>,<span class="keyword">sizeof</span>(vis));</span><br><span class="line"> &#125;</span><br><span class="line">&#125;</span><br><span class="line"><span class="function"><span class="keyword">int</span> <span class="title">main</span><span class="params">()</span></span>&#123;</span><br><span class="line"> <span class="built_in">cin</span>&gt;&gt;n&gt;&gt;m;</span><br><span class="line"> <span class="keyword">for</span>(<span class="keyword">long</span> <span class="keyword">long</span> i=<span class="number">1</span>;i&lt;=m;i++)&#123;</span><br><span class="line"> <span class="keyword">long</span> <span class="keyword">long</span> u,v,c;</span><br><span class="line"> <span class="built_in">cin</span>&gt;&gt;u&gt;&gt;v&gt;&gt;c;</span><br><span class="line"> adde(u,v,c);</span><br><span class="line"> adde(v,u,<span class="number">0</span>);</span><br><span class="line"> &#125;</span><br><span class="line"> solve();</span><br><span class="line"> <span class="built_in">cout</span>&lt;&lt;ans&lt;&lt;<span class="built_in">endl</span>;</span><br><span class="line"> <span class="keyword">return</span> <span class="number">0</span>;</span><br><span class="line">&#125;</span><br></pre></td></tr></table></figure><h3 id="带优化">带优化</h3><figure class="highlight cpp"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br><span class="line">35</span><br><span class="line">36</span><br><span class="line">37</span><br><span class="line">38</span><br><span class="line">39</span><br><span class="line">40</span><br><span class="line">41</span><br><span class="line">42</span><br><span class="line">43</span><br><span class="line">44</span><br><span class="line">45</span><br><span class="line">46</span><br><span class="line">47</span><br><span class="line">48</span><br><span class="line">49</span><br><span class="line">50</span><br><span class="line">51</span><br><span class="line">52</span><br><span class="line">53</span><br><span class="line">54</span><br><span class="line">55</span><br><span class="line">56</span><br><span class="line">57</span><br><span class="line">58</span><br><span class="line">59</span><br><span class="line">60</span><br><span class="line">61</span><br><span class="line">62</span><br><span class="line">63</span><br><span class="line">64</span><br><span class="line">65</span><br><span class="line">66</span><br><span class="line">67</span><br><span class="line">68</span><br><span class="line">69</span><br><span class="line">70</span><br><span class="line">71</span><br><span class="line">72</span><br><span class="line">73</span><br><span class="line">74</span><br><span class="line">75</span><br></pre></td><td class="code"><pre><span class="line"><span class="meta">#<span class="meta-keyword">include</span> <span class="meta-string">"iostream"</span></span></span><br><span class="line"><span class="meta">#<span class="meta-keyword">include</span> <span class="meta-string">"cstring"</span></span></span><br><span class="line"><span class="meta">#<span class="meta-keyword">include</span> <span class="meta-string">"queue"</span></span></span><br><span class="line"><span class="keyword">using</span> <span class="keyword">namespace</span> <span class="built_in">std</span>;</span><br><span class="line"><span class="keyword">const</span> <span class="keyword">int</span> maxn=<span class="number">1e6</span>+<span class="number">5</span>;</span><br><span class="line"><span class="keyword">const</span> <span class="keyword">int</span> INF=<span class="number">9999999</span>;</span><br><span class="line"><span class="keyword">int</span> s,t,tot,n,m;<span class="comment">//s为源点,t为汇点 </span></span><br><span class="line"><span class="class"><span class="keyword">struct</span> <span class="title">node</span>&#123;</span></span><br><span class="line"> <span class="keyword">int</span> next;</span><br><span class="line"> <span class="keyword">int</span> v;</span><br><span class="line"> <span class="keyword">int</span> c;</span><br><span class="line">&#125;e[maxn];</span><br><span class="line"><span class="keyword">int</span> h[maxn],q[maxn],deth[maxn];</span><br><span class="line"><span class="function"><span class="keyword">void</span> <span class="title">add</span><span class="params">(<span class="keyword">int</span> u,<span class="keyword">int</span> v,<span class="keyword">int</span> c)</span></span>&#123;</span><br><span class="line"> e[tot].v=v;</span><br><span class="line"> e[tot].next=h[u];</span><br><span class="line"> e[tot].c=c;</span><br><span class="line"> h[u]=tot++;</span><br><span class="line">&#125;</span><br><span class="line"><span class="function"><span class="keyword">int</span> <span class="title">bfs</span><span class="params">()</span></span>&#123;</span><br><span class="line"> <span class="built_in">memset</span>(deth,<span class="number">0</span>,<span class="keyword">sizeof</span>(deth));<span class="comment">//deth记录深度 </span></span><br><span class="line"> <span class="built_in">queue</span>&lt;<span class="keyword">int</span>&gt; q;</span><br><span class="line"> <span class="keyword">while</span>(!q.empty())q.pop();</span><br><span class="line"> deth[s]=<span class="number">1</span>;</span><br><span class="line"> q.push(s);</span><br><span class="line"> <span class="keyword">do</span>&#123;</span><br><span class="line"> <span class="keyword">int</span> u=q.front();</span><br><span class="line"> q.pop();</span><br><span class="line"> <span class="keyword">for</span>(<span class="keyword">int</span> i=h[u];i;i=e[i].next)&#123;</span><br><span class="line"> <span class="keyword">if</span>(e[i].c&gt;<span class="number">0</span>&amp;&amp;deth[e[i].v]==<span class="number">0</span>)&#123;</span><br><span class="line"> deth[e[i].v]=deth[u]+<span class="number">1</span>;</span><br><span class="line"> q.push(e[i].v);</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> <span class="keyword">while</span>(!q.empty());</span><br><span class="line"> <span class="keyword">if</span>(deth[t]!=<span class="number">0</span>)<span class="comment">//当汇点的深度不存在时,说明不存在分层图,同时也说明不存在增广路</span></span><br><span class="line"> <span class="keyword">return</span> <span class="number">1</span>;</span><br><span class="line"> <span class="keyword">else</span> <span class="keyword">return</span> <span class="number">0</span>;</span><br><span class="line">&#125;</span><br><span class="line"><span class="function"><span class="keyword">int</span> <span class="title">dfs</span><span class="params">(<span class="keyword">int</span> u,<span class="keyword">int</span> dist)</span><span class="comment">//u是当前节点,dist为当前流量 </span></span></span><br><span class="line"><span class="function"></span>&#123;</span><br><span class="line"> <span class="keyword">if</span>(u==t)<span class="keyword">return</span> dist;</span><br><span class="line"> <span class="keyword">for</span>(<span class="keyword">int</span> i=h[u];i;i=e[i].next)&#123;</span><br><span class="line"> <span class="keyword">if</span>((deth[e[i].v]==deth[u]+<span class="number">1</span>)&amp;&amp;(e[i].c!=<span class="number">0</span>))&#123;</span><br><span class="line"> <span class="keyword">int</span> di=dfs(e[i].v,min(dist,e[i].c)); </span><br><span class="line"> <span class="keyword">if</span>(di&gt;<span class="number">0</span>)&#123;</span><br><span class="line"> e[i].c-=di;</span><br><span class="line"> e[i^<span class="number">1</span>].c+=di;</span><br><span class="line"> <span class="keyword">return</span> di;</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> <span class="keyword">return</span> <span class="number">0</span>;<span class="comment">//否则说明没有增广路,返回0</span></span><br><span class="line">&#125;</span><br><span class="line"><span class="function"><span class="keyword">int</span> <span class="title">dinic</span><span class="params">()</span></span>&#123;</span><br><span class="line"> <span class="keyword">int</span> ans=<span class="number">0</span>;</span><br><span class="line"> <span class="keyword">while</span>(bfs())&#123;</span><br><span class="line"> <span class="keyword">while</span>(<span class="keyword">int</span> di=dfs(s,INF))&#123;</span><br><span class="line"> ans+=di;</span><br><span class="line"> &#125;</span><br><span class="line"> &#125;</span><br><span class="line"> <span class="keyword">return</span> ans;</span><br><span class="line">&#125;</span><br><span class="line"><span class="function"><span class="keyword">int</span> <span class="title">main</span><span class="params">()</span></span>&#123;</span><br><span class="line"> <span class="built_in">cin</span>&gt;&gt;n&gt;&gt;m&gt;&gt;s&gt;&gt;t;</span><br><span class="line"> <span class="keyword">int</span> x,y,z;</span><br><span class="line"> <span class="keyword">for</span>(<span class="keyword">int</span> i=<span class="number">0</span>;i&lt;m;i++)&#123;</span><br><span class="line"> <span class="built_in">cin</span>&gt;&gt;x&gt;&gt;y&gt;&gt;z;</span><br><span class="line"> add(x,y,z);</span><br><span class="line"> add(y,x,<span class="number">0</span>);</span><br><span class="line"> &#125;</span><br><span class="line"> <span class="built_in">cout</span>&lt;&lt;dinic()&lt;&lt;<span class="built_in">endl</span>;</span><br><span class="line"> <span class="keyword">return</span> <span class="number">0</span>;</span><br><span class="line">&#125;</span><br></pre></td></tr></table></figure><h2 id="我踩的坑">我踩的坑</h2><p>注意当实现朴素算法时,代码中的 dfs 不能回溯,只需要一条路走到底。</p></div><div><ul class="post-copyright"><li class="post-copyright-author"> <strong>Post author:</strong> Alex</li><li class="post-copyright-link"> <strong>Post link:</strong> <a href="https://schtonn.github.io/blog/posts/dinic/" title="网络最大流-Dinic">https://schtonn.github.io/blog/posts/dinic/</a></li><li class="post-copyright-license"> <strong>Copyright Notice:</strong> All articles in this blog are licensed under<a href="https://creativecommons.org/licenses/by-nc-sa/4.0/" rel="noopener" target="_blank"><i class="fa fa-fw fa-creative-commons"></i> BY-NC-SA</a> unless stating additionally.</li></ul></div><footer class="post-footer"><div class="post-tags"><a href="/blog/tags/graph/" rel="tag"><i class="fa fa-tag"></i> graph</a></div><div class="post-nav"><div class="post-nav-item"><a href="/blog/posts/min-span-tree/" rel="prev" title="最小生成树"><i class="fa fa-chevron-left"></i> 最小生成树</a></div><div class="post-nav-item"> <a href="/blog/posts/fibonacci/" rel="next" title="斐波那契数列">斐波那契数列<i class="fa fa-chevron-right"></i></a></div></div></footer></article></div><div class="comments" id="valine-comments"></div><script>
  2. window.addEventListener('tabs:register', () => {
  3. let { activeClass } = CONFIG.comments;
  4. if (CONFIG.comments.storage) {
  5. activeClass = localStorage.getItem('comments_active') || activeClass;
  6. }
  7. if (activeClass) {
  8. let activeTab = document.querySelector(`a[href="#comment-${activeClass}"]`);
  9. if (activeTab) {
  10. activeTab.click();
  11. }
  12. }
  13. });
  14. if (CONFIG.comments.storage) {
  15. window.addEventListener('tabs:click', event => {
  16. if (!event.target.matches('.tabs-comment .tab-content .tab-pane')) return;
  17. let commentClass = event.target.classList[1];
  18. localStorage.setItem('comments_active', commentClass);
  19. });
  20. }
  21. </script></div><div class="toggle sidebar-toggle"><span class="toggle-line toggle-line-first"></span><span class="toggle-line toggle-line-middle"></span><span class="toggle-line toggle-line-last"></span></div><aside class="sidebar"><div class="sidebar-inner"><ul class="sidebar-nav motion-element"><li class="sidebar-nav-toc"> Table of Contents</li><li class="sidebar-nav-overview"> Overview</li></ul><div class="post-toc-wrap sidebar-panel"><div class="post-toc motion-element"><ol class="nav"><li class="nav-item nav-level-2"><a class="nav-link" href="#前置知识"><span class="nav-number">1.</span> <span class="nav-text">前置知识</span></a></li><li class="nav-item nav-level-2"><a class="nav-link" href="#引入"><span class="nav-number">2.</span> <span class="nav-text">引入</span></a></li><li class="nav-item nav-level-2"><a class="nav-link" href="#概念"><span class="nav-number">3.</span> <span class="nav-text">概念</span></a></li><li class="nav-item nav-level-2"><a class="nav-link" href="#朴素算法"><span class="nav-number">4.</span> <span class="nav-text">朴素算法</span></a></li><li class="nav-item nav-level-2"><a class="nav-link" href="#dinic优化"><span class="nav-number">5.</span> <span class="nav-text">Dinic优化</span></a></li><li class="nav-item nav-level-2"><a class="nav-link" href="#代码"><span class="nav-number">6.</span> <span class="nav-text">代码</span></a><ol class="nav-child"><li class="nav-item nav-level-3"><a class="nav-link" href="#朴素"><span class="nav-number">6.1.</span> <span class="nav-text">朴素</span></a></li><li class="nav-item nav-level-3"><a class="nav-link" href="#带优化"><span class="nav-number">6.2.</span> <span class="nav-text">带优化</span></a></li></ol></li><li class="nav-item nav-level-2"><a class="nav-link" href="#我踩的坑"><span class="nav-number">7.</span> <span class="nav-text">我踩的坑</span></a></li></ol></div></div><div class="site-overview-wrap sidebar-panel"><div class="site-author motion-element" itemprop="author" itemscope itemtype="http://schema.org/Person"><p class="site-author-name" itemprop="name">Alex</p><div class="site-description" itemprop="description">blog</div></div><div class="site-state-wrap motion-element"><nav class="site-state"><div class="site-state-item site-state-posts"> <a href="/blog/archives"><span class="site-state-item-count">35</span> <span class="site-state-item-name">posts</span></a></div><div class="site-state-item site-state-tags"> <a href="/blog/tags/"><span class="site-state-item-count">8</span> <span class="site-state-item-name">tags</span></a></div></nav></div><div class="links-of-author motion-element"><span class="links-of-author-item"><a href="https://github.com/schtonn" title="GitHub → https:&#x2F;&#x2F;github.com&#x2F;schtonn" rel="noopener" target="_blank"><i class="fa fa-fw fa-github"></i> GitHub</a></span><span class="links-of-author-item"><a href="mailto:schtonn@163.com" title="E-Mail → mailto:schtonn@163.com" rel="noopener" target="_blank"><i class="fa fa-fw fa-envelope"></i> E-Mail</a></span></div><div class="links-of-blogroll motion-element"><div class="links-of-blogroll-title"><i class="fa fa-fw fa-link"></i> Links</div><ul class="links-of-blogroll-list"><li class="links-of-blogroll-item"> <a href="https://yonghong.github.io/" title="https:&#x2F;&#x2F;yonghong.github.io" rel="noopener" target="_blank">Yonghong</a></li><li class="links-of-blogroll-item"> <a href="https://source.unsplash.com/random/1600x900" title="https:&#x2F;&#x2F;source.unsplash.com&#x2F;random&#x2F;1600x900" rel="noopener" target="_blank">Background</a></li></ul><iframe width="400" height="300" frameborder="0" src="https://cdn.abowman.com/widgets/treefrog/index.html?up_bodyColor=2d2d2d&up_pattern=0&up_patternColor=000000&up_footColor=2d2d2d&up_eyeColor=3a3a3a&up_bellySize=50&up_backgroundColor=222222&up_tongueColor=2b2d2d&up_flyColor=3a3a3a&up_releaseFly=0"></iframe></div></div></div></aside><div id="sidebar-dimmer"></div></div></main><footer class="footer"><div class="footer-inner"><div class="copyright"> &copy; 2019 – <span itemprop="copyrightYear">2023</span><span class="with-love"><i class="fa fa-user"></i></span> <span class="author" itemprop="copyrightHolder">Alexander</span></div></div></footer></div><script src="/blog/lib/anime.min.js"></script><script src="//cdn.jsdelivr.net/npm/lozad@1/dist/lozad.min.js"></script><script src="/blog/lib/velocity/velocity.min.js"></script><script src="/blog/lib/velocity/velocity.ui.min.js"></script><script src="/blog/js/utils.js"></script><script src="/blog/js/motion.js"></script><script src="/blog/js/schemes/muse.js"></script><script src="/blog/js/next-boot.js"></script><script>!function(){var t=document.createElement("script"),e=window.location.protocol.split(":")[0];t.src="https"===e?"https://zz.bdstatic.com/linksubmit/push.js":"http://push.zhanzhang.baidu.com/push.js";var s=document.getElementsByTagName("script")[0];s.parentNode.insertBefore(t,s)}()</script><script>
  22. if (typeof MathJax === 'undefined') {
  23. window.MathJax = {
  24. loader: {
  25. load: ['[tex]/mhchem'],
  26. source: {
  27. '[tex]/amsCd': '[tex]/amscd',
  28. '[tex]/AMScd': '[tex]/amscd'
  29. }
  30. },
  31. tex: {
  32. inlineMath: {'[+]': [['$', '$']]},
  33. packages: {'[+]': ['mhchem']},
  34. tags: 'ams'
  35. },
  36. options: {
  37. renderActions: {
  38. findScript: [10, doc => {
  39. document.querySelectorAll('script[type^="math/tex"]').forEach(node => {
  40. const display = !!node.type.match(/; *mode=display/);
  41. const math = new doc.options.MathItem(node.textContent, doc.inputJax[0], display);
  42. const text = document.createTextNode('');
  43. node.parentNode.replaceChild(text, node);
  44. math.start = {node: text, delim: '', n: 0};
  45. math.end = {node: text, delim: '', n: 0};
  46. doc.math.push(math);
  47. });
  48. }, '', false],
  49. insertedScript: [200, () => {
  50. document.querySelectorAll('mjx-container').forEach(node => {
  51. let target = node.parentNode;
  52. if (target.nodeName.toLowerCase() === 'li') {
  53. target.parentNode.classList.add('has-jax');
  54. }
  55. });
  56. }, '', false]
  57. }
  58. }
  59. };
  60. (function () {
  61. var script = document.createElement('script');
  62. script.src = '//cdn.jsdelivr.net/npm/mathjax@3/es5/tex-mml-chtml.js';
  63. script.defer = true;
  64. document.head.appendChild(script);
  65. })();
  66. } else {
  67. MathJax.startup.document.state(0);
  68. MathJax.texReset();
  69. MathJax.typeset();
  70. }
  71. </script><script>
  72. NexT.utils.loadComments(document.querySelector('#valine-comments'), () => {
  73. NexT.utils.getScript('https://cdn.jsdelivr.net/npm/valine@1/dist/Valine.min.js', () => {
  74. var GUEST = ['nick', 'mail', 'link'];
  75. var guest = 'nick,mail';
  76. guest = guest.split(',').filter(item => {
  77. return GUEST.includes(item);
  78. });
  79. new Valine({
  80. el : '#valine-comments',
  81. verify : false,
  82. notify : false,
  83. appId : 'BmologYYnRqCv0SLHDeDdA17-gzGzoHsz',
  84. appKey : 'w9mVebFMdCmY6Nh9vfcBGaGt',
  85. placeholder: "Comment...",
  86. avatar : 'mp',
  87. meta : guest,
  88. pageSize : '10' || 10,
  89. visitor : false,
  90. lang : 'en' || 'zh-cn',
  91. path : location.pathname,
  92. recordIP : true,
  93. serverURLs : ''
  94. });
  95. }, window.Valine);
  96. });
  97. </script></body></html>