[SCOI2015]小凸玩矩阵-题解

【题目地址-loj】

题意见题面。

我们容易看出，答案肯定是单调的，所以便想到二分，对于判断能否选到，我们发现行和列的限制就是二分图匹配的模型，所以直接套个二分图就好啦。

注意：是第$K$大的数！！，所以判断时为$n-K+1$，因为这个wa了好几遍QAQ

#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;
const int N=300,M=1e6+1,inf=1e9;
int n,m,K,s,t;
int mat[N][N],ls[M],tot;

struct ss{
	int to,last,cap;
	ss(){}
	ss(int a,int b,int c):to(a),last(b),cap(c){}
}g[M<<1];
int head[N+N],cnt;
void clear(){memset(head,0,sizeof(head));cnt=1;}
void add(int a,int b,int c){
	g[++cnt]=ss(b,head[a],c);head[a]=cnt;
	g[++cnt]=ss(a,head[b],0);head[b]=cnt;
}
int dep[N+N],cur[N+N],que[M],p,q;
bool bfs(){
	for(int i=s;i<=t;i++){
		cur[i]=head[i];dep[i]=0;
	}	dep[s]=1;que[p=q=0]=s;
	for(int a;p<=q;p++){
		a=que[p];
		for(int i=head[a];i;i=g[i].last){
			if(!dep[g[i].to]&&g[i].cap>0){
				dep[g[i].to]=dep[a]+1;
				if(g[i].to==t) return 1;
				que[++q]=g[i].to;
			}
		}
	}
	return 0;
}
int dfs(int a,int cap){
	if(a==t||cap<=0) return cap;
	int tot=0;
	for(int &i=cur[a];i;i=g[i].last){
		if(dep[g[i].to]==dep[a]+1&&g[i].cap>0){
			int now=dfs(g[i].to,min(cap,g[i].cap));
			tot+=now;cap-=now;
			g[i].cap-=now;g[i^1].cap+=now;
			if(cap<=0) break;
		}
	}
	if(!tot)dep[a]=0;
	return tot;
}
int maxflow(){
	int ans=0;
	while(bfs())ans+=dfs(s,inf);
	return ans;
}
bool check(int val){
	clear();
	for(int i=1;i<=n;i++)add(s,i,1);
	for(int i=n+1;i<=n+m;i++)add(i,t,1);
	for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)if(mat[i][j]<=val)add(i,j+n,1);
	return maxflow()>=n-K+1;//!!!
}
int main(){
	freopen("matrix.in","r",stdin);
	freopen("matrix.out","w",stdout);
	scanf("%d%d%d",&n,&m,&K);
	s=0;t=n+m+1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			scanf("%d",&mat[i][j]);
			ls[++tot]=mat[i][j];
		}
	}
	sort(ls+1,ls+tot+1);
	tot=unique(ls+1,ls+tot+1)-ls-1;
	int l=1,r=tot,ans=tot;
	while(l<=r){
		int mid=l+r>>1;
		if(check(ls[mid]))r=mid-1,ans=mid;
		else l=mid+1;
	}
	printf("%d\n",ls[ans]);
	fclose(stdin);
	fclose(stdout);
	return 0;
}