• <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 此最相思 閱讀(144) 評論(0)  編輯 收藏 引用

            中文字幕久久精品| 三级片免费观看久久| 久久人人青草97香蕉| 精品人妻伦一二三区久久| 久久免费视频6| 久久精品国产福利国产秒| 青青热久久综合网伊人| 久久影院亚洲一区| 麻豆一区二区99久久久久| 久久精品免费网站网| 91精品国产综合久久精品| 一本久久a久久精品综合香蕉| 色99久久久久高潮综合影院| 亚洲中文字幕无码久久精品1| 国产福利电影一区二区三区久久久久成人精品综合 | 66精品综合久久久久久久| 久久人妻AV中文字幕| 久久99精品久久久久久齐齐| 久久精品桃花综合| 国产69精品久久久久99| 国产精品久久久久蜜芽| 性高湖久久久久久久久AAAAA | 伊人久久综合热线大杳蕉下载| 久久精品国产亚洲AV香蕉| 国产福利电影一区二区三区,免费久久久久久久精 | 国产精品久久久久a影院| 免费观看久久精彩视频| 久久久久久精品成人免费图片 | 精品国产乱码久久久久久浪潮| 少妇内射兰兰久久| 国产成人无码精品久久久性色| 精品国产91久久久久久久a | 欧美激情精品久久久久| 久久SE精品一区二区| 日韩久久久久中文字幕人妻| 久久精品国产亚洲精品| 久久99热这里只有精品国产| 欧美久久精品一级c片片| 99热精品久久只有精品| 日本免费久久久久久久网站| 久久综合丁香激情久久|