青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

考慮這樣一種網絡流問題:需要對同一個圖求多次最大流。則在每次求最大流之前,需要將所有邊的容量全部恢復到初始值(求最大流的過程中,邊的容量f值被改變了)。不過這還不算最猥瑣的,有的時候,我們需要在每次求最大流之前都刪去圖中的一些點或一些邊,或者改變某些原有的邊的容量,特別是需要刪點或刪邊的情況爆難搞。因為,一般的邊表中邊類型定義如下:
struct edge {
        
int a, b, f, next;
} ed[MAXM 
+ MAXM];
表示這條邊起點為a,終點為b,容量為f,鄰接邊(就是下一條起點為a的邊)的編號為next。
如果要求多次最大流,那么在每次求最大流之前就要把所有邊的容量恢復到初始值,解決方法是引入一個“初始容量”域fs,在這條邊剛剛被加入的時候將它的fs值和f值都設為初始給它的容量,在恢復時,只要將所有邊的f值恢復到fs即可。若要改變邊容量,則將該邊的fs值和f值都設為改變后的容量即可。

下面來分析需要刪邊或刪點的情況。在這種情況下,如果采用只有next的單向鏈表,則刪除時next域極難處理,而且,在一般的邊表中,還設立了表頭hd,hd[a]表示邊表中起點為a的第一條邊的編號。因此,若刪除的邊<a, b, f>是起點為a的第一條邊,還會影響hd[a]的值,使情況變得更為復雜。因此,必須采用雙向鏈表!還記得Dancing Link么?邊表其實也可以搞成Dancing Link,方法如下:
設圖中有N個點,M條邊(注意,這M條邊只包括正向邊,不包括反向邊。由于每條正向邊<a, b, f>都對應一條反向邊<b, a, 0>,因此邊表中邊的數目其實是M+M)。首先把邊表ed的0~N這(N+1)個位置(下標)空出來,作表頭(表頭不是邊,因此在遍歷邊的時候不會遍歷到它們)。其中,ed[0]為總表頭,用于遍歷ed[1..N]中每個未被刪去的點;ed[1..N]為單點表頭,ed[i]用來遍歷圖中所有以i為起點的邊(和DLX中的二維DL驚人相似)。然后,若N為偶數,則空一個位置(也就是將ed[N+1]丟棄不用),這是因為我們在增廣過程中需要引用到一條邊對應的逆向邊(正向邊對應反向邊,反向邊對應正向邊),一般認為編號為p的邊對應的逆向邊是p ^ 1,這樣,就要求圖中所有正向邊的編號都是偶數,所有反向邊的編號都是奇數(否則會造成混亂)。因此當N為偶數時,(N+1)為奇數,不能放置第一條正向邊,需要從ed[N+2]開始放置正向邊。若N為奇數則不用空位。
接下來就是邊類型了。在這里,邊類型一共需要包括6個域:a, b, fs, f, pre, next,表示這條邊起點為a,終點為b,初始容量為fs,當前容量為f,上一條起點為a的邊編號為pre,下一條起點為a的邊編號為next。注意,和DL一樣,整個鏈表是循環的,也就是我們認為表中最后一條起點為a的邊的下一條鄰接邊編號就是a(表頭),同樣,a的上一條鄰接邊也就是這條邊。
struct edge {
    
int a, b, fs, f, pre, next;
} ed[MAXM 
+ MAXM];
接下來就是幾個重要過程了。
(1)初始化表頭:
void init_d()
{
    re1(i, n) {ed[i].a 
= i; ed[i].b = -1; ed[i].f = 0; ed[i].pre = ed[i].next = i;}
    
if (n % 2) m = n + 1else m = n + 2;
}
這里n是圖中的點數(相當于N),m是邊的編號指針(相當于DLX中的nodes)
(2)添加新邊:
void add_edge(int a, int b, int f)
{
    ed[m].a 
= a; ed[m].b = b; ed[m].fs = ed[m].f = f; ed[m].pre = ed[a].pre; ed[m].next = a; ed[a].pre = m; ed[ed[m].pre].next = m++;
    ed[m].a 
= b; ed[m].b = a; ed[m].fs = ed[m].f = 0; ed[m].pre = ed[b].pre; ed[m].next = b; ed[b].pre = m; ed[ed[m].pre].next = m++;
}
這個和DLX類似,不解釋了囧……

