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

            最新隨筆

            搜索

            最新評論

            獨立博客: 哲學與程序
            国内精品久久久久久久亚洲 | 9久久9久久精品| 狠狠色综合网站久久久久久久高清| 久久男人中文字幕资源站| 人妻中文久久久久| 久久精品中文字幕一区| 国产精品99久久99久久久| 国内精品欧美久久精品| 久久这里都是精品| MM131亚洲国产美女久久| 久久国产综合精品五月天| 亚洲精品国产综合久久一线| 久久久无码精品亚洲日韩蜜臀浪潮| 国产成人综合久久综合| 亚洲精品NV久久久久久久久久| 久久久久久久久无码精品亚洲日韩| 国产香蕉97碰碰久久人人| 午夜精品久久久久久久| 久久强奷乱码老熟女网站| 久久婷婷五月综合97色一本一本 | 国产精品久久久久久久午夜片 | 国产三级久久久精品麻豆三级| 九九久久精品无码专区| 久久久久女人精品毛片| 亚洲中文字幕伊人久久无码 | 中文字幕一区二区三区久久网站| 国内精品伊人久久久影院 | 狠狠精品久久久无码中文字幕| 97精品国产97久久久久久免费| 久久久久久毛片免费播放| 亚洲人成网站999久久久综合 | 99精品国产综合久久久久五月天| 青青青国产精品国产精品久久久久 | 国产精品美女久久久久av爽| 国产麻豆精品久久一二三| 一本一本久久A久久综合精品 | 久久久久一本毛久久久| 久久国产精品-国产精品| 国产亚洲欧美成人久久片 | 精品熟女少妇a∨免费久久| 亚洲国产精品无码久久一区二区|