• <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

            最新隨筆

            搜索

            最新評論

            獨立博客: 哲學與程序
            国内高清久久久久久| 亚洲国产精品久久| 久久久久久精品成人免费图片 | 亚洲精品乱码久久久久久蜜桃图片| 少妇久久久久久被弄到高潮| 久久久久亚洲精品日久生情| 国产精品免费看久久久| 久久精品国产亚洲7777| 久久综合色老色| 国产精品免费久久| 久久精品无码午夜福利理论片| 国产午夜精品久久久久九九电影| 少妇熟女久久综合网色欲| 久久er热视频在这里精品| 伊人久久综合精品无码AV专区 | 国产欧美一区二区久久| 欧美性大战久久久久久| 蜜臀久久99精品久久久久久小说| 久久久久九国产精品| 99久久精品费精品国产一区二区| 久久免费视频6| 香蕉久久一区二区不卡无毒影院| 久久久久久国产精品无码下载| 精品综合久久久久久88小说| 亚洲国产欧美国产综合久久 | 99久久精品免费看国产一区二区三区| 亚洲嫩草影院久久精品| 久久精品国产亚洲AV无码娇色| 婷婷国产天堂久久综合五月| 久久本道综合久久伊人| 久久精品www| 亚洲一区中文字幕久久| 国产福利电影一区二区三区,免费久久久久久久精 | 欧美大香线蕉线伊人久久| 少妇无套内谢久久久久| 久久亚洲精品无码aⅴ大香| 亚洲国产成人久久综合区| 色老头网站久久网| 久久精品人人做人人爽电影| 久久精品国产久精国产一老狼| 99久久香蕉国产线看观香|