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

            misschuer

            常用鏈接

            統計

            積分與排名

            百事通

            最新評論

            poj 3189 Steady Cow Assignment

            //題意 貌似是求滿意度的范圍最小.
            //順便貼個以前找到的, 效率還算可以的最大流模板

            #include 
            <cstdio>
            #include 
            <iostream>
            #include 
            <cstring>
            using namespace std;
            const int MAXV = 1210;
            const int MAXE = 1000010;
            const int inf = 1 << 29;

            int head[MAXV], N;

            struct Edge {
                
                
            int v, next, w;
            }edge[MAXE];

            int n, m, s, t, cnt, ans;

            void add(int u, int v, int w) {
                
                edge[cnt].v 
            = v;
                edge[cnt].w 
            = w;
                edge[cnt].next 
            = head[ u ];
                head[ u ] 
            = cnt++;
                
                edge[cnt].v 
            = u;
                edge[cnt].w 
            = 0;
                edge[cnt].next 
            = head[ v ];
                head[ v ] 
            = cnt++;
            }

            int sap() {
                
                
            int pre[MAXV], cur[MAXV], dis[MAXV], gap[MAXV];
                
            int flow = 0, aug = inf, u;
                
            bool flag;
                
                
            for(int i = 0; i <= N; ++ i) {
                    
                    cur[ i ] 
            = head[ i ];
                    gap[ i ] 
            = dis[ i ] = 0;
                }
                
                gap[ s ] 
            = N;
                u 
            = pre[ s ] = s;
                
                
            while(dis[ s ] < N) {
                    
                    flag 
            = 0;
                    
                    
            for(int &= cur[ u ]; j != -1; j = edge[ j ].next) {
                        
                        
            int v = edge[ j ].v;
                        
                        
            if(edge[ j ].w > 0 && dis[ u ] == dis[ v ] + 1) {
                            
                            flag 
            = 1;
                            
                            
            if(edge[ j ].w < aug) aug = edge[ j ].w;
                            
                            pre[ v ] 
            = u;
                            u 
            = v;
                            
                            
            if(u == t) {
                                
                                flow 
            += aug;
                                
                                
            while(u != s) {
                                    
                                    u 
            = pre[ u ];
                                    edge[cur[ u ]  ].w 
            -= aug;
                                    edge[cur[ u ]
            ^1].w += aug;
                                }
                                aug 
            = inf;
                            }
                            
            break;
                        }
                    }
                    
                    
            if(flag) continue;
                    
                    
            int mindis = N;
                    
                    
            for(int k = head[ u ]; k != -1; k = edge[ k ].next) {
                        
                        
            int v = edge[ k ].v;
                        
                        
            if(edge[ k ].w > 0 && dis[ v ] < mindis) {
                            
                            mindis 
            = dis[ v ];
                            cur[ u ] 
            = k;
                        }
                    }
                    
                    
            if( (--gap[dis[ u ]]) == 0break;
                    
                    dis[ u ] 
            = mindis + 1;
                    gap[dis[ u ]] 
            ++;
                    u 
            = pre[ u ];
                }
                
            return flow;
            }

            int cap[30], i, j;
            int mat[MAXV][30];

            int main() {
                
                
            while(scanf("%d %d"&m, &n) != EOF) {
                    
                    
            for(i = 1; i <= m; ++ i) {
                        
                        
            for(j = 1; j <= n; ++ j) {
                            
                            scanf(
            "%d"&mat[ i ][ j ]);
                        }
                    }
                    
                    
            for(i = 1; i <= n; ++ i) {
                        
                        scanf(
            "%d"&cap[ i ]);
                    }
                    
                    
            int left = 1, right = 1int mmax = n;
                    
                    
            while(left <= n && right <= n) {
                        
                        
                        s 
            = 0; t = m + n + 1; cnt = 0; N = t + 1;
                        memset(head,
            -1,sizeof(head));
                        
                        
            for(i = 1; i <= m; ++ i) {
                            
                            add(s, i, 
            1);
                        }
                        
                        
            for(i = 1; i <= m; ++ i) {
                            
                            
            for(j = 1; j <= n; ++ j) {
                                
                                
            if(left <= j && j <= right) {
                                    
                                    add(i, mat[ i ][ j ]
            +m, 1);
                                }
                            }
                        }
                        
                        
            for(j = 1; j <= n; ++ j) {
                            
                            add(m
            +j, t, cap[ j ]);
                        }
                        
                        ans 
            = sap();
                        
                        
            if(ans == m) {
                            
                            
            if(mmax > right - left + 1) mmax = right - left + 1;
                            left 
            ++;
                        }
                        
            else right ++;
                    }
                    
                    printf(
            "%d\n", mmax);
                }
                
            return 0;
            }

            posted on 2011-03-24 11:05 此最相思 閱讀(148) 評論(0)  編輯 收藏 引用

            国产午夜电影久久| 久久狠狠高潮亚洲精品| 欧美成a人片免费看久久| 久久天天躁狠狠躁夜夜av浪潮 | 伊人久久综合无码成人网| 国产成人精品综合久久久| 四虎国产永久免费久久| 久久精品国产男包| 久久影院亚洲一区| 国内精品久久久久| 狠狠色噜噜色狠狠狠综合久久 | 国内精品九九久久久精品| 久久国产精品免费一区二区三区| 久久乐国产综合亚洲精品| 精品久久久久久中文字幕| 久久久久久久精品妇女99| 九九久久精品国产| 99久久精品国内| 三上悠亚久久精品| 亚洲国产视频久久| 国产农村妇女毛片精品久久| 欧洲人妻丰满av无码久久不卡| 一本久道久久综合狠狠躁AV| 91超碰碰碰碰久久久久久综合| 熟妇人妻久久中文字幕| 久久这里都是精品| 欧美日韩精品久久久久| 日产久久强奸免费的看| 狠狠人妻久久久久久综合蜜桃| 国产精品天天影视久久综合网| 久久精品夜夜夜夜夜久久| 午夜天堂精品久久久久| 亚洲av成人无码久久精品| 囯产精品久久久久久久久蜜桃| 一97日本道伊人久久综合影院| 99久久免费国产精品特黄| 亚洲美日韩Av中文字幕无码久久久妻妇 | 久久99精品久久久久久噜噜| 久久久久一本毛久久久| 日批日出水久久亚洲精品tv| 久久SE精品一区二区|