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); }
