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

考慮這樣一種網絡流問題:需要對同一個圖求多次最大流。則在每次求最大流之前,需要將所有邊的容量全部恢復到初始值(求最大流的過程中,邊的容量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>
            99精品国产在热久久| 亚洲午夜精品久久| 亚洲欧美激情在线视频| 亚洲乱码国产乱码精品精天堂 | 亚洲美女福利视频网站| 欧美一区三区二区在线观看| 久久深夜福利免费观看| 亚洲剧情一区二区| 欧美一区二区国产| 国产精品久久婷婷六月丁香| 亚洲国产成人av好男人在线观看| 久久国产婷婷国产香蕉| 欧美精品18+| 国产一区二区黄色| 欧美一区日本一区韩国一区| 国产日韩专区| 销魂美女一区二区三区视频在线| 久久精品99国产精品日本| 91久久久久久久久久久久久| 欧美va天堂在线| 国语自产精品视频在线看8查询8 | 一区二区视频欧美| 欧美一级片在线播放| ●精品国产综合乱码久久久久| 亚洲欧美国产精品桃花| 91久久精品国产91久久性色| 亚洲人成在线观看| 亚洲精品影院| 亚洲国产一区二区三区在线播| 亚洲视频一区在线| 在线日韩av永久免费观看| 欧美午夜电影在线观看| 欧美日韩国产成人| 国产麻豆9l精品三级站| 激情视频一区二区| 日韩一级大片| 久久久视频精品| 亚洲经典三级| 一区二区三区鲁丝不卡| 欧美一区二区三区在线免费观看| 久久精品综合一区| 欧美日韩精品在线| 狠狠入ady亚洲精品| 一区二区三区欧美| 久久av最新网址| 久久婷婷丁香| 欧美精品激情在线观看| 国产一区二区三区久久精品| 亚洲综合色在线| 一区二区动漫| 欧美日韩国产综合在线| 99视频超级精品| 久久久久久久精| 激情欧美国产欧美| 欧美精品一区二区三区蜜桃| 男女视频一区二区| 亚洲深夜福利| 国产小视频国产精品| 久久精品最新地址| 一区二区三区视频观看| 老鸭窝亚洲一区二区三区| 亚洲精品一区二区在线| 一区二区三区欧美激情| 亚洲一区二区高清| 91久久久久久久久| 亚洲激情视频在线| 亚洲伦伦在线| 午夜在线精品偷拍| 精品不卡在线| 亚洲欧美激情视频在线观看一区二区三区 | 午夜精彩视频在线观看不卡 | 国产一区二区三区久久精品| 国产一区二区三区成人欧美日韩在线观看 | 欧美午夜精品理论片a级按摩 | 久久夜色精品| 美女网站久久| 亚洲国产激情| 亚洲精品网站在线播放gif| 日韩午夜在线视频| 午夜日韩av| 美女精品一区| 国产精品激情| 在线看片成人| 中文精品视频一区二区在线观看| 欧美一级视频| 亚洲第一视频| 99热在这里有精品免费| 欧美在线免费观看亚洲| 免费观看在线综合| 国产精品爽黄69| 亚洲黄色毛片| 亚洲欧美在线视频观看| 免费91麻豆精品国产自产在线观看| 91久久国产综合久久蜜月精品 | 99视频超级精品| 亚洲午夜免费福利视频| 亚洲欧美日韩精品| 免费欧美电影| 亚洲影院在线| 欧美国产综合视频| 国产一区二区久久久| aa级大片欧美| 免费不卡在线观看av| 国产精品99久久久久久久久久久久 | 亚洲欧美日本国产有色| 欧美激情久久久| 精东粉嫩av免费一区二区三区| 亚洲综合欧美| 欧美激情一级片一区二区| 亚洲欧美日韩精品久久奇米色影视 | 国产亚洲精品aa| 亚洲欧美日韩一区在线观看| 亚洲福利精品| 久久精品视频在线| 国产精品专区一| 亚洲午夜在线| 国产区精品视频| 一本色道久久综合亚洲91| 蜜臀91精品一区二区三区| 午夜视频一区二区| 欧美天堂亚洲电影院在线观看 | 国产精品高潮呻吟久久av黑人| 亚洲精品影院在线观看| 欧美二区不卡| 另类av导航| 亚洲激情av| 亚洲激情社区| 欧美成人有码| 亚洲人久久久| 欧美激情亚洲另类| 噜噜噜久久亚洲精品国产品小说| 国产综合久久久久久| 久久久久国产成人精品亚洲午夜| 性色av一区二区三区| 国产视频久久| 久久久免费av| 久久久久国色av免费观看性色| 狠色狠色综合久久| 欧美刺激性大交免费视频| 久热这里只精品99re8久| 亚洲福利精品| 亚洲国产另类久久精品| 欧美华人在线视频| 欧美成人r级一区二区三区| 欧美成人高清视频| 免费成年人欧美视频| 日韩视频三区| 一本色道久久| 国产亚洲成精品久久| 免费不卡中文字幕视频| 久久综合综合久久综合| 亚洲免费高清| 一区二区国产精品| 国产一区二区三区在线观看视频 | 亚洲女性裸体视频| 久久精彩免费视频| 亚洲欧洲精品一区| 中文欧美日韩| 亚洲大胆人体视频| 一本色道综合亚洲| 国产综合久久久久影院| 亚洲黄色av一区| 国产伦精品一区二区三区| 久久免费99精品久久久久久| 欧美成人在线网站| 性感少妇一区| 欧美精品日韩精品| 久久精品理论片| 欧美精品粉嫩高潮一区二区 | 99在线|亚洲一区二区| 亚洲精品一区二区三区四区高清 | 欧美影院午夜播放| 夜夜嗨av一区二区三区| 亚洲女性裸体视频| 亚洲黄一区二区| 午夜免费在线观看精品视频| 亚洲国产婷婷综合在线精品| 在线一区观看| 亚洲精品在线视频观看| 午夜精品一区二区三区电影天堂 | 亚洲激情视频在线播放| 亚洲欧美日韩一区二区在线 | 欧美一级黄色网| 欧美喷潮久久久xxxxx| 久久一综合视频| 国产精品v欧美精品v日韩| 亚洲国产免费看| 国产伊人精品| 久久免费黄色| 久久国产主播| 亚洲小少妇裸体bbw| 欧美www在线| 美女视频黄免费的久久| 国产欧美日韩亚洲一区二区三区 | 欧美在线一二三区| 午夜精品久久久久影视| 欧美日韩免费观看一区二区三区| 久久综合图片| 精品动漫3d一区二区三区免费版| 亚洲小说欧美另类社区|