最后進入最核心的部分——到底如何處理刪邊或刪點?有了DL型邊表就爆好搞了:刪去一條邊,只要直接刪去該邊在DL中的位置即可:
void deledge(int No)
{
    ed[ed[No].pre].next 
= ed[No].next; ed[ed[No].next].pre = ed[No].pre;
}
恢復一條已刪去的邊:
void resuedge(int No)
{
    ed[ed[No].pre].next 
= ed[ed[No].next].pre = No;
}
需要刪點的情況類似,對單點表頭處理即可。

【具體題目】PKU1815
這題就是求有向圖的字典序最小的最小點割集問題,方法是先求出最小點連通度(有關最小點連通度的求法見圖的連通度問題的求法),然后按編號遞增順序枚舉每個點,若刪去該點(其實是刪去建成的新圖中該點i'到該點附加點i''之間的邊)后圖的最小點連通度減小,則應刪去該點,否則不應刪去該點。刪去后,繼續枚舉下一個點,直到求出點割集為止。
注意,本題只有刪邊,沒有刪點,因此總表頭可以不需要,直接從ed[0]開始作單點表頭。此時,關于是否空位就剛好反過來了:如果N是奇數就要空位,N是偶數不空位(不過這題里由于建出的網絡流圖中有2*N0個結點,總是偶數,可以不管,不過本沙茶還是管這個了)。

代碼(神犇不要鄙視):
#include <iostream>
#include 
<stdio.h>
using namespace std;
#define re(i, n) for (int i=0; i<n; i++)
const int MAXN = 501, MAXM = 100000, INF = ~0U >> 2;
struct edge {
    
int a, b, fs, f, pre, next;
} ed[MAXM 
+ MAXM];
int n0, n, m = 0, s, t, start[MAXN], pr[MAXN], hs[MAXN], lev[MAXN], q[MAXN], now, flow, reslen = 0, res[MAXN];
bool vst[MAXN];
void init_d()
{
    re(i, n) {ed[i].a 
= i; ed[i].b = -1; ed[i].f = 0; ed[i].pre = ed[i].next = i;}
    
if (n % 2) m = n + 1else m = n;
}
void add_edge(int a, int b, int f)
{
    ed[m].a 
= a; ed[m].b = b; ed[m].fs = ed[m].f = f; ed[m].pre = ed[a].pre; ed[m].next = a; ed[a].pre = m; ed[ed[m].pre].next = m++;
    ed[m].a 
= b; ed[m].b = a; ed[m].fs = ed[m].f = 0; ed[m].pre = ed[b].pre; ed[m].next = b; ed[b].pre = m; ed[ed[m].pre].next = m++;
}
void init()
{
    
int x;
    scanf(
"%d%d%d"&n0, &s, &t);
    n 
= n0 << 1; s--; t += n0 - 1; init_d();
    re(i, n0) re(j, n0) {
        scanf(
"%d"&x);
        
if (i == j) {
            
if (i == s || i == t - n0) add_edge(i, i + n0, INF); else add_edge(i, i + n0, 1);
        } 
else if (x) add_edge(i + n0, j, INF);
    }
}
void aug()
{
    
int z = hs[t], i = t, p; flow += z;
    
while (i != s) {
        hs[i] 
-= z; p = pr[i]; ed[p].f -= z; ed[p ^ 1].f += z; i = ed[p].a;
        
if (!ed[p].f) now = i;
    }
}
bool dfs()
{
    re(i, n) vst[i] 
= 0; vst[s] = 1; q[0= s; lev[s] = 0;
    
int i, j, f0;
    
for (int front=0, rear=0; front<=rear; front++) {
        i 
= q[front];
        
for (int p=ed[i].next; p != i; p=ed[p].next) if (ed[p].f) {
            j 
= ed[p].b;
            
if (!vst[j]) {vst[j] = 1; q[++rear] = j; lev[j] = lev[i] + 1;}
        }
    }
    
if (!vst[t]) return 0;
    now 
= s; re(i, n) {start[i] = ed[i].next; vst[i] = 0;} hs[now] = INF;
    
bool fd;
    
while (!vst[s]) {
        
if (now == t) aug();
        fd 
= 0;
        
for (int p=start[now]; p != now; p=ed[p].next) {
            j 
= ed[p].b; f0 = ed[p].f;
            
if (lev[now] + 1 == lev[j] && !vst[j] && f0) {
                fd 
= 1; start[now] = pr[j] = p; hs[j] = hs[now] <= f0 ? hs[now] : f0; now = j; break;
            }
        }
        
if (!fd) {
            vst[now] 
= 1;
            
if (now != s) now = ed[pr[now]].a;
        }
    }
    
return 1;
}
void deledge(int No)
{
    ed[ed[No].pre].next 
= ed[No].next; ed[ed[No].next].pre = ed[No].pre;
}
void resuedge(int No)
{
    ed[ed[No].pre].next 
= ed[ed[No].next].pre = No;
}
void resu_all()
{
    re(i, n) 
for (int p=ed[i].next; p != i; p=ed[p].next) ed[p].f = ed[p].fs;
}
void solve()
{
    flow 
= 0while (dfs()) ; int f_ = flow;
    
if (!flow) {reslen = -1return;}
    re(i, m) 
if (ed[i].a + n0 == ed[i].b && ed[i].a != s && ed[i].b != t) {
        deledge(i); deledge(i 
^ 1); resu_all();
        flow 
= 0while (dfs()) ;
        
if (flow < f_) {res[reslen++= ed[i].a + 1; f_--;} else {resuedge(i ^ 1); resuedge(i);}
    }
}
void pri()
{
    
if (reslen == -1) puts("0"); else if (reslen) {
        printf(
"%d\n", reslen);
        re(i, reslen 
- 1) printf("%d ", res[i]); printf("%d\n", res[reslen - 1]);
    } 
else puts("NO ANSWER!");
}
int main()
{
    init();
    solve();
    pri();
    
return 0;
}

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲精品免费一二三区| 午夜视频在线观看一区二区三区 | 亚洲一区二区三区午夜| 亚洲精品一区久久久久久| 亚洲人成网站色ww在线| 亚洲国产精品一区二区尤物区 | 久久精品在线观看| 老司机一区二区| 女女同性精品视频| 亚洲理论电影网| 午夜精品福利视频| 麻豆国产精品va在线观看不卡 | 国产在线欧美| 亚洲国产精品黑人久久久| 一区二区三区四区五区视频 | 免费成人高清| 亚洲国产日韩综合一区| 在线一区亚洲| 久久成人18免费观看| 欧美成人dvd在线视频| 欧美午夜免费电影| 激情欧美丁香| 亚洲午夜av| 欧美成人性网| 午夜精品电影| 欧美日韩亚洲另类| 一区在线观看| 亚洲男人影院| 亚洲国产导航| 久久av在线看| 国产精品国产三级国产普通话蜜臀 | 91久久久精品| 久久精品国产999大香线蕉| 欧美乱在线观看| 国产综合激情| 欧美一区91| 99精品国产在热久久婷婷| 久久久久久九九九九| 亚洲韩日在线| 欧美在线观看一区二区| 欧美视频精品一区| 亚洲激情国产| 麻豆亚洲精品| 欧美一区二区三区四区高清| 欧美午夜宅男影院| 亚洲精品中文字幕女同| 蜜臀av在线播放一区二区三区| 亚洲无线观看| 欧美三级免费| 一本色道久久综合亚洲精品婷婷 | 免费不卡中文字幕视频| 国产一区二区三区观看 | 亚洲日韩欧美视频| 免费不卡在线视频| 亚洲春色另类小说| 免费观看亚洲视频大全| 久久精品五月| 在线不卡a资源高清| 久久琪琪电影院| 久久福利一区| 伊人狠狠色丁香综合尤物| 久久永久免费| 久久久国产成人精品| 在线成人免费观看| 免费久久久一本精品久久区| 久久午夜精品| 亚洲黄色在线视频| 日韩视频免费看| 欧美色欧美亚洲另类二区| 亚洲一区在线观看视频 | 久久久蜜桃精品| 久久深夜福利免费观看| 亚洲高清久久网| 亚洲人成毛片在线播放| 欧美日韩精品一区二区三区四区| 亚洲美女色禁图| 亚洲精品日产精品乱码不卡| 欧美视频免费在线观看| 欧美一区二区三区另类| 久久精品一二三区| 亚洲日韩视频| 亚洲免费在线观看视频| 曰本成人黄色| 日韩亚洲视频| 国产在线拍揄自揄视频不卡99| 欧美 日韩 国产精品免费观看| 欧美黄免费看| 欧美亚洲免费| 欧美α欧美αv大片| 亚洲欧美日韩直播| 久久视频在线看| 国产精品日韩欧美一区二区| 午夜久久美女| 看欧美日韩国产| 亚洲视频在线二区| 久久久久久久久久久久久女国产乱 | 欧美日韩情趣电影| 久久精品国产亚洲一区二区三区| 久久午夜羞羞影院免费观看| 亚洲一区视频在线观看视频| 欧美一区二区女人| 一区二区高清在线| 久久久久.com| 亚洲欧美国产日韩中文字幕| 久久久欧美一区二区| 亚洲欧美日产图| 老司机久久99久久精品播放免费| 亚洲综合日韩| 欧美α欧美αv大片| 午夜一区不卡| 欧美日韩一区二区三区免费看| 久久婷婷久久一区二区三区| 欧美亚州在线观看| 欧美二区在线| 国产自产精品| 亚洲欧美卡通另类91av| 9色国产精品| 欧美电影电视剧在线观看| 久久亚裔精品欧美| 国产毛片精品国产一区二区三区| 亚洲精品免费看| 亚洲精品日韩精品| 久久婷婷成人综合色| 久久精品91久久香蕉加勒比| 国产精品都在这里| 99伊人成综合| 一区二区免费在线观看| 免费视频最近日韩| 老司机久久99久久精品播放免费| 国产日韩精品在线| 正在播放日韩| 中日韩美女免费视频网址在线观看| 麻豆精品在线视频| 免费观看欧美在线视频的网站| 国产日韩av在线播放| 亚洲一二三区视频在线观看| 一区二区三区四区五区精品| 欧美激情久久久久久| 亚洲国产成人精品久久| 亚洲国产精品激情在线观看| 久久久精彩视频| 久久大香伊蕉在人线观看热2| 国产精品久久久99| 艳妇臀荡乳欲伦亚洲一区| 国产精品伦子伦免费视频| 中文av一区二区| 亚洲欧美欧美一区二区三区| 欧美特黄一区| 亚洲一区图片| 久久久国际精品| 韩国女主播一区二区三区| 久久综合久久综合久久| 欧美激情第二页| 欧美激情视频一区二区三区免费| 欧美国产欧美综合| 日韩亚洲欧美一区二区三区| 欧美日韩成人免费| 一区二区久久久久| 销魂美女一区二区三区视频在线| 国产精品sm| 亚洲欧美久久久| 久久手机精品视频| 亚洲娇小video精品| 欧美日韩一区二区三区视频 | 久久精品国产99精品国产亚洲性色| 久久精品二区三区| 亚洲第一福利在线观看| 欧美77777| 亚洲天堂av在线免费| 久久综合影音| 正在播放亚洲一区| 国产日韩在线播放| 欧美高清不卡| 亚洲视频精选| 欧美jizzhd精品欧美巨大免费| 日韩视频在线你懂得| 国产精品―色哟哟| 久久综合图片| 亚洲午夜激情免费视频| 免费观看日韩| 欧美在线观看视频一区二区三区 | 欧美国产日韩二区| 亚洲在线观看视频网站| 久久综合综合久久综合| 一区二区av| 在线日韩中文字幕| 国产精品久久久久999| 久久五月婷婷丁香社区| 一区二区三区欧美在线| 欧美成人一区二区三区| 先锋影音久久久| 一二三区精品福利视频| 在线观看视频免费一区二区三区| 欧美日韩一区二区免费视频| 久久午夜精品| 亚洲欧美日韩在线| 亚洲一级在线观看| 亚洲人成免费| 欧美肥婆在线| 巨乳诱惑日韩免费av|