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

            Drolca

            Apologize To Drolca
            隨筆 - 28, 文章 - 1, 評論 - 6, 引用 - 0
            數(shù)據(jù)加載中……

            有一個悲劇...最小點(diǎn)割

             

            #include <iostream>
            using namespace std;
            const int maxn=200;
            const int INF=1000000;
            int g[maxn][maxn];
            int f[maxn][maxn];
            int r[maxn][maxn];

            int Edmonds_Karp(int n,int g[][maxn],int s,int t,int f[][maxn]) 

                
            int i,j,k,c,head,tail,flow=0 ; 
                
            int prev[maxn],visit[maxn],q[maxn]; 
                
            for(i=0;i<n;i++)for(j=0;j<n;j++
                

                    f[i][j]
            =0 ; 
                    r[i][j]
            =g[i][j]; 
                }
             
                
            while(1
                

                    head
            =tail=0 ; 
                    memset(visit,
            0,sizeof(visit)); 
                    q[tail
            ++]=s ; 
                    prev[s]
            =-1 ; 
                    visit[s]
            =1 ; 
                    
            while(head<tail) 
                    

                        k
            =q[head++]; 
                        
            for(i=0;i<n;i++)
                            
            if(!visit[i]&&r[k][i]>0
                            

                                visit[i]
            =1 ; 
                                prev[i]
            =k ; 
                                
            if(i==t)goto next ; 
                                q[tail
            ++]=i ; 
                            }
             
                    }
             
                    next : 
                    
            if(!visit[t])break ; 
                    
            for(c=INT_MAX,j=t;j!=s;j=i) 
                    

                        i
            =prev[j]; 
                        
            if(c>r[i][j])c=r[i][j]; 
                    }
             
                    
            for(j=t;j!=s;j=i) 
                    

                        i
            =prev[j]; 
                        f[i][j]
            +=c ; 
                        f[j][i]
            =-f[i][j]; 
                        r[i][j]
            =g[i][j]-f[i][j]; 
                        r[j][i]
            =g[j][i]-f[j][i]; 
                    }
             
                    flow
            +=c ; 
                }
             
                
            return flow ; 
            }


            void Floyd(int n,int map[][maxn])
            {
                
            int i,j,k;
                
            for(k=0;k<n;k++){
                    
            for(i=0;i<n;i++){
                        
            for(j=0;j<n;j++){
                            
            if(map[i][k]+map[k][j]<map[i][j])
                                map[i][j]
            =map[i][k]+map[k][j];
                        }

                    }

                }

            }


            int mat[maxn][maxn];
            int tmp[maxn][maxn];
            int mp[maxn][maxn];
            int make[maxn][maxn];
            int bus[maxn];

            void build(int id)
            {
                
            int i,j;
                
            for(i=0;i<id;i++){
                    make[i][i
            +id]=1;
                    
            for(j=0;j<id;j++){
                        
            if(mp[i][j])
                            make[i
            +id][j]=INF;
                    }

                }

            }


            int main()
            {
                
            int n,m,p;
                
            while(scanf("%d%d%d",&n,&m,&p)!=EOF)
                
            {
                    
            if(n==0&&m==0&&p==0)break;
                    
            int i,j;
                    
            for(i=0;i<n;i++){
                        
            for(j=0;j<n;j++){
                            mat[i][j]
            =INF;
                        }

                        mat[i][i]
            =0;
                    }

                    
            for(i=0;i<m;i++){
                        
            int u,v;
                        scanf(
            "%d%d",&u,&v);
                        
            --u;--v;
                        mat[u][v]
            =1;
                    }

                
                    
            for(i=0;i<n;i++){
                        
            for(j=0;j<n;j++){
                        tmp[i][j]
            =mat[i][j];
                        }

                    }
                
                    Floyd(n,tmp);

                    
            if(p<tmp[0][n-1])
                    
            {
                        printf(
            "0\n");
                        
            continue;
                    }

                    
            int id=0;
                    
            for(i=0;i<n;i++){
                        
            if(tmp[0][i]+tmp[i][n-1]<=p){
                            bus[id
            ++]=i;
                        }

                    }

                    
            for(i=0;i<id;i++){
                        
            for(j=0;j<id;j++){
                            mp[i][j]
            =0;
                            
            if(mat[bus[i]][bus[j]]==1)
                                mp[i][j]
            =1;
                        }

                    }

                    build(id);
                    printf(
            "%d\n",Edmonds_Karp(2*id,make,id,id-1,f));
                }

                
            return 0;
            }

            posted on 2009-11-21 21:08 Drolca 閱讀(241) 評論(0)  編輯 收藏 引用


            只有注冊用戶登錄后才能發(fā)表評論。
            網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


            avtt天堂网久久精品| 国产真实乱对白精彩久久| 久久er国产精品免费观看8| 色成年激情久久综合| 午夜精品久久影院蜜桃| 久久精品国产亚洲AV久| 69久久精品无码一区二区| 狠狠色综合久久久久尤物| 久久天天躁狠狠躁夜夜2020一 | 日本精品久久久中文字幕| 国产精品久久久天天影视香蕉 | 美女写真久久影院| 久久久久久av无码免费看大片| 思思久久精品在热线热| 99久久精品免费| 日产精品久久久一区二区| 久久无码AV中文出轨人妻| 精品久久久久久无码专区不卡 | 麻豆国内精品久久久久久| 日本强好片久久久久久AAA| 久久人人爽人人澡人人高潮AV | 精品午夜久久福利大片| AV无码久久久久不卡蜜桃| 久久精品无码免费不卡| 97精品国产97久久久久久免费| 亚洲精品美女久久777777| 久久精品国产欧美日韩99热| 国产精品一区二区久久精品无码 | 久久国产精品偷99| 欧美伊香蕉久久综合类网站| 国内精品久久久久久99| 久久亚洲国产成人精品性色| 久久精品成人欧美大片| 久久午夜无码鲁丝片秋霞 | 欧美久久久久久午夜精品| 亚洲欧美精品伊人久久| 日韩亚洲欧美久久久www综合网| 久久精品国产99国产电影网| 国产一区二区三区久久精品| 99精品伊人久久久大香线蕉| 精品国产一区二区三区久久蜜臀|