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

考慮這樣一種網絡流問題:需要對同一個圖求多次最大流。則在每次求最大流之前,需要將所有邊的容量全部恢復到初始值(求最大流的過程中,邊的容量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>
            久久亚洲综合色| 亚洲免费一在线| 亚洲男女自偷自拍| 久久精品国产亚洲aⅴ| 亚洲第一在线综合在线| 国产精品视频一二| 欧美va亚洲va香蕉在线| 亚洲三级视频在线观看| 性18欧美另类| 欧美黄免费看| 亚洲国产老妈| 黄色在线一区| 欧美人与性动交cc0o| 久久男人资源视频| 激情综合在线| 性欧美大战久久久久久久免费观看| 国产精品狼人久久影院观看方式| 国产精品日韩在线播放| 午夜精品久久一牛影视| 亚洲国产专区| 日韩天堂在线视频| 亚洲精品一级| 亚洲欧美一区二区三区在线| 亚洲高清毛片| 亚洲国产欧美日韩精品| 日韩亚洲精品视频| 美女黄网久久| 久久久久国产成人精品亚洲午夜| 国产日韩欧美麻豆| 亚洲第一成人在线| 亚洲天堂av高清| 亚洲一二三区在线| 亚洲欧美日韩在线一区| 亚洲电影免费在线| 亚洲人体影院| 揄拍成人国产精品视频| 亚洲欧美视频一区| 日韩天堂在线视频| 欧美日韩免费在线视频| 欧美一级久久久久久久大片| 亚洲视频综合| 亚洲高清久久久| 亚洲精品乱码久久久久久蜜桃麻豆 | 这里只有精品电影| 欧美1区2区视频| 91久久精品一区二区三区| 美日韩精品视频| 亚洲精品一级| 午夜精品久久久久99热蜜桃导演| 久久久精品一品道一区| 欧美黄色一区| 午夜久久99| 久久精品久久99精品久久| 国产女精品视频网站免费| 免费欧美日韩| 亚洲视频欧美在线| 亚洲自拍电影| 国产亚洲精品激情久久| 久久亚洲国产精品一区二区| 欧美黑人在线观看| 午夜日韩在线观看| 亚洲美女中出| 一区国产精品| 亚洲美女中文字幕| 在线观看国产成人av片| 亚洲国产精品尤物yw在线观看| 亚洲欧洲日本mm| 99国产一区| 国产亚洲精品久久久| 狠狠色综合一区二区| 亚洲男人的天堂在线aⅴ视频| 日韩视频在线观看一区二区| 国产精品久久久久久超碰| 一本色道久久88综合亚洲精品ⅰ | 国产亚洲福利社区一区| 欧美一区二区精品在线| 欧美日本在线播放| 黄色成人在线| 亚洲免费在线视频| 亚洲区中文字幕| 亚洲视频一二| 一区二区久久久久| 久久精品亚洲一区二区| 久久久99免费视频| 欧美午夜精品久久久久久浪潮| 国产午夜精品久久久久久久| 日韩亚洲精品在线| 亚洲二区免费| 黄色成人免费观看| 久久一二三四| 国产精品萝li| 蜜桃av一区| 亚洲欧美国产日韩天堂区| 国产亚洲欧美日韩美女| 亚洲小说春色综合另类电影| 亚洲精选一区二区| 欧美美女bbbb| 欧美激情91| 国产精品久久久久久久久久免费| 亚洲人体一区| 久久九九久精品国产免费直播| 久久久久久网址| 国产伦理一区| 欧美一区二区三区久久精品| 欧美高清在线视频观看不卡| 久久久国产一区二区| 久久九九国产精品| 久久精品一区四区| 欧美成人一品| 亚洲精品中文字幕女同| 久久久视频精品| 亚洲第一福利视频| 亚洲经典在线看| 美女视频一区免费观看| 亚洲精品一区在线| 夜夜嗨av一区二区三区网站四季av| 亚洲综合久久久久| 国产精品九九| 欧美一区中文字幕| 亚洲国产精品成人va在线观看| 久久成人免费电影| 亚洲欧美日韩国产综合在线| 国产亚洲视频在线| 久久精品国产免费观看| 亚洲精品1234| 亚洲一区二区三区在线观看视频 | 亚洲电影网站| 裸体一区二区三区| 一区二区三区高清不卡| 久久久99精品免费观看不卡| 欧美中文字幕在线| 精品成人一区二区| 久久av老司机精品网站导航| 欧美激情亚洲一区| 欧美日韩国产在线播放| 欧美www在线| 亚洲午夜精品一区二区| 国产精品剧情在线亚洲| 美日韩在线观看| 99riav1国产精品视频| 午夜影院日韩| 日韩视频一区二区| 国产精品性做久久久久久| 免费在线一区二区| 亚洲视频电影图片偷拍一区| 亚洲一区二区三区视频播放| 国产精品99久久久久久久久久久久 | 久久综合狠狠综合久久激情| 日韩视频在线一区| 久久久999精品免费| 国产真实乱子伦精品视频| 欧美成人资源网| 亚洲欧美日韩在线播放| 91久久午夜| 欧美在线视频一区二区三区| 影音先锋亚洲精品| 欧美中文在线字幕| 一本色道久久综合狠狠躁篇的优点 | 亚洲人成小说网站色在线| 国产精品亚洲片夜色在线| 你懂的视频一区二区| 一区二区三区欧美成人| 亚洲一区亚洲二区| 日韩一级黄色大片| 国产一区二区日韩精品| 欧美午夜欧美| 噜噜噜噜噜久久久久久91| 亚洲视频视频在线| 99re8这里有精品热视频免费| 欧美影院在线播放| 午夜久久99| 久久激情综合网| 亚洲一区二区三区乱码aⅴ蜜桃女| 亚洲国产精品va在看黑人| 国产精品视频导航| 国产精品久久毛片a| 欧美激情精品久久久久久免费印度| 久久夜色精品| 欧美在线观看视频一区二区| 欧美一区二区观看视频| 亚洲网站在线| 免费精品视频| 亚洲成人资源网| 国产一区观看| 国产日本欧美视频| 欧美视频在线免费看| 欧美日韩在线影院| 欧美激情精品久久久久久黑人| 欧美二区在线看| 蜜臀久久99精品久久久久久9| 亚洲视频在线观看三级| 亚洲精品永久免费| 欧美综合二区| 亚洲欧洲日本一区二区三区| 99精品国产热久久91蜜凸| 欧美专区中文字幕| 欧美一区二区视频网站| 噜噜噜91成人网| 开心色5月久久精品| 亚洲国产成人精品女人久久久 |