• <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>
            獨立博客: 哲學與程序

            哲學與程序

            判斷最小割是否唯一 ZOJ@2587

            ZOJ@2587
            題意:判斷最小割是否唯一。
            思路(轉):判斷最小割是否唯一,先求一次最大流,然后在殘留網絡中分別從源匯開始dfs一次,找出最小割[S,T],如果SUT不包含所有點,那么最小割不唯一。假設點i不被SUT包含,那么殘留網絡中s不能到達i,i不能到達t,即進入i的邊和從i出去的邊都滿流,假設某條進入i的邊x滿流,這些流量從若干條邊y流出i,那么,如果選x為割邊,或者選所有對應的y為割邊,不會影響最大流,即最小割容量不變,最小割也就不唯一。
            // 2390377      2011-01-21 16:52:22        Accepted      2587      C++      150      3324      redsea
            //  Dinic最大流 O(V^2 * E) 
            #include<cstdio>
            #include
            <algorithm>
            #include
            <cstring>
            using namespace std;
            #define N 802
            #define E 200050
            #define typec int                   // type of cost 
            const typec inf = 30000000;       // max of cost 
            struct edge { 
                
            int x, y, nxt; typec c;
            }bf[E]; 
            int ne, head[N], cur[N], ps[N], dep[N]; 

            void addedge(int x, int y, typec c) 
            {  
                
            // add an arc(x -> y, c); vertex: 0 ~ n-1; 
                  bf[ne].x = x; bf[ne].y = y; bf[ne].c = c; 
                  bf[ne].nxt 
            = head[x]; head[x] = ne++
                  bf[ne].x 
            = y; bf[ne].y = x; bf[ne].c = 0
                  bf[ne].nxt 
            = head[y]; head[y] = ne++

            typec flow(
            int n, int s, int t) 

                  typec tr, res 
            = 0;
                
            int i, j, k, f, r, top; 
                
            while (1) { 
                    memset(dep, 
            -1, n * sizeof(int)); 
                       
            for (f = dep[ps[0= s] = 0, r = 1; f != r; ) 
                        
            for (i = ps[f++], j = head[i]; j; j = bf[j].nxt) 
                        { 
                             
            if (bf[j].c && -1 == dep[k = bf[j].y]){ 
                                  dep[k] 
            = dep[i] + 1; ps[r++= k; 
                                  
            if (k == t) { f = r; break; } 
                            } 
                           } 
                    
            if (-1 == dep[t]) break
                    memcpy(cur, head, n 
            * sizeof(int)); 
                    
            for (i = s, top = 0; ; ) { 
                        
            if (i == t) { 
                            
            for (k = 0, tr = inf; k < top; ++k) 
                                 
            if (bf[ps[k]].c < tr) 
                                    tr 
            = bf[ps[f = k]].c; 
                            
            for (k = 0; k < top; ++k) 
                                bf[ps[k]].c 
            -= tr, bf[ps[k]^1].c += tr; 
                               res 
            += tr;  i = bf[ps[top = f]].x; 
                        } 
                        
            for  (j=cur[i]; cur[i]; j = cur[i] = bf[cur[i]].nxt) 
                               
            if  (bf[j].c && dep[i]+1 == dep[bf[j].y])break
                        
            if (cur[i]) { 
                               ps[top
            ++= cur[i]; 
                               i 
            = bf[cur[i]].y; 
                           } 
                           
            else { 
                             
            if (0 == top) break
                               dep[i] 
            = -1; i = bf[ps[--top]].x; 
                           } 
                       } 
                } 
                
            return res; 
            }
            int cnts, cntt;
            int flag1[N], flag2[N];

            void dfs1(int v)
            {
                flag1[v] 
            = 1;
                cnts
            ++;
                
            for(int i = head[v]; i != 0; i = bf[i].nxt){
                    
            if(flag1[bf[i].y]==0 && bf[i].c)
                        dfs1(bf[i].y);
                }
            }
            void dfs2(int v)
            {
                flag2[v] 
            = 1;
                cntt
            ++;
                
            for(int i = head[v]; i != 0; i = bf[i].nxt){
                    
            if(flag2[bf[i].y]==0 && bf[i^1].c)
                        dfs2(bf[i].y);
                }
            }
            int main()
            {
                
            int n, m, s, t, x, y, c;
                
            while(scanf("%d%d%d%d",&n,&m,&s,&t), n||m||s||t)
                {
                    memset(head,
            0,sizeof(head));
                    ne 
            = 2;
                    
            for(int i = 0; i < m; i++){
                        scanf(
            "%d%d%d",&x,&y,&c);
                        x
            --,y--;
                        addedge(x,y,c);
                        addedge(y,x,c);
                    }
                    s
            --;t--;
                    flow(n,s,t);
                    memset(flag1,
            0,sizeof(flag1));
                    memset(flag2,
            0,sizeof(flag2));
                    cnts 
            = cntt = 0;
                    dfs1(s);
                    dfs2(t);
                    
            if(cnts+cntt!=n)
                        printf(
            "AMBIGUOUS\n");
                    
            else printf("UNIQUE\n");
                }
                
            return 0;
            }


            posted on 2011-01-21 17:11 哲學與程序 閱讀(805) 評論(0)  編輯 收藏 引用 所屬分類: Algorithm

            導航

            公告

            歡迎訪問 http://zhexue.sinaapp.com

            常用鏈接

            隨筆分類(37)

            隨筆檔案(41)

            Algorithm

            最新隨筆

            搜索

            最新評論

            獨立博客: 哲學與程序
            99久久国产综合精品网成人影院| 亚洲国产天堂久久综合| 久久久老熟女一区二区三区| 久久久av波多野一区二区| 久久国产高潮流白浆免费观看| 99久久久精品免费观看国产| 国内精品久久国产大陆| 婷婷综合久久中文字幕| 久久久黄色大片| 久久99精品久久久久子伦| 久久狠狠一本精品综合网| 欧美日韩精品久久久免费观看| 2020久久精品国产免费| 久久天天躁狠狠躁夜夜不卡| 囯产极品美女高潮无套久久久| 久久狠狠色狠狠色综合| 性做久久久久久久久老女人| 久久超乳爆乳中文字幕| 大香伊人久久精品一区二区| www久久久天天com| 久久精品国产乱子伦| 精品乱码久久久久久夜夜嗨| 久久99国产精品尤物| 欧美亚洲国产精品久久| 久久亚洲av无码精品浪潮| 久久国产精品-国产精品| 亚洲国产精品高清久久久| 久久久久18| 国内精品伊人久久久久网站| 97久久天天综合色天天综合色hd| 国产精品久久久久久五月尺| 亚洲国产精品婷婷久久| 99久久精品午夜一区二区| 色综合久久久久综合体桃花网| 国产精品久久新婚兰兰| 久久久久久免费视频| 久久乐国产精品亚洲综合| 国产精品美女久久久免费| 88久久精品无码一区二区毛片 | 99精品国产99久久久久久97| 久久青青草原精品国产软件|