[Scoi2015]小凸玩矩阵

http://www.lydsy.com/JudgeOnline/problem.php?id=4443

二分+最大流
 二分一个x，矩阵中的数i若<=x，在这个点和行、列间各连一
条边
 源点向每一行连一条流量为1的边，列向汇点连一条流量为1
的边
 最后判断流量是否>=n-k+1，大于则下调上界， 否则上调下
界
 TLE 优化？
 矩阵中的数i若<=x，就从数i的行向列连一条流量为1的边

#include<cstdio> 
#include<cstring> 
#include<queue> 
#include<algorithm> 
using namespace std; 
int n,m,k,src,decc,l,r,tot=1,ans; 
int a[251][251],b[251000],front[63010],nextt[251000],to[251000],cap[251000],cnt[63010],lev[63010]; 
queue<int>q; 
void add(int u,int v,int w) 
{ 
    to[++tot]=v;cap[tot]=w;nextt[tot]=front[u];front[u]=tot; 
    to[++tot]=u;cap[tot]=0;nextt[tot]=front[v];front[v]=tot; 
} 
bool bfs() 
{ 
    for(int i=0;i<=n+m+1;i++) {cnt[i]=front[i];lev[i]=-1;} 
    while(!q.empty()) q.pop(); 
    q.push(src);lev[src]=0; 
    while(!q.empty()) 
    { 
        int now=q.front();q.pop(); 
        for(int i=front[now];i!=0;i=nextt[i]) 
        { 
            int t=to[i]; 
            if(cap[i]>0&&lev[t]==-1) 
            { 
                q.push(t); 
                lev[t]=lev[now]+1; 
                if(t==decc) return true; 
            } 
        } 
    } 
    return false; 
} 
int dinic(int now,int flow) 
{ 
    if(now==decc) return flow; 
    int delta,rest=0; 
    for(int & i=cnt[now];i!=0;i=nextt[i]) 
    { 
        int t=to[i]; 
        if(lev[t]==lev[now]+1&&cap[i]>0) 
        { 
            delta=dinic(t,min(cap[i],flow-rest)); 
            if(delta) 
            { 
                cap[i]-=delta;cap[i^1]+=delta; 
                rest+=delta;if(rest==flow) break; 
            } 
        } 
    } 
    if(rest!=flow) lev[now]=-1; 
    return rest; 
} 
bool check(int u) 
{ 
    int tmp=0;tot=1; 
    memset(front,0,sizeof(front)); 
    for(int i=1;i<=n;i++)  
    { 
        add(src,i,1); 
        for(int j=1;j<=m;j++)  
        { 
            if(a[i][j]>u) continue; 
            add(i,n+j,1); 
        }  
    } 
    for(int i=1;i<=m;i++) 
     add(n+i,decc,1); 
    while(bfs()) tmp+=dinic(src,63000); 
    if(tmp>=n-k+1) return true; 
    return false;  
} 
int main() 
{ 
    scanf("%d%d%d",&n,&m,&k); 
    for(int i=1;i<=n;i++) 
     for(int j=1;j<=m;j++) 
       { 
          scanf("%d",&a[i][j]); 
          b[(i-1)*m+j]=a[i][j]; 
       } 
    decc=n+m+1; 
    sort(b+1,b+n*m+1); 
    l=1,r=n*m; 
    while(l<=r) 
    { 
        int mid=l+r>>1; 
        if(check(b[mid])) { ans=b[mid];r=mid-1;} 
        else l=mid+1; 
    } 
    printf("%d",ans); 
}