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

                最大流最小割定理介紹:把一個流網絡的頂點集劃分成兩個集合S和T,使得源點s ∈S且匯點t ∈T,割(S,T)的容量C(S,T) = ∑ Cuv, 其中u∈S且v∈T。
            從直觀上看,截集(S,T)是從源點s到匯點t的必經之路,如果該路堵塞則流從s無法到達t。于是我們可以得到最大流最小割定理:任意一個流網絡的最大流量等于該網絡的最小的割的容量。
                最小割問題:ZOJ@2788
                題意:在一棟樓中有若干房間,一個房間有若干扇門通往別的房間,現在房間中有一個secure room。由于有一些敵人已經潛入某一些房間,現在問關閉最少的一些門,使得敵人們到達不了secure room。對于一扇門連通A-B,若控制面板CP在A方,則總是可以從A->B,但是若關上門,則B不能到達A。
                解法:最大流。首先通過floyd算法判斷是否敵人總是可以到達secure room,如果能直接輸出結果;如果不能,則通過最大流算法求解即可。對于一扇門A-B,CP在A方則從A連一條有向邊到B,容量為無窮大,從B連一條邊到A,容量為1。兩個房間之間可能有多扇門,邊的容量累加即可。增加一個源點,從源點連邊到有敵人的房間,容量為1。對圖求最大流即為結果。
            代碼如下:
            // 2389732      2011-01-20 15:06:16        Accepted      2788      C++      0      184      *******
            #include<cstdio>
            #include
            <algorithm>
            #include
            <cstring>
            #define MAXN 22
            #define inf 1000000
            int map[MAXN][MAXN];
            int flow[MAXN][MAXN];
            int max_flow(int n,int mat[][MAXN],int source,int sink,int flow[][MAXN]){ 
                
            int pre[MAXN],que[MAXN],p,q,t,i,j;
                
            int d[MAXN];
                
            int tmp;
                
            if (source==sink) return inf; 
                
            for (i=0;i<n;i++
                    
            for    (j=0;j<n;flow[i][j++]=0); 
                
            for (;;){ 
                    
            for (i=0;i<n;pre[i++]=0); 
                    pre[t
            =source]=source+1,d[t]=inf; 
                    
            for (p=q=0;p<=q&&!pre[sink];t=que[p++]) 
                           
            for (i=0;i<n;i++
                            
            if (!pre[i]&&(tmp=mat[t][i]-flow[t][i])) 
                                 pre[que[q
            ++]=i]=t+1,d[i]=d[t]<tmp?d[t]:tmp; 
                            
            else if (!pre[i]&&(tmp=flow[i][t])) 
                             pre[que[q
            ++]=i]=-t-1,d[i]=d[t]<tmp?d[t]:tmp; 
                    
            if (!pre[sink]) break
                    
            for (i=sink;i!=source;) 
                           
            if (pre[i]>0
                            flow[pre[i]
            -1][i]+=d[sink],i=pre[i]-1
                           
            else 
                            flow[i][
            -pre[i]-1]-=d[sink],i=-pre[i]-1
                } 
                tmp 
            = 0;
                
            for (i=0;i<n;tmp+=flow[source][i++]); 
                
            return tmp; 
            }
            int main()
            {
                
            int T, n, m, k, y;
                
            char str[3];
                
            int d[MAXN];
                scanf(
            "%d",&T);
                
            while(T--)
                {
                    scanf(
            "%d%d",&m,&n);
                    n
            ++;
                    memset(map,
            0,sizeof(map));
                    memset(flow,
            0,sizeof(flow));
                    
            for(int i = 1; i <= m; i++){
                        scanf(
            "%s %d",str, &k);
                        d[i] 
            = (strlen(str) == 1);
                        
            if(d[i])map[0][i] = inf;
                        
            for(int j = 0; j < k; j++){
                            scanf(
            "%d",&y);
                            y
            ++;
                            flow[i][y] 
            = 1;
                            map[i][y] 
            += inf;
                            map[y][i] 
            += 1;
                        }
                    }
                    
            // Checking always exists a path to secure room. Output "PANIC ROOM BREACH"!
                    for(k = 1; k <= m; k++)  
                        
            for(int i = 1; i <= m; i++)
                            
            for(int j = 1; j <= m; j++)
                                
            if(flow[i][k]&&flow[k][j])
                                    flow[i][j] 
            = 1;
                    
            for(k = 1; k <= m && !(d[k] && flow[k][n]); k++);
                    
            if(k <= m)
                        printf(
            "PANIC ROOM BREACH\n");
                    
            else
                        printf(
            "%d\n", max_flow(m+1,map,0,n,flow));
                }
                
            return 0;
            }


            posted on 2011-01-20 15:52 哲學與程序 閱讀(359) 評論(0)  編輯 收藏 引用 所屬分類: AlgorithmC & C++

            導航

            公告

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

            常用鏈接

            隨筆分類(37)

            隨筆檔案(41)

            Algorithm

            最新隨筆

            搜索

            最新評論

            獨立博客: 哲學與程序
            韩国免费A级毛片久久| 色综合久久天天综线观看| 久久天堂AV综合合色蜜桃网| 久久久久久久久久久精品尤物 | 久久久久久久久久久| 亚洲AV日韩精品久久久久久| 久久精品国产精品青草app| 香蕉aa三级久久毛片| 久久精品国产亚洲AV麻豆网站| 亚洲综合婷婷久久| 日韩人妻无码一区二区三区久久99| 亚洲精品乱码久久久久久| 99久久精品免费国产大片| 中文字幕久久波多野结衣av| 99久久精品免费国产大片| 午夜精品久久久久久久久| 久久国产乱子伦精品免费午夜| 2020久久精品亚洲热综合一本| 国产精品丝袜久久久久久不卡| 蜜臀久久99精品久久久久久小说| 精品久久久久国产免费| 久久精品国产亚洲欧美| 五月丁香综合激情六月久久| 一极黄色视频久久网站| 久久93精品国产91久久综合| 国产精品久久久久久久久免费| 精品无码久久久久国产动漫3d| 深夜久久AAAAA级毛片免费看| 好久久免费视频高清| 99999久久久久久亚洲| 精品国产一区二区三区久久久狼 | 麻豆精品久久久一区二区| 精品国产青草久久久久福利| 久久精品视频一| 久久久久亚洲AV成人网人人网站| 国产精品内射久久久久欢欢| 2020最新久久久视精品爱| 久久91亚洲人成电影网站| 亚洲国产精品成人久久| 蜜臀av性久久久久蜜臀aⅴ| 久久久久亚洲AV成人片|