• <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 哲學與程序 閱讀(802) 評論(0)  編輯 收藏 引用 所屬分類: Algorithm

            導航

            公告

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

            常用鏈接

            隨筆分類(37)

            隨筆檔案(41)

            Algorithm

            最新隨筆

            搜索

            最新評論

            獨立博客: 哲學與程序
            少妇无套内谢久久久久| 久久综合丁香激情久久| 久久青青草原亚洲av无码| 久久久久一级精品亚洲国产成人综合AV区| 九九久久精品国产| 久久精品桃花综合| 国内精品久久久久影院优 | 国产美女久久久| 国产精品99久久久久久董美香 | 欧美成a人片免费看久久| 色播久久人人爽人人爽人人片aV| 久久久久久精品久久久久| 久久久91精品国产一区二区三区| 久久精品亚洲福利| 久久国产欧美日韩精品| 久久夜色精品国产亚洲av| 国内精品伊人久久久久av一坑| 一个色综合久久| 久久se精品一区精品二区国产| 伊人久久大香线蕉av一区| 久久99精品免费一区二区| 国产91色综合久久免费| 久久综合亚洲色一区二区三区| 国产成人久久精品二区三区| 久久久久久久精品妇女99| 久久久久久久综合日本| 91久久九九无码成人网站| 99久久精品国产高清一区二区| 久久婷婷色综合一区二区| 亚洲一级Av无码毛片久久精品| 狠狠人妻久久久久久综合| av无码久久久久不卡免费网站| 久久久久亚洲av综合波多野结衣| 欧美国产精品久久高清| 久久青青草原亚洲av无码| 久久久综合香蕉尹人综合网| 国产精品免费久久久久影院| 国产成人精品久久亚洲| 精品多毛少妇人妻AV免费久久| 国产日韩欧美久久| 欧美日韩中文字幕久久久不卡